CASIMIR EFFECg: IN SYSTEMS IN AND 



OUT OF EquEibrium - Efecto 



Casimir en sis^emas en equilibrio 
y fuera gel equilibrio 



> 

^ By 

o^ 

Pablo Rodriguez Lopez 
o 

memoriaj^ue presenta a la 
Universidad ^^mplutense de Madrid 
para optar al grado^e doctor en clencias fisicas 
Departamento de FfsiCA Aplicada I (Termologia) 

FACULTAD de ClENCIAS FiSICAS 
DiRIGIDA POR EL DR. RICARDO BRITO LOPEZ 



© Pablo Rodriguez Lopez, 2012. 
Typeset in WF^2£. 



Agradecimientos 



No hay palabras para agradecer suficientemente a Ricardo Brito, con el que he com- 
partido los anos que ha Uevado producir esta tesis. Ricardo ha sido un mentor y un 
amigo, Ueno de una paciencia y comprension infinitas. For supuesto, tambien quiero 
agradecer al departamento de Fisica ApUcada I (Termologia) de la UCM (Universidad 
Complutense de Madrid), por la inmejorable acogida y por los esfuerzos hechos para 
que tuviera un espacio donde realizar esta tesis, no hay manera de sentirse mas acogido. 

Lo mismo se puede decir de Thorsten Emig, desde que Uegue a Colonia, ha sido un 
placer trabajar y aprender a su lado. Tambien quiero agradecer al Instituto de Fisica 
Teorica de la Universidad de Colonia, Alcmania (Institut fuer Theoretische Physik, 
Koeln) su hospitalidad durante las distintas estancias que realice alli. 

Quiero agradecer a Rodrigo Soto, Adolfo Grushin y Alberto Cortijo el tiempo de- 
dicado, con ellos es un placer hablar de Fisica, siempre se aprende algo nuevo. 

Agradczco tambien a mis compaiieros del GISC por todo su apoyo, cspecialmente 
a Ricardo, Juan, Luis, Edgar y Jordan, con los que he compartido conferencias, caiias 
y char las, formales e informales, habeis hecho que este camino valga la pena. 

No existen palabras para expresar el agradecimento y cariiio que siento hacia mis 
padres, German y Carmiiia. Gracias a vosotros, a vuestro constante apoyo y confianza 
sin limite, todo esto ha sido posible. Tambien quiero agradecer a mi hermana Maria 
y a su marido, Antoine, los sabios consejos que siempre me han sabido dar, y a mi 
sobrinita Ines, por la alegria que rebosa de la que nos hace a todos participes. 

No puedo evitar acordarme aqui de mis abuelos, Vicente, por desgracias del destino 

iii 



iv 



Agradecimientos 



nunca Uegue a conocerte. Eduardo, Mercedes, Jexuxa, no pasa un dia sin que me 
acuerde de vosotros, y ahora que ya no estais, me doy cuenta de lo mucho que os hecho 
de menos. Espcro que csteis donde quiera que esteis, os podais scntir orguUosos de mi. 
Yo lo estoy de vosotros. 

Por supuesto que me acuerdo del resto de mi familia: mis padrinos, mis primos, 

tios, sobrinos La lista es afortunadamente interminable, muchas gracias a todos 

vosotros. 

Y que decir de mis amigos, Javier, Eduardo y Jacobo especialmente, de la gente de 
la Facultad y de Hypatia en particular. Y de toda la gente del Calasancio (esas clases 
de baile de Alfredo), vale la pena cada minuto que se pasa a vuestro lado. 

Por ultimo, mi mas sincero agradecimiento a ti, Bea. El motivo ultimo de esta tesis 
ha sido la busqueda de la felicidad, y a medida que la he ido escribiendo he aprendido 
que esta felicidad no esta completa si no es contigo. Gracias por acompaiiarme en este 
camino. 

Mientras que estaba escribiendo esta tesis, mi abuela Jexuxa fallecio tras una larga 
enfermedad. Esta tesis se la dedico especialmente a ella. Descanse en Paz. 

Durante el desarroUo de esta tesis he recibido financiacion del programa de becas 
de Formacion de Personal Universitario (FPU), niimero de referenda de la beca: 2006- 
01310 del Ministerio de Educacion y Ciencia (MEG). Dicho programa tambien financio 
una de las estacias en Golonia. Tambien he recibido financion de los proyectos: Grupo 
Interdisciphnar de Sistemas Gomplejos (MODELIGO-GM), Grupo Interdisciphnar de 
Sistemas Gomplejos. Modelizacion y Garacterizacion (Ref: GR58/08), Fluctuaciones 
en sistemas fisicos complejos (Ref: GGG07/UGM/ESP-2870), Fluctuaciones en Sis- 
temas Fisicos de no equilibrio. Aplicacion a Medios Granulares (Ref: PR34/07-15859), 
Modehzacion, Simulacion y Analisis de Sistemas Gomplejos (MOSAIGO) y Orden y 
Fluctuaciones en Sistemas Gomplejos (Ref: FIS2004/271). 



Acknowledgements 



I do not have the words to thank my supervisor Ricardo Brito with whom I have 
shared the long years it took to produce this thesis. Ricardo has been a mentor and 
a friend, always full of patience and comprehension. Of course, I also thank to the 
department of Fisica Aplicada I (Termologia) (Applied Physics I (Thermo logy)) of 
the UCM (Universidad Complutense de Madrid), for the excellent hospitality and the 
efforts made to have a place to carry out this thesis, there is no way to feel most 
welcome. 

The same can be said about Thorsten Emig, since 1 arrived in Cologne, it has been 
a pleasure to work and learn at your side. 1 also thank the Institute of Theoretical 
Physics at the University of Cologne, Germany (Institut fuer Theoretische Physik, 
Koeln) their hospitality during the fellowships that I made there. 

I want to thank to Rodrigo Soto, Adolfo Grushin and Alberto Cortijo the pleasure 
to talk about physics with them, you always learn something new with them. 

I also thank my colleagues of CISC for their support, especially Ricardo, Juan, Luis, 
Edgar and Jordan, with whom I shared conferences, drinks and talks, both formal and 
informal, you have made this journey worthwhile. 

No words can express the acknowledgments and affection 1 feel toward my par- 
ents, German and Carmina. Thanks to them, their continued support and confidence 
without limitation, all this has been possible. I also thank my sister Maria and her 
husband, Antoine, the wise advices that they have always been able to give me, and 
my niece Ines, she always makes us participants of her full joy. 

V 



vi 



Acknowledgements 



I can not avoid thinking of my grandparents, as my grandfather Vicente, whom for 
misfortunes of fate I never got to meet him. Eduardo, Mercedes, Jexuxa, not a day 
goes by without remembering you, and now that you arc not here with us, I rcahze 
how much I miss you. I hope that wherever you are, you can feel proud of me. I am 
proud of you. 

Of course I remember the rest of my family, my godparents, my cousins, uncles, 
nephews, .... Fortunately the list is endless, thank you very much to all of you. 

And what about my friends, Javier, Eduardo and Jacobo especially, the people of 
the Faculty and Hypatia in particular. And all the people of Calasancio (those dance 
classes with Alfredo), every minute spent at your side is worth. 

Finally, my deepest thanks to you, Bea. The reason of this dissertation has been 
the pursuit of happiness, and as I was writing I learned that happiness is incomplete 
without you. Thanks for joining me in this way. 

While I was writing this thesis, my grandmother Jexuxa died after a long illness. 
This thesis is dedicated specifically to her. Rest in Peace. 

During the development of this thesis I was supported by a FPU grant (Becas de 
Formacion de Personal Universitario), the reference number of the scholarship: 2006- 
01310 of the Ministry of Education and Science (MEC). This program also funded one 
of the fellowships in Cologne. I have also received financing from the projects: Grupo 
Interdisciphnar de Sistemas Complejos (MODELICO-CM), Grupo Interdisciphnar de 
Sistemas Complejos. Modelizacion y Caracterizacion (Ref: GR58/08), Fluctuaciones 
en sistemas fisicos complejos (Ref: CCG07/UCM/ESP-2870), Fluctuaciones en Sis- 
temas Fisicos de no equilibrio. Aplicacion a Medios Granulares (Ref: PR34/07-15859), 
Modelizacion, Simulacion y Analisis de Sistemas Complejos (MOSAICO) y Orden y 
Fluctuaciones en Sistemas Complejos (Ref: FIS2004/271). 



Contents 



Title i 

Agradecimientos iii 

Acknowledgements v 

Table of contents xi 

1 Introduction 1 

1.1 Introduction and objectives 1 

1.2 Structure of tlie tliesis dissertation 6 

Part I: Casimir effect in systems subject to noise 11 

2 Dynamical approach to the Casimir effect 13 

2.1 Equilibrium and nonequilibrium fluctuations 16 

2.2 Casimir forces from the average stress tensor 20 

2.3 Computation of Casimir forces 22 

2.4 Reaction-diffusion systems 26 

2.4.1 Uncorrelated noise in time and space 28 

2.4.2 Temporally correlated noise 30 

2.4.3 Maximally spatially correlated noise 32 

2.5 Liquid crystals 32 

vii 



viii 



Contents 



2.5.1 Uncorrelated noise in time and space 34 

2.5.2 Temporally correlated noise 35 

2.5.3 Maximally spatially correlated noise 37 

2.6 Two non hermitian fields system 38 

2.7 Discussion and conclusions 41 

2.8 Appendix 43 

2.8.1 Appendix A: Elizalde zeta function 43 

2.8.2 Appendix B: Derivation of Eq. (2.33) 44 

3 Stochastic quantization and Casimir forces. 47 

3.1 Parisi-Wu formalism 49 

3.2 Casimir forces from the average stress tensor 51 

3.3 Perfect metal pistons of arbitrary cross section 53 

3.3.1 Short distance limit 56 

3.3.2 Large distance limit 57 

3.4 Thermal Casimir force 58 

3.5 Fluctuation of Casimir forces 59 

3.6 Equivalence between dynamical approach to Casimir effect, Stress Ten- 
sor formalism and Partition function formalism 60 

3.7 Conclusions 64 

3.8 Appendices 64 

3.8.1 Appendix A: Derivation of averaged stress tensor over the piston 
surface 64 

3.8.2 Appendix B: Derivation of the fluctuations of the force 67 

4 Extension of Langevin formalism to higher temporal derivatives and 
its application to Casimir effect. 75 

4.1 Langevin equation with higher temporal derivatives 76 

4.2 Green function of the Generalized Langevin Equation 77 

4.3 Correlation function in the steady state 79 

4.4 Casimir force between infinite parallel plates 81 



Contents ix 



4.5 Discussion 83 

Part II: Multiscattering formalism of Casimir effect 85 

5 Introduction to the multiscattering formalism of Casimir effect 87 

5.1 Partition function description of the problem 88 

5.2 Free energy of the system 90 

6 Three-body Casimir effects and non-monotonic forces 95 

6.1 General approach 97 

6.2 Two atoms 99 

6.3 Two macroscopic metallic spheres 103 

6.3.1 Analytical results 103 

6.3.2 Numerical results 104 

6.4 Discussion 106 

7 Pairwise Summation Approximation of Casimir energy 109 

7.1 Introduction 110 

7.2 Diluted limit at zero temperature Ill 

7.2.1 Purely electric and purely magnetic Casimir energy 114 

7.2.2 Coupled electromagnetic Casimir energy 116 

7.2.3 Asymptotic Casimir energy 117 

7.3 Diluted limit at any temperature T 118 

7.3.1 Purely electric and purely magnetic energy 119 

7.3.2 Coupled magnetic - electric Casimir energy terms 120 

7.3.3 Low and zero temperature limit 121 

7.3.4 High temperature and classical limit 121 

7.4 PSA for three bodies system 122 

7.4.1 PSA for general N bodies system 123 

7.5 Second order expansion of diluted limit 126 

7.6 PSA for objects with general magneto-electric coupling 128 

7.7 PSA for Topological Insulators 129 



X 



Contents 



7.7.1 Casimir energy between spheres and plates 136 

7.8 Conclusions 138 

7.9 Appendix A: Obtention of the matricial Green function 139 

7.10 Appendix B: Dielectric polarizability in the diluted limit 142 

8 Casimir energy and entropy in the sphere— sphere Geometry 145 

8.1 Multiscattering formalism of the Casimir energy 147 

8.2 Perfect metal spheres 148 

8.2.1 Casimir energy in the large separation limit 148 

8.2.2 Casimir entropy in the large separation limit 150 

8.2.3 Numerical study at smaller separations 152 

8.2.4 Casimir force in the large separation limit 154 

8.3 Plasma model 156 

8.3.1 Casimir energy in the large separation limit 156 

8.3.2 Casimir entropy in the large separation limit 157 

8.4 Drude model 157 

8.4.1 Casimir energy and entropy in the large separation limit .... 158 

8.5 Thermodynamical consequences 160 

8.6 Conclusions 162 

8.7 Appendix A: Krein formula from multiscattering formalism 164 

9 Casimir energy between non parallel cylinders 169 

9.1 Geometry of the system under study 171 

9.2 Translation matrices between non-parallel cylindrical multipole basis . 173 

9.2.1 Scalar waves 174 

9.2.2 Electromagnetic waves 175 

9.3 Asymptotic Casimir energy 177 

9.3.1 Scalar field 177 

9.3.2 Electromagnetic Casimir energy 179 

9.4 Numerical results 182 

9.5 Conclusions 184 



Contents xi 



9.6 Appendix A: Derivation of translation matrices for scalar cylindrical 
multipoles between non-parallel system of coordinates 187 

9.7 Appendix B: Translation matrices for scalar cylindrical multipoles be- 
tween parallel system of coordinates 190 

9.8 Appendix C: Proximity Force Approximation of Casimir energy .... 191 

10 Conclusion 195 

10.1 Specific outcomes 197 

List of publications 201 



Bibliography 



203 



xii Contents 



"Once you start asking questions, innocence is gone. " 

Mary Astor 



1 

Introduction 



1.1 Introduction and objectives 

Historical introduction 

The field of Casimir effect, or fluctuation induced forces, was started by two seminal 
papers published in 1948. In the first one [1], Casimir and Polder calculated the energy 
between two polarizable neutral atoms in the retarder limit (considering a finite speed 
of light c). In the second one [2], Casimir studied the energy between two perfect 
metal plates facing each other. These calculations were the generalization of London 
energies [3] when a finite velocity of light c is taken into account. In addition, these 
forces were predicted by Van der Waals in the context of the equation of state of gases 
which has his name [4]. Since them, this effect has been studied along the years with 
an exponential growth in the number of papers published. 
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For a historical review of this field of knowledge, see [5] [6] , and see [7] [8] [9] [10] [11] [12] 
[13] [14] as a short list of recent reviews in the field. 

The description of this effect is as follows: Two neutral polarizable atoms in vacuum 
placed at a distance d from each other, attract because the vacuum fluctuations induce 
polarizabilities in them, as a consequence, the induced polarizations interact and an 
attractive force between the neutral atoms appears. When effects due to flnite speed 
of light c are not important, the force is proportional to where h is the reduced 
Planck constant, and ujr is the frequency of resonance of the atoms. This is the result 
obtained by London [3]. When the effect of flnite c becomes important, we are in 
the so called retarded limit, where the decay of the force with the distance is more 
pronounced (the force decays with instead (i"''), and is proportional to hcaia2-, 
where is the electric polarizability of the i— atom. As this force depends on the 
induced polarizations, it is expected a dependence with the temperature. In fact, 
when At = tends to zero (it is the so called classical limit, or high temperature 
regime), the force is proportional to ksT and reduces the decay to rf^'' again. 

Casimir energy had a theoretical importance in the early times of Quantum Field 
Theory (QFT) as a consequence of the reality of the zero point energy [7]. 

Lifshitz extended the Casimir calculus of perfect metal plates at zero temperature 
to the more general case of dielectric plates at any flnite temperature [15], which is 
one of the most famous and more used formulas of Casimir effect even nowadays. In 
contrast to the description of Casimir effect as a consequence of the zero point energy, 
Lifshitz proposed that this effect comes from fluctuating sources of the electromagnetic 
fleld inside the interacting objects. 

There were another important contributions in the fleld, the result of Boyer [16] 
of the Casimir energy of an spherical shell has also an historical importance, as the 
derivation of Balian and Duplantier [17] [18], just to cite a couple of them. 

In Ref. [19], the multiscattering formalism was derived for the flrst time. It was 
later used by Kardar [20], and nowadays it leads to its modern form [21] [22] [23] [24] [25], 
which accounts with dielectrics of any geometry and nature at any given temperature 
T. The multiscattering formalism will be briefly introduced in Chap. 5, because it will 
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be widely used in the part II of this Thesis. 

In the second part of this Thesis, we center our study in the multiscattering for- 
mahsm of the EM Casimir effect. 



Numerical and Analytical Approximations to Casimir Effect 

Several approaches have been proposed to obtain reliable results of Casimir effect. The 
Pairwise Summation Approximation (PWS in the literature, PSA here) [15] [26] [27], 
which consists of the assumption of a pairwise behavior of Casimir energy. Then, the 
Casimir energy between two objects is obtained as the sum of the relative Casimir 
energy between their constituents. It can be demonstrated that PSA is a good ap- 
proximation in the diluted limit [28], but it gives spurious infinities when it is applied 
to metals. The Proximity Force Approximation (PFA) [29] [30] [31] [32], or Derjaguin 
Approximation, consists of assuming that the surface of each arbitrary body consists 
of infinitesimal plates. Then the Casimir energy between the objects is the surface 
integral of these infinitesimal Casimir results. It is a valid approximation when objects 
are close together and their curvature radius are small. The great disadvantage of 
PFA is that it is an uncontrolled approximation. It is not clear at a quick glance the 
orientation of the infinitesimal plates, there are several criteria which lead to different 
results [33]: (i) Plates perpendicular to the union axis between the center of the two 
objects; (ii) Minimum distance between the points; or (iii) plates perpendicular to the 
union axis between the center of one object with the nearest point of the surface of 
the other. Semiclassical approximations to Casimir effect lead to PFA results [30] , and 
Optical Paths approximation [34] [35] lead also to results of the same nature. From the 
stress-tensor formalism of Casimir effect, several numerical methods have been pro- 
posed [36], but the modern formulation of Multiscattering formalism has led to several 
new methods [37] [33]. For a modern review in numerical calculations of Casimir effect, 
see [13]. 
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Experimental Results 

Early experiments carried out in the late 50's provided at best a qualitative support 
for an attractive force [38], but they were quantitative inconclusive, because of an 
experimental error of 100%. The difficulty of the experiments has multiple causes. 
The first of all is that Casimir force is weak, and it has a microscopic nature (in 
general), because its proportionality to he. The second one is that the experimentalists 
have to take into account a plethora of corrections associated to the phenomenon, as 
the rugosity of materials, their dielectric properties, finite temperature effects, . . . and, 
finally, there is the difficulty to obtain the desired geometry of two microscopic parallel 
plates near to contact, but without being in touch. 

In 1997, Lamoreaux [39] with a torsion pendulum, and in 1998, U. Mohideen and 
A. Roy [40] measured the Casimir force between a plate and a sphere with a cantilever. 
Several experiments followed to that one and started a new era of interest in Casimir 
effect experiments and theory p. e. [41]. 

Several more experiments have been performed in agreement with the theory: 
Casimir levitation [42], Casimir effect in critical systems [43], . . . 

Ubiquity of Casimir Effect 

Even from the early times, Casimir effect has been considered a consequence of the 
quantum fluctuations of the EM fleld. Casimir himself claimed that, speaking with 
Bohr about this phenomena, Bohr claimed that this effect must have something to do 
with the zero-point energy [5]. 

In fact, this effect has always been related with fluctuating fields and with the Van 
der Waals force, proposed by Van der Waals himself to justify the existence of Phase 
transitions in classical Statistical Physics [44]. 

After that, it has been proposed the existence of Casimir force between objects 
immersed in fluctuating flelds as: QED, critical systems [45], systems in their ther- 
modynamical equilibrium [46], fluids [47], granular media [48], out of equihbrium sys- 
tems [49] [50] [51], avalanches [52], capilar Casimir effect [53] [54], colloids [43], . . . 
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All these systems share a property, there are fluctuating fields which induce an 
effective force between intrusions. 

In the first part of this Thesis, we propose a dynamical formulation of Casimir effect 
from Langevin PDE's which models the fluctuations of the background fleld. This 
approach is not new in the sense that it has been used previously to model the Casimir 
effect in systems in thermodynamic equilibrium, as the case of liquid crystals [46]. Our 
formulation allows us to identify what these systems of different nature have in common, 
which is the existence of fluctuations, and thus it make explicit the relationship of the 
Casimir force with them. In fact, it allows to deflne the Casimir force as the response 
of a fluctuanting medium to the breakdown of the translation symmetry because of the 
presence of intrusions in that medium. 

Related Physical Phenomena 

In this Thesis, our study is centered in the Casimir effect which appears between static 
intrusions in a surrounding field, which is placed in a steady state induced by a noise. 

This is equivalent to assume quasi-static systems in Thermodynamics, but it does 
not have to be true when the intrusions move, because the assumption that we have 
placed in a steady state is broken. The simple non-free motion of a particle in a 
medium can break the status of thermodynamic equilibrium. 

Therefore, we do not study the case when the intrusions have a non zero velo- 
city [55], when the steady state changes with time [56], or about other possible out of 
equilibrium situations [57]. 

Casimir effect is not the only physical phenomena associated to fiuctuations. 

In fact, there exists the Dynamical Casimir effect [58] [59], which consists of the 
generation of thermal photons of accelerated intrusions (recently measured in [60]), it 
is also proposed a vacuum friction of intrusions at constant velocity from each other [55] . 

There also exists the Brownian motion of intrusions, also originated by fluctuations 
of the surrounding fleld, but independent of the Casimir effect. In fact, we could say 
that, if Brownian particles interact, it is due to the Casimir effect, because the fluctu- 
ations of the surrounded fleld can generate the erratic motion of the Brownian particle 
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at the same time that the constriction of such fluctuations induce a Casimir force be- 
tween the particles. It also appears induced conductivity because of the proximity of 
objects [57], and a plethora of other phenomena associated to fluctuations, as the prob- 
lem of Dark energy, constrictions to modifications of Newton's Law of gravitation. . . 

As we will see in the second part of this Thesis, Casimir effect is a generalization 
of Statistical Physics when we include the effect of boundaries in our system. 

1.2 Structure of the thesis dissertation 

The research presented here is structured in two blocks that conform the two parts of 
this Thesis: 

• Part I: The dynamic formalism of Casimir effect is presented and studied here. 
This model complements the already known models in equilibrium systems, but 
this let us extend the study of Casimir effect to non-equilibrium systems. Then it 
allows us to define the Casimir force as the response of a fluctuanting medium to 
the breakdown of the translation symmetry because of the presence of intrusions 
in that medium. 

— In Chapter 2, we present the dynamical formalism of Casimir forces between 
intrusions originated from the fluctuations of the surrounding media, which 
is necessary to extend the study of Casimir effect to non-equilibrium sys- 
tems. In fact, it is explained why, in general non equilibrium steady states, 
we need such dynamical formalism and we define the stochastic variable 
Stochastic Casimir force over a given a body, which let us, by the use of 
stochastic calculus, define the mean value of this stochastic variable as the 
Casimir force itself. As an advantage, such formalism avoids the appear- 
ance of spurious divergences in the theory, because quantities independent 
of the distance between objects are automatically removed. As a particular 
example, we study the Casimir force between parallel plates immersed in 
a reaction-difussion media, and the case of parallel plates immersed in a 
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nematic liquid crystal with an attenuation term because of the presence of 
an external magnetic field, in equilibrium and out of equilibrium cases. We 
also study the Casimir force between parallel plates immersed in a system 
whose fluctuations follow a non Hermitian evolution, something forbidden in 
equilibrium systems. The theory recovers the classical equilibrium Casimir 
force as a particular case, and admits generalizations to non-linear dynamics 
of the background (critical phenomena, Ginzburg-Landau, KPZ, . . . ) and 
to more general noises (multiplicative noises. Levy noises, . . . ). 

— In Chapter 3, by the use of the Stochastic formulation of QFT, also called 
Parisi-Wu formalism, we extend the dynamic formalism of Casimir effect to 
the important case of electromagnetic field subject to quantum fluctuations. 
As an example, we use this formalism to obtain the Casimir force between 
two perfect metal pistons of arbitrary given section at any temperature, 
and the limits of short and long distances. These results are consistent with 
results in certain limits of the literature. The use of the dynamical formalism 
also let us to obtain the variance of the Casimir force between pistons at 
a any given temperature. This result is more difficult to obtain from other 
formulations of the theory of Casimir effect. As a conclusion, this Chapter 
show us how Casimir effect appears in QED because of quantum fluctuations 
of EM field on an explicit form, and it is nothing else but a particular 
case of fiuctuation-induced forces. In addition to that, we obtained the 
explicit equivalence between the Stress-Tensor formalism of Casimir effect, 
the partition function formalism, and the dynamic formalism of Casimir 
effect in the particular case of equilibrium presented in this Thesis. 

— In Chapter 4, we present another generalization of the dynamical formal- 
ism. Already in Chapter 2 it was shown how to generalize the dynamical 
formalism to non-linear fields and to other kind of more general noises, but 
the generalization presented here is in the temporal evolution of the fiuctu- 
ations. We want to know what happens if a steady state is reached from a 
dynamic whose order of the temporal dynamic is greater than one and how 
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would be modified the Casimir force in sucli non equilibrium steady state. 
As a result, we obtained the two point correlation function of the field in 
such steady states, and therefore we were able to obtain the Casimir force 
between parallel plates in those systems. It is important to note that this 
is not a useless result, when fluctuations are of an external origin, they can 
affect to fields whose temporal evolution was described by temporal deriva- 
tives of order greater than one, as it is the case of classical Electromagnetic 
field subject to random currents described by a white noise. Such case is 
qualitative different of the Parisi-Wu formalism studied in Chapter 3. 

• Part II: The multiscattering formalism of EM Casimir effect is studied here in 
order to obtain results with a physical importance. 

— In Chapter 5, we give a short introduction to the multiscattering formalism, 
required to understand how this formalism is derived and how we use it in 
the following Chapters. 

— In Chapter 6, we study the superposition principle in Casimir effect and the 
appearance of nonmonotonicities in the Casimir force between two objects 
because of the presence of a third one in the system. Is it the Casimir force 
over an object equal to the sum of pairwise Casimir force of this object 
with the rest — 1 objects of the system?. The answer is not. Casimir 
force does not follow a superposition principle, neither to usual EM for 
example. In particular, with a modified multiscattering formalism (which 
takes into account the images of the objects), we obtained how the Casimir 
force between two neutral atoms and between two perfect metal spheres 
varies when a perfect metal plate is placed in the system. As a result, the 
Casimir force between compact objects, and between a compact object and 
the plate not only depends on the distance with the third object, but also 
it is nonmonotonous with such distance. Sometimes, force increases with 
the distance (but it do not reverse sign) because of the presence of a third 
object in the system. 
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— In Chapter 7, we derive the Pairwise Summation Approximation (PSA here, 
PWS in the hterature) for the EM field from the muhiscattering formahsm. 
PSA is the tree level of a Born expansion of the Casimir effect in the dielec- 
tric properties of the objects. This approach allows us to reobtain several 
results of the theory, for example the limits of high and low temperature of 
the Casimir energy between two objects. But we also obtain new results, as 
the energy of two objects at any given temperature and the magnetic con- 
tribution to PSA. We have also studied the superposition principle, which is 
obtained at the tree level of PSA, but broken by its first perturbation. From 
this formalism we do not only reobtain the point potential formalism of 
Power and Thirunamachandran, but we also are able to extend such formal- 
ism to magnetic couplings and finite temperature cases. Last, but not least, 
we study the Casimir effect between dielectric with magneto-electric cou- 
pling in general, and between Topological Insulators (TI) in special. These 
new state of matter materials leads to attractions, repulsions and even to 
the appearance of equilibrium distance points. By the use of PSA, we obtain 
a good approximation of the Casimir energy between objects with a general 
magnetoelectric coupling, in particular for TI, and we derive a criteria for 
the appearance of all the behaviors described here. 

— In Chapter 8, we study the entropy for a system of two metal spheres 
with three different models: perfect metal model, plasma model and Drude 
model. This work is motivated by the appearance of negative entropies 
because of Casimir effect between two Drude plates and between a sphere 
and a plate, both perfect metal. As a result, we obtained that an interval 
of negative entropies appears between perfect metal spheres and for some 
parameters of the plasma model, for a range of distances and temperatures. 
In the case of perfect metal spheres, there exists a minimum critical distance 
from which this anomalous behavior disappears. For Drude spheres and for 
plasma model with large enough penetration length this effect simply does 
not appear. A thermodynamical study of this system is performed and, as 
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a consequence, we claim that the only consequence of this range of negative 
entropies is a nonmonoticity of the Casimir force with the temperature. 

— In Chapter 9, we obtain the Casimir energy for a system not studied before, 
it is the case of non-parallel cylinders. To do so, we obtain the translation 
matrix of non parallel and translated cylindrical coordinates basis, for the 
scalar and EM fields. In this system, the energy does not longer scale with 
the length L of the cylinders, but with the cosecant of the angle 6 between 
their axis. For perfect metal cylinders in the high T limit, we find a case 
when the theorem of existence of a trace-class operator generator of Casimir 
energy is not applicable and is not fulfilled [61]. Anyway, we are able to avoid 
this anomalous result evaluating the force instead of the energy in this case. 

The main results of this Thesis are summarized in the conclusions Section. 
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"The same equations have the same solutions." 

Richard P. Feynman 



2 

Dynamical approach to the Casimir effect 



In this Chapter we present a mathematical model for the evaluation of fluctuation 
induced forces which appears in objects immersed in fluctuating media. We discuss 
the interest and the generality of these systems. They arc ubiquitous in nature, as well 
as previous work in the literature. This will lead to understand why we call to the forces 
acting between bodies mediated by a fluctuating medium as Casimir forces. We will see 
that in the case where we have an environment in thermodynamic equilibrium, we could 
use all the tools of Statistical Physics and Thermodynamics to study this phenomenon, 
so partition functions, free energies and thermodynamic formalism will be used, whereas 
if we consider a medium in a steady state but outside the thermodynamic equilibrium, 
thermodynamics potentials are meaningless. The partition function, characteristic of 
the problem, may be useless or even unsolvable, so we only have the option of a dynamic 
system formalism. 

13 
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There are many systems in nature which are subjected to fluctuations, of thermal 
or quantum origin. For such systems, under certain physical conditions, Casimir forces, 
created by the confinement of fluctuations, exist and have been calculated (see, e.g., 
[12]). The usual way to obtain the Casimir forces uses equilibrium techniques and is 
therefore valid only for systems in thermodynamic equilibrium. This means that the 
fluctuations must satisfy a fluctuation-dissipation theorem that guarantees the exis- 
tence of an equilibrium state, as discussed below. Casimir forces for these systems are 
calculated in the spirit of the original work of H. G. Casimir for the electromagnetic 
case [2]. The method takes as a starting point the Hamiltonian of the system, from 
which the partition function Z = j exp(— is calculated, either directly or using 
functional integration [21]. In the calculation of the partition function one must take 
into account the boundary conditions, that is, the macroscopic bodies which are im- 
mersed in the system. The partition function of the system will have different values 
for different configurations, e.g., different separations of the objects. Once the partition 
function has been obtained, its logarithm provides the free energy F. The final step 
required to obtain the Casimir force is the calculation of the pressure as the difference 
in the free energy when the configurations of the macroscopic bodies change (for ex- 
ample, changing their position, distance or sizes). For instance, in the usual Casimir 
case of forces between two flat parallel plates at separation L, the force per unit area 
is given by Fc/A = —dF/dL. 

The second approach also takes as a starting point the Hamiltonian of the sys- 
tem. However, in this approach the Casimir force is derived not from the free energy 
but from the stress tensor T, which is integrated over the surface of the macroscopic 
bodies and then averaged over the thermal Boltzmann distribution of the associated 
Hamiltonian exp(— /3if). The approach based on the stress tensor has been taken by 
several authors [62] [49] [63] [64] . In fact, both approaches are equivalent and valid for 
equilibrium systems only, as will be demonstrated in Sect. 3.6. The reason is that both 
are based on properties which are only valid in equilibrium situations. The former uses 
the thermodynamic relation for the pressure as the derivative of the free energy with 
respect to the volume, and the latter uses the Boltzmann distribution function, which 
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is only valid for systems in equilibrium. 

On the other hand, other authors have developed a dynamical approach [46] [49] [65] 
[56]. Here the starting point is an evolution equation for the considered field(s), supple- 
mented with a noise source term, so that the evolution of the field takes the form of a 
Langevin equation. Once this equation is solved, the field is inserted into the expression 
for the pressure and the average over the noise is taken. As we will see in the next Sec- 
tion, if the noise is of internal origin, say thermal or quantum, this description reduces 
to the equilibrium one, because of the fluctuation-dissipation theorem. However, the 
noise does not necessarily have to be internal but can have an external origin [66], for 
instance, a system in a fluctuating temperature gradient [67], subjected to external en- 
ergy injection such as vibration [68] or electrically driven convection [69], light incident 
on a photosensitive medium [70], or spatially and/or temporally correlated noise, as 
considered in Ref. [49]. Recently, Ref. [71] generalized this method to a nonequilibrium 
temperature gradient. In none of these cases can the equilibrium approach be applied. 
Also, the internal dynamics cannot satisfy the condition of detailed balance, and there- 
fore the internal noise is not described by the fluctuation-dissipation relation. In both 
cases, it is only possible to calculate Casimir forces via the dynamical approach. In all 
these nonequilibrium cases a common feature shared with the equilibrium Casimir force 
is that the origin is the limitation of the fluctuation spectrum at large wavelengths. 
They are, therefore, conceptually different from other fluctuation-induced phenomena 
such as ratchets or Brownian motors that act at small length scales. 

The plan of this Chapter is as follows. We start in Sect. 2.1 by presenting the 
Langevin equation subjected to a general noise. We stress the differences between 
the cases when the Langevin equation derives from an energy functional or not, and 
discuss the implications of the fluctuation-dissipation theorem. Section 2.2 derives the 
Casimir force from the stress tensor, while Sect. 2.3 calculates the actual Casimir force 
by substituting the solution of the Langevin equation into the stress tensor. 

The subsequent Sections 2.4, 2.5, and 2.6 are devoted to the application of the 
formalism to different physical systems and different nonequilibrium conditions, that is, 
different ways of violating the fluctuation-dissipation theorem. In particular. Sect. 2.4 
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studies a reaction-diffusion system with three types of noise: (1) a noise uncorrelated 
in space and time, (2) a noise exponentially correlated in time, and (3) a spatially 
homogeneous noise, only fluctuating in time. Section 2.5 is devoted to the study of a 
liquid crystal, with an equilibrium noise, satisfying the fluctuation-dissipation theorem 
and therefore in an equilibrium situation. We continue with a temporally correlated 
noise and finish the Section with a maximally correlated noise. Finally, to illustrate the 
power of the method, we apply it to a two-field system where the evolution equation is 
non-Hermitian in Sect. 2.6. Usual approaches, based on equilibrium properties, have 
no applicability in this case. We finish with some conclusions. 

The contents of this Chapter is based on the work published in [72]. 

2.1 Equilibrium and nonequilibrium fluctuations 

The most widely used tool to study the dynamics of fluctuations is the Langevin equa- 
tions and its related Fokker-Planck equation. There is a wide literature on this subject, 
in particular using Langevin equations; see, for example, Refs. [73] [74] [75] [76] [77] . 
Let us consider a linear stochastic differential equation for the field 0(r,t), 

9t0=-M0 + e(r,t), (2.1) 

which is a generalization of the Langevin equation to spatially extended systems. In 
this equation, M is an operator (usually differential) that can be Hermitian or non- 
Hermitian. The operator does not depend on the field (p, so the Langevin equation 
(2.1) is linear. To simplify notation, we have assumed Langevin equations without 
memory, but the generalization to memory kernels is direct. 

For the important case of critical fluctuations, which leads to the so called Critical 
Casimir effect, the correct associated Langevin equation is a Ginzburg-Landau equation 
with a white additive noise term, then the hypothesis of linearity of Ai operator fails in 
this case. It has consequences in the study of critical Casimir effect with this formalism, 
but the linearity of Ai is enough general to restrict us to this hypothesis, as will be 
shown. 
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The term ^(r,t) is a Gaussian noise that represents the random or stochastic force 
acting over the field 0, and therefore it is the source of fiuctuations for 0. It is customary 
to assume that the noise is Gaussian, and its averages are 

(^(r,t)) = 0, (2.2) 
(e(r,t)e(r',tO) = Q6ir-r')5it-t') = hir-r')6it-t'), 

where Q is a Hermitian operator that can contain differential and integral terms. Dif- 
ferential terms characterize noise of conserved quantities, and integral terms noises 
with spatial correlations. The application of this operator to the Dirac delta function 
produces the spatial correlation distribution h. The noise here is uncorrelated in time, 
although temporal correlations will also be considered below. 

Equation (2.1) admits a solution for an initial condition 0°(r) as 

0(r, t) = e-^V°(r) + e"^* T dr e^^^{r, r). (2.3) 

Jo 

In the limit t oo, <f){r,t) reaches a stationary state iff e~"^*</)° — t- . This implies 
that the eigenvalues of A4 must have positive real parts. 

From the Langevin equation, one can construct a functional Fokker-Planck equation 
for the probability distribution P of the field 0. The technique is standard (see, e.g., 
[74]), and its solution (which is not normalizable) is a Gaussian of the form 

P[0] = — ^e-^/*^'^^, (2.4) 



where the Hermitian operator /C is the solution of the functional Lyapunov equation [77] 

MIC-^ + 1C-^M+ = Q, (2.5) 

where denotes the adjoint of A^. The probability distribution P[0] depends both 
on the matrix Ai and also on the intensity of the fiuctuations Q via Eq. (2.5). 

Unfortunately, although it is known that the Lyapunov equation (2.5) is solvable 
for finite dimensional systems [78], to the infinite dimensional case is unknown by us 
whether a general solution of equation (2.5) exists, or even if can be solved analytically. 
While this question were without an answer, we will not be able to solve the equation 
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(2.5) analytically, so that we are not capable to obtain an analytical closed form for 
the /C operator as a function of M and Q. 

The partition function, or the normalization constant of the probability distribution 
P[4>], is defined as the integral over the whole space of definition of our random variable, 
which in the case of the field can be written as the functional integral 



For all the above, we have a situation in which, for steady states in general, we know 
there exists a well defined partition function, but despite knowing that this is solvable 
as a Gaussian functional integral, it is not useful to us, because in general we do not 
know the functional /C. 

What happens now if the system is at equilibrium? For this case there exists an 
energy functional F (that can be either the entropy [73], a Lyapunov functional [79], 
the Hamiltonian, or a free energy [80],- ■ ■ ) which is an integral over space of a local 
functional J-", which depends on the field and its gradients. The evolution equation 
for can be obtained by generalization to the continuum of the Thermodynamics 
of irreversible process (see, e.g.. Chapter Vll of [73]). This theory relates the time 
evolution of the fields with its conjugated variables $, or the so-called thermodynamic 
forces, as 



Here C is the (symmetric) Onsager operator, which is a generalization to continuum 
of the Onsager matrix of transport coefficients. It is also called the dissipation matrix. 
The second law of Thermodynamics requires that the real parts of the Onsager operator 
eigenvalues are positive in order to guarantee local increase of entropy. The fields $ 
appearing in (2.7) are the conjugated variables of and can be derived from the 
functional F as 



If the fiuctuations are small and the system is far from a phase transition, we can 
assume that 6F/6(j) is linear in the field 0. This implies that F[0] is bilinear in the field 




(2.6) 



9i0 = -£$ + e(r,t). 




5£ 
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0: 

F[<P] = j dv ^(0, V(P,...) = ^j dr ^g^, (2.9) 

and the equation above defines the operator Q. Although the definition of J-' may not 
be unique, Q is unique. For instance the two local energy functions J^i = |V0p/2 and 
J-2 = —(f)A(f)/2 yield the same Q operator: Q = A. In addition to that, Q must be 
positive definitive in order to guarantee the existence of a minimum of the free energy. 
It can also be chosen to be Hermitian, because the antisymmetric part (if there is 
any) does not contribute to the free energy. In this way, the corresponding Langevin 
equation far away of a Phase Transition is linear in 0. Therefore, we can write that 




(2.10) 



Combination of Eqs. (2.7) and (2.10) allows us to write the evolution equation for 
(j) as 

dt<P = -Cg<P + ar,t). (2.11) 

If the system is at equilibrium, the well-known Fluctuation-Dissipation Theorem [81] 
imposes that the intensity of the noise, given by Q, must be related to the Onsager 
operator £ by 

Q = kBT{C + C+). (2.12) 

Equation (2.11) is formally equal to Eq. (2.1), with the operator Ai given by 
Ai = CQ. As both Q and {C + £^) are Hermitian and definitive positive, it can 
be shown that the eigenvalues of Ai have positive real parts, even though Ai can 
be non-Hermitian or undefined (as in the case of the linear hydrodynamic equations) 
(see [73], Chapter V) [82]. In both equilibrium and nonequilibrium dynamics we will 
assume that the real part of the spectrum of Ai is strictly positive, i.e., there are no 
neutral modes as happens when there is continuous symmetry breaking [83] or critical 
phenomena [84] [85]. 

In equilibrium, the fiuctuation-dissipation relation has drastic consequences for the 
solution of the Langevin equation associated to Eq. (2.11). The equation (2.5) is now 
written 

£gic-^ + ic-^gc+ = kBT{c + c+), (2.13) 
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which admits the solution /C = PQ. Once substituted into Eq. (2.4), the probabihty 
distribution is given by the exponential of the functional F multiplied by (3 = {ksT)'^. 
More precisely, 

P[0] = ie-^^['^], (2.14) 

where Z is the partition function or the normalization constant of P[4>]- Given this 
probability distribution P[4>], we can now calculate the average of any dynamical vari- 
able A{(f)) as 

{A) = J P</./l(</.)P[0]. (2.15) 

In particular, the average of the functional F can be calculated as 

, , d\n(Z) 

(F) = (2-16) 

Equations (2.14) and (2.16) are only valid for equilibrium systems, for which an energy 
functional exists and the fluctuation-dissipation theorem is valid. However, if the 
system is out of equilibrium, the probability distribution is not the exponential of F 
and therefore its average is not given in terms of the partition function Z. 



2.2 Casimir forces from the average stress tensor 

How is this discussion related to the calculation of Casimir forces? The Casimir force 
is normally calculated for equilibrium situations, that is, when the noise is of thermal 
origin and the fluctuation-dissipation theorem is satisfied. One way to calculate the 
Casimir force is by the evaluation of the stress tensor Tij. From the functional F, the 
stress tensor is calculated as [86] 

J.. = _ Wid)^^ - 2V ' 



T = r[0,0,r], (2.17) 



which allows the definition of the symmetric bilinear stress tensor operator f. For 
isotropic systems, the local stress is simply given by the diagonal components of the 
stress tensor, or by one-third of its trace. Usual forms of T are A0(r)^ times the 
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identity matrix or a tensorial product of gradient as in liquid crystals, but it can be 
also nonlocal, as in [87]. 

Because of the intrinsically fluctuating nature of the flelds, the stress tensor has to 
be averaged over the random flelds ^(r,t) or the probability distribution (2.4). Once 
we have the averaged stress tensor, the Casimir force over a body a of surface is 
obtained as the sum over all the points of the surface of the local pressure over each 
point because the external fluctuanting medium, it is 



where the integral extends over the surface of the embedded bodies and the vector n 
is a unit vector normal to the surface, pointing inward the body. 

As in the original Casimir calculation, the geometry that will be considered through- 
out this paper consists of two parallel, inflnite plates, perpendicular to the x-axis, sep- 
arated by distance (Fig. 3.2). In this geometry, the Casimir force per unit area on 
the plates is then the difference between the normal stress on the interior and exterior 
side, where the latter is obtained by taking the limit L'^ — > oo. The force per unit area 
on the left plate is 



The interpretation is that, if FcjA is negative, the plates repel each other, while if 
it is positive, an attraction between the plates appears. Note that this is not the 
conventional interpretation of signs. 

Let us discuss Eq. (2.19) for equilibrium and nonequilibrium situations. In the 
former case, i.e., a system in equilibrium, the average of the stress tensor can be taken 
in two ways: as an average over the probability distribution given by Eq. (2.14), or as 
an average over the fluctuating term ,^(r,t). Equilibrium Thermodynamics guarantees 
that both averages are the same. In contrast, in a system out of equilibrium, we are 




(2.18) 



In components it is written as 




(2.19) 




(2.20) 
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Figure 2.1: Parallel-plate geometry used to compute the Casimir force. The system is 
confined between plates located at x = and x = L^- Additional plates are located at 
distances L'^ from these plates, and finally the limit — )• oo is taken to mimic an infinite 
system. 

left with one option, the average over the noise ^(r, t), because K. cannot be obtained 
in general. As mentioned, the system can be out of equilibrium if the fluctuation- 
dissipation relation is not satisfied. In this case, there still exists a functional F (from 
which the Langevin equation is constructed), and the stress tensor can be defined via 
Eq. (2.17). Then, the average in (2.19) has to be taken over the noise. Finally, a 
more complex situation is when the Langevin equation is in its most general form, i.e., 
Eq. (2.1), without Ai deriving from an Onsager matrix and a functional F. In this 
case, the stress tensor cannot be constructed from Eq. (2.17), and one must appeal to 
other considerations in order to construct a stress tensor. One can use a microscopic 
analysis of momentum transfer, kinetic theory, or invoke, for instance, the existence of 
a hydrostatic pressure from which the Casimir force can be derived. We will assume, 
therefore, that it will always be possible to build the stress tensor operator X and, 
then, compute the Casimir force. 

2.3 Computation of Casimir forces 

In this Section we will develop a formalism, valid for both equilibrium and nonequi- 
librium systems, that allows us to compute the average stress tensor and therefore 
the Casimir force. We will assume that the dynamics close to the stationary state 
is described by the dynamical equation (2.1), where the noise term is assumed to be 
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Gaussian with vanishing mean. We assume that the noise has temporal and spatial 
correlations, but no cross-correlations, 

(e(r,t)r(r',t')) = Mr-r')c(t-t')- (2.21) 

Note that we have assumed a dynamical model whose deterministic part is local in time 
(no memory) but that the noise can have some memory. This possibility is not allowed 
by the fluctuation-dissipation theorem, and therefore the system is automatically put 
out of equilibrium. A necessary condition to recover equilibrium is a local correlation 
in time, although this condition is not sufficient, as shown in Sect. 2.1. 

To solve (2.1) we construct the left and right eigenvalue problems of Ai with the 
appropriate boundary conditions over the immersed bodies. Although we will consider 
the case of two parallel plates, the formalism developed in this Section is completely 
general. The left and right eigenvalue problems read 

Mfn{r) = /i„/n(r), (2.22) 
M+gn{r) = filgnir), (2.23) 

with the boundary conditions provided by Ai (which are the same as those of C if 
the dynamics derives from a free energy functional). The left and right eigenfunctions 
are orthogonal under the scalar product, i.e., {gn\fm) = / drg^{r)f,n{r)] that is, under 
appropriate normalization, {gn\fm) = ^nm- We can project the field and the noise over 
the left eigenvalues 

0(r,t) = 5^0„(t)/„(r), e(r,t) = 5^e«(t)/n(r), (2.24) 

n n 

where 0„(t) = {gn\(p(t)) and ^„(t) = {gn\^(t)). By inserting these expressions (2.24) 
into the evolution equation (2.1) we get the evolution equation of each mode 0„(t) as 

dtMt) = -l^nMt) + Ut)- (2.25) 
In order to solve this equation, we obtain its correspondent Green function 



idt + fin)Git,t') = 6it-t'), 



(2.26) 
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which is 

G{t, f) = e-(*-*')^"e(t - f), (2.27) 
then the analytical solution of (2.25) is 



oo 

Mt) = I dt'Git,t')Ut'). (2.28) 

^— oo 

Mt)= f rfre-(*-^)'^"en(r) = e-*^"0„(O) + f dre'^'-^^^-Ur)- (2.29) 

J -oo Jo 

The first term e~^"*(/)„,(0) is a transient term that vanishes for times longer than t ^ 
^^^-^ y so that the average of each mode over the noise ^ is zero in this limit. 

To compute the average stress tensor at each point, we need to compute {'T'[(f), 0, r]). 
Expanding on the eigenvalue basis we get that, in the steady state, 

lim(T(r,t)) = hm V(0„(t)0;,(t))r„„(r), (2.30) 

t— >oo t— >oo ' 

m,n 

where T„m(r) = T[fn, f^,r]. 

The cross-average of the mode amplitudes is obtained from (2.29) and in the sta- 
tionary regime [t ^ l/Re(/i„), l/Re(/im)] can be written as 

lim (MtWmit)) = lim e-('^"+'^^)* /* dn f dr^e^"^^^^'-^^ {Uri)Cm{r2)) • (2.31) 

J — oo J —oo 

Therefore, we need to calculate the correlation of the n and m components of the noise 

{Uri)C^{r^)) = jdr, j dv^gl{v^)g^^{v^) (^(n, ri)r(r2, Ts)) . (2.32) 

Substituting Eq. (2.21) into Eq. (2.32) and the result into Eq. (2.31), we obtain (see 
Appendix 2.8.2), 

lim {UtWmit)) = hn J^^'-^^Y"'^ ^ (2-33) 



where 



hnm = j dvi j dv2gliri)h{ri - r2)fi'„(r2) = {gn\Qgm) (2.34) 



and c is the Laplace transform of c. Finally, the local average of the stress tensor in the 
stationary regime, where transients have been eliminated and the value of the stress 
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tensor is independent of time, is given by 



(T(r)) = J2 



nm 




(2.35) 



This expression is generally divergent when summed over all eigenfunctions. This di- 
vergence comes from the highest eigenvalues (corresponding to small wavelengths) and 
is due to consider the mesoscopic dynamics given by Eq. (2.1), valid for all wavelengths. 
However, it is only valid above a certain minimal distance (the atomic or molecular 
length, for example). There are some techniques to avoid this divergence. For instance, 
a short- wavelength cutoff could be introduced as in Ref . [50] , but here we will use regu- 
larization techniques similar to the Riemann zeta function used in the electrodynamic 
case [2]. 

Using the previous expression, the conditions under which Casimir forces exist 
in an equilibrium system can be deduced. As mentioned above, if the dynamics is 
local in time, the fluctuation-dissipation theorem implies that the noise terms must 
not have memory either, then c(t — t') = 6{t — t') and therefore c(/i) = 1/2. Also, 
the equilibrium relation between noise autocorrelation and Onsager matrices given in 
Eq. (2.12) implies that h^m = f^BT{fin + fJ'm) {9n\G~^gm), and therefore the equilibrium 
average stress tensor simplifies to 



If the free energy functional F depends only on (f) but not on its derivatives, the stress 
tensor operator ITnmi''^) turns out to be isotropic and is given by Tnmij^) = /n^/m^axs 
(2.17). Then, thanks to the completeness of the basis, the stress tensor can be further 
simplified to (Teq(r)) = kBTS{r). This expression, once properly regularized, gives a 
stress that is independent of system size; that is, the stress is not renormalized by the 
fluctuations in a size-dependent way and therefore no Casimir force can be developed. 
On the contrary, if the stress tensor is not isotropic, as in the case of liquid crystals, 
the result is not trivial and Casimir forces can develop, as shown in [88]. 

All these equations provide expressions for the average fields and fluctuations, ex- 
pressed in terms of the eigenvalues and eigenvectors of the problem, which encode the 
information of the evolution equation together with the boundary conditions. 




(2.36) 
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To summarize, in this Section we have proven that the Casimir force over a body 
is given by 



which is the main result of this Thesis. It shows how to derive the Casimir force from 
the dynamical equations for the field subjected to any kind of noise. It is obtained by 
diagonalizing the evolution operator of the field, and projecting the noise correlation 
and the stress tensor over the set of eigenfunctions. This approach provides the Casimir 
force for both equilibrium and nonequilibrium systems. 

Equation (2.37) shows the well-known nonadditive character of the Casimir force: 
neither the eigenvalues nor the eigenfunctions for different boundary conditions are 
easily related. They cannot be written as a sum of the eigenvalues and eigenfunctions 
of each different problem. This behavior is studied in Chapters 6 and 7 of this Thesis 
for the electromagnetic case. 

The rest of the Chapter deals with applications of Eq. (2.37) to different physical 
systems, in both equilibrium and nonequilibrium situations. 



To show how this formalism works, we calculate the Casimir pressure between two 
plane, infinite plates separated by distance L^. immersed in a medium described by a 
quadratic free energy 



The multiplicative constant /o can be absorbed into 0, and we will eliminate it in 
what follows. The dynamics is described by two transport processes: relaxation and 
diffusion; that is, the Onsager operator is £ = A — DV^, where A and D are the 
transport coefficients (and consequently, positive) associated with the two irreversible 
processes of relaxation and diffusion, respectively. The resulting equation is 




(2.37) 
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'dt 



A0 + DV'0 + e(r,t). 



(2.39) 
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Fluctuation-dissipation is satisfied if the noise is delta correlated in time and the space 
correlation function is 

h{r) = 2kBT{X - DV'^)6{r). (2.40) 

As the energy functional for this system is 0^/2, without terms with spatial derivatives, 
the stress tensor is identical to the local energy functional. Also, the dynamic operator 
is Hermitian, implying that eigenvalues are real and that there is no need to distinguish 
between left and right eigenfunctions. In order to obtain Casimir forces the appropriate 
boundary conditions are of Neumann type, as Dirichlet boundary condition would 
imply trivial vanishing forces. 

We need to solve the eigenfunction problem for the spatial part of the dynamics 
given by Eq. (2.22) in order to calculate the average of the fields that will lead to the 
Casimir pressure over the plates. So, we have to solve the eigenfunction problem given 
by Eq. (2.22) with Ai = X — DV^ obeying Neumann boundary conditions (no-fiux 
boundary conditions), dx(p{0,y,z) = dx4>{Lx,y, z) = 0. 

The normalized eigenfunctions are characterized by three indices Ux, Uy, and n^, 
denoted whole by n, and their form is 

/n(r) = V^e^*^""" if = 

/n(r) = cos (Kx) e'^rn if > 1. (2.41) 

Here, ry = yy + zz and ky = kyj + fc^z. The eigenvalues are 

^i^ = D [kl + kl + kl + kl) , (2.42) 

where k^ = -f-nx, ky = j^Uy, and kz = j^nz, with Ux = 0,1,2, .. . and {ny,nz) E I?. 
The quantity k^^ = ^JD /\ is the characteristic correlation length of the system. 
The average stress is then 

(T..(r)) = ^ E ^Tfl^j ^~"tf fr.{v)f:^{v). (2.43) 

This expression needs to be regularized, otherwise it is divergent. The divergence, 
as explained in [50] [48], is due to the application of the mesoscopic model (2.39) up 
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to very large wavevectors. Conceptually, the stress could be regularized by consid- 
ering generalized hydrodynamic models valid for high wavevectors leading to finite 
stresses, but as the Casimir forces have their origin in the limitation of the fiuctuation 
at small wavevectors, this is not necessary and other procedures are available. There 
are various regularization methods that allow the isolation of the divergent term that 
is independent of the plate separation and therefore cancels out in the computation of 
the Casimir force. The regularization method used in this manuscript is based on the 
Elizalde function detailed in the Appendix 2.8.1. 

To obtain quantitative predictions, we consider specific cases for the noise correla- 
tions. 



2.4.1 Uncorrelated noise in time and space 

We first consider the case of a noise with vanishing correlation time and length, and 
intensity F, i.e., 

(e(r, t)e(r', t')) = T6{r - r')6{t - t'). (2.44) 

This noise correlation, without the — V^5(r) term, automatically puts the system out 
of equilibrium (2.12), as Q 7^ £ + £+. The addition of such a term would have led to 
a stress that was independent of plate separation, not producing a Casimir force [65]. 
Therefore, we consider the effect of the nonequilibrium noise (2.44) on Casimir forces. 
In this case, /i„m(c(/i„) + c(Aim)) = ^^nm and the double sum in (2.43) is reduced. On 
the surface of a plate and applying the limit Ly, L.^ — )■ 00, the stress is given by 

1 



"p noo poo 



dkk ^ , (2.45) 



SirL.D Jo ^^^^ j + A;2 + A;, 



where polar coordinates in the y- and 2;-components have been used. Note that the 
original sum over in (2.43) runs only over N, but the form of the normalizations of 
the eigenfunctions (2.41) allows extension of the sum over Z with a prefactor of 1/2. 
This expression is divergent, so it must be regularized. In order to do so, we use 
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the Chowla-Selberg expression shown in Eq. (2.100). The parameters are s = 1 , 
a = ir/Lx, and = /c^ + Zcq. The first term in the sum of (2.100) equals L^/ ^/k"^ + k"^, 
which combined with the prefactor in Eq. (2.45) yields a term which is independent of 
Lj,., and therefore its contribution to the stress tensor cancels in virtue of Eq. (2.20). It 
must be remarked that the size- independent term is actually divergent if the continuous 
model is assumed to be valid for any wavevector. Then, we are left with the infinite 
sum of modified Bessel functions This sum can be performed analytically, with 

the result 

Fc A = / dk , ; = -■ 2.46 

47rDyo v^PT^e^V^^- - 1 ^o^- 

Let us note that, because the divergence was eliminated, we could have interchanged 
the integral with the summation of the modified Bessel functions to obtain the same 
result. At distances long compared with the correlation length, that is L^, ^ k^^, the 
force decays as 

Fc/A=—^e-'''^'^-. (2.47) 

In the opposite limit, when the plates are at a distance much smaller than the corre- 
lation length, or Lx <^ kQ^ , the force is 

Fc/A = - log (koL,) . (2.48) 

Equation (2.46), as well as Eq. (2.48), shows that the Casimir force diverges if the 
correlation length tends to infinity, i.e., if k^ — )■ 0. This result was obtained in [50] 
using a regularizing kernel technique. The force, as well as short and long distance 
limits are plotted in Fig. 2.2. 

Eq. (2.46) shows that Casimir force diverges if correlation length tends to infinity 
(fco — 0). This is a general property of Casimir forces in Fluctuation-Reaction systems, 
and its origin is the fact that, when ko — )■ 0, does not exist any counter-term which could 
dampens the field fluctuations in (2.39). The total mass, defined as M{t) = J dr(f){r,t) 
perform an unbounded random walk. This problem disappears if we add higher order 
terms to the dynamics, as a 0^ and/or a |V0|^ term in the energy functional, as done 
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0.5 1.0 1.5 2.0 2.5 3.0 

X = koL 

Figure 2.2: Casimir force in units of Fq = in dimensionless units of distance x = k^L 
over a plate immersed in a reaction-diffusion media subject to white noise of intensity F in 
the presence of another plate at a distance L. The exact result is the black curve. The short 
(Eq. (2.69)) and long (Eq. (2.68)) distance limits are the red and blue curve respectively. 

in Ginzburg Landau theory of Critical Phenomena [80]. In this case, it also would be 
needed Renormalization Group techniques [89]. 

2.4.2 Temporally correlated noise 

We next consider the case of a noise that is delta correlated in space but has exponential 
correlation in time 

(r, t) i (r', t')) = T5{v- r') (l + |) e-^\'-'\ (2.49) 

where the factor (1 + |) allows both the white noise limit (a — t- oo) and the quenched 
noise limit (a — > 0) to be taken. Again, the delta correlation in space leads to a term 
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5nm that eliminates one summation in the stress at the plate, which is then given by 

(T..(0)) = -^m^ (2-50) 



If a > 0, we can factorize the quotient as 

" " 1 1 



(2.51) 



with /i„ = kl + ky + kl + /cq as before. We note that this stress is the difference between 
the Casimir stress of two systems with a white temporal noise (2.45) of intensity (1 + 
a)T/a, the first one with = ^ and the second one with kf = ^ + j)- Then, the 
stress on the plate is 

The Casimir force per unit surface on the plate is given just by this expression, because 
the stress on the unbounded side vanishes [as shown by taking the limit — )■ oo in 
Eq. (2.52)]. Finally, we can reobtain the white noise limit if a — )■ oo. 

The case a — )■ corresponds to the quenched limit, where static sources of noise 
are randomly distributed in space. The average normal stress at the wall is 

(T..(0)) = ^E1- (2-53) 
Taking the limit L^, — )■ oo and using polar coordinates. 

Although this expression is finite and does not require a regularization procedure, the 
size-independent contribution can be eliminated using the same regularization proce- 
dure as before, using Eq. (2.100) with s = 2. The result is 

Fr/A - r^^^fc e^-"-(2c.L. + l)-l 

where co = \Jk'^ + k^. After carrying out the integral, we obtain 
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We remark that this system is not dynamically fluctuating, because the noise is quenched 
and the transients have been eliminated. Nevertheless, it creates a Casimir force whose 
origin is the same as previously considered in the sense that the presence of the sec- 
ond plate limits the spectrum of possible fluctuations, and therefore the renormalized 
stresses on the two sides of the plate are different. 

2.4.3 Maximally spatially correlated noise 

As a flnal case, we consider the situation in which the medium is perturbed externally 
by a spatially homogeneous noise, with vanishing correlation time. This could be the 
case when a rapidly fluctuating external fleld is applied to the medium. 

{ar,t)^{v',t'))=T5it-t'). (2.57) 

Applying the same computation procedure as in the other cases, the average local stress 
on each side of the plates is 

r 1 r 

(T^^(O)) = = — , (2.58) 

\ xxK ji 2D2kl 4A' ^ ' 

which is independent of the plate separation. Therefore, the Casimir force vanishes in 
this case. 



2.5 Liquid crystals 

The existence of Casimir forces in liquid crystals has been known for some time now [88] . 
In this Section we apply the formalism presented in Sect. 2.1 to a nematic crystal, 
obtaining the known Casimir force for an equilibrium situation, and expressions for 
the force for some nonequilibrium conditions. The free energy functional of a nematic 
liquid crystal [90] can be written in terms of a planar fleld (f) as 



dr 



y0^ + f(V0)^ , (2.59) 



where we have assumed that the director vector is written in terms of the fleld (f) as 
h = (sin0, 0, cos 0), together with the one-constant approximation (proportional to 
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^2). The first term in Eq. (2.59) comes from a magnetic field directed along the 2;-axis, 
whose intensity is absorbed into ni. 

The simplest dynamical model is obtained with a single relaxational transport co- 
efficient, with the Onsager operator £ = A, leading to 

^ = -A/€i0 + A/caVV + (2.60) 
at 

which is identical in form to (2.39), but with three main differences: the form of the 
ffuctuation-dissipation relation to be in equilibrium, the stress tensor, and the possible 
boundary conditions that produce Casimir forces. Fluctuation-dissipation is realized, 
according to Sect. 2.1, if the noise satisfies 

(e(r, t)e(r', t')) = 2kBT\5{v - v')5{t - t'), (2.61) 

i.e., is purely nonconservative. 

Due to the presence of the gradient terms in the free energy functional (2.59), the 
stress tensor is not isotropic, and therefore even in equilibrium Casimir forces can 
appear. Using Eq. (2.17) the xx component of the stress tensor is 

^- = ^^ +yUJ +yUJ ~^\d-x) ■ ^2.62) 

It is then possible to develop Casimir forces by imposing either Dirichlet or Neumann 
boundary conditions. Dirichlet boundary conditions are equivalent to the strong an- 
choring conditions, i.e., = over the surfaces, and will be the case studied here. 
Casimir forces with Neumann boundary conditions can be easily extracted from the 
Dirichlet ones. 

In this case of Dirichlet boundary conditions, the eigenfunctions of the operator 
M. = \ [ki — K2V^] are given by /n(r) = \J~^ sin {kxx) e*''ii'''ii with eigenvalues 

^ln = \'^2{kl + kl + kl + kl), (2.63) 

where = -f-rix^ ky = j^riy, k^ = j^n^, and ko = 1/^, with indices G N and 
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Because of the boundary conditions, the xx component of the stress tensor at the 
plates is simply given by 

T.M^-'fi^)'. (2.4) 

As in the case of the reaction-diffusion system, we will consider different types of 
noise correlations that, as will be shown below, produce Casimir forces of different 
character. 



2.5.1 Uncorrelated noise in time and space 

We consider an uncorrelated noise as described by (2.44). This case can be considered 
as in equilibrium with a temperature given by F = 2A;bTA. Again the double sum in 
Eq. (2.37) can be reduced, and the stress tensor on the surface of a plate is given by 

Applying the limit Ly, L~, — > oo, we obtain 

2 



— r r°° r°° \ i — 
= Tfi^^ / / E , • (2-66) 

ibVr L^A ^^^^ ^ j ^ ^2 ^ ^2 + ^2 

Using polar coordinates and regularizing the resulting expression using Eqs. (2.103) 
and (2.100) with s = 1, 

12 

,2^fc2+fc2L, _ ^ 

^ 16^ [Li3(e-2'=»^^) + 2A;oL,.Li2(e-2'=«^-) + 2klLlU,{e-^''^'^-)\ (2.67) 

where Lis(z) = ^ is the polylogarithm function. This force, as well as the short 

and long distance limits are plotted in Fig 2.3. At distances long compared with the 
correlation length, that is ^ k^^, the force decays as 

FcM = -^e-''=«^^ (2.68) 
8nXL^ 

In the opposite limit, when the plates are at a distance much smaller than the corre- 
lation length, or <^ /cq \ the force is 



47rA Jo 2y/k^+klL^ _ 



2.5 Liquid crystals 



35 




0.5 1.0 1.5 2.0 2.5 3.0 

X = koL 



Figure 2.3: Casimir force in units of Fq = in dimensionless units of distance x = k^L 
over a plate immersed in a liquid crystal subject to white noise of intensity V in the presence 
of another plate at a distance L. The exact result is the black curve. The short (Eq. (2.69)) 
and long (Eq. (2.68)) distance limits are the red and blue curve respectively. 



This result has already been obtained in the context of liquid crystals in [46] if we 
use the fluctuation-dissipation theorem. It is also the high-temperature limit of the 
electromagnetic Casimir force between two perfect metal plates [14]. 



2.5.2 Temporally correlated noise 

We consider the temporally correlated noise described in Eq. (2.49). By using the 
eigenfunctions of the Dirichlet problem, the stress tensor over a plate takes the value 



V ^ /i„ + a /i„ 
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For any a 7^ 0, the same factorization method as used in Sect. 2.4.2 can be used, leading 
to a Casimir force per unit surface equal to 



fi + -) r 

^' SnaXLl 



Li^^g-2fcoL.) ^ 2A;oL^.Li2(e-2^o^-) + 2^0^1^1(6-^^°^^) 
.Li3(e-2feiL.) - 2A;iL,Li2(e-2'=i^-) - 2kfLlLh{e-^''^^-) 

(2.71) 

where k?: = — and A;? = — + In the limit of infinite correlation length we have 
/co — and ki — )■ ^j^, from which we obtain 

Pc/A = ^l^!}^ [C(3) - Li3(e-2'=^^^) - 2fciL.Li2(e-2^^^^) - 2A;2L2Lii(e-2'=^^^)] . 

(2.72) 

The presented result should be compared with [49], where the same system was studied, 
but a different answer was given [91]. At long distances it decays as in the case of white 
noise (2.69) with a prefactor (1 + |)/a. 

For a — )■ 0, we obtain the quenched limit of the stress tensor at the plates 

/ \ 2 

-r r ^ [t 



Txx = , ^ / dkk V 2.73 

This expression is regularized using Eqs. (2.103) and (2.100) with s = 2, resulting in 

FcA = ——— dk-^ ^ 2.74 

where uj = ^k"^ + k^. This integral can be carried out to obtain the Casimir force as 

^^/^=47rA2^«:2L.e2^o°fli- (2.75) 
In the limit of infinite correlation length we have 

and in the limit of small correlation length {k^L^ ^ 1) we obtain 

Vk 

Fc/A = -^l^e-'^o^^ (2.77) 
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2.5.3 Maximally spatially correlated noise 



As a final case we consider a noise that is rapidly fluctuating in time but that is 
homogeneous in space, described by the correlation (2.57). In this case, hnm is not 
diagonal but is given by 

Jl_(_l)n.][l_(_l)m.] 



The stress on the plates is then given by 



(T..(0)) 



2r g [1 -(-!)"-][!- (-1)"-] 



Avr 



nx,mx=l 



n: 



+ m2 + 2(^)^ 



(2.78) 



(2.79) 



As in the double summation above only odd values of n and m are summed, it can be 
expressed in terms of the Elizalde zeta function over odd numbers, defined as 

1 



{a^{2n + 1)2 + /32(2m + 1)2 + uj^) 



(2.80) 



which can be written in terms of four Elizalde zeta functions. Using the asymptotic 
expansion of the Elizalde zeta functions given in Eq. (2.93) with p = 2 and s = 1, 
the Casimir force can be expressed as an infinite sum of Bessel functions Kq{x) with 
different values of x. The divergent terms, given by the first term in Eq. (2.93), are 
independent of L^,., so the final expression is finite and given by 

Ko {2V2koL,yM + 



(n,m)^(0,0) 



Ko (2V2koL, 



+ 



[Ko (2V2koL,xR 



4 "T 4 



(2.81) 



Two limiting cases can be considered to clarify this result. First, in the limit of long 
distances koL^ ^ 1, the Casimir force is given by 



Fc/A 



-y/2koL 



2\y koL 

whereas in the opposite limit of long correlation length koL ^ 1, the result is 



Fc/A = -ajlog{koL) 



(2.82) 



(2.83) 
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with a = ^{7+ 12y/5 - 20V2) ^ 0.346784, obtained ClS Sb Taylor expansion of the 
associated two dimensional Abel-Plana formula of Eq. (2.81). In the limit of infi- 
nite correlation length this result diverges. Eq (2.81), and the limits of long a short 
correlation length are represented in Fig. 2.4. 




0.1 1 

X = koL 



Figure 2.4: Casimir force in a liquid crystal subject to a spatially correlated noise in 
dimensionless units of distance x = k^L over a plate immersed in a liquid crystal subject to 
a spatially correlated noise of intensity F in the presence of another plate at a distance L. 
The exact result is the black curve (Eq. (2.81)). The short (Eq. (2.83)) and long (Eq. (2.82)) 
distance limits are the red and blue curve respectively. 



2.6 Two non hermitian fields system 

In the two systems we have considered so far (reaction-diffusion 2.4 and liquid crystals 
2.5), the dynamics is described by a Hermitian operator and therefore the potential 
of the method described herein is not fully evident. In this Section, we build a more 
complex system, described by a model with two fields (which could be temperature 
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and concentration, for example) coupled in a non-symmetric way. For simplicity and 
to be concrete we will consider that the fields ipi and 1^2 are scalar, subject to Neumann 
boundary conditions, and with eigenfunctions described by the Fourier modes (2.41). 
In Fourier space, the dynamic equation is 



dt 



,^'2/ \/3k 7k/ 1^2/ \C: 



2 , 



The noises are assumed to be white with different correlation intensities (allowing one 
of them to be set equal to zero later) and no cross-correlation 

(^i(r,t)ei(r',t')) = ri5(r-r')5(t-0, 
(6(r,t)6(r',t')) = r25(r-r')5(t-t'), 

(^i(r,t)e2(r',O) = 0. (2.85) 

Finally, the stress tensor is assumed to be isotropic, depending only on the fields as 

= ftiV^i + «2^2- (2.86) 

As the dynamic matrix is non-Hermitian, the left and right eigenmodes are different, 
being given by 

/i,k(r) = a/I cos (A;.x) 6^*^11-11 ( 



,"k-7k/ 



/2,k(r) = A/^cos(A;,x)e^''i|-^ii '° 



^i,k(r) = A/:^cos(/c^x)e'''i|-'"' 

V \ 



^2,k(r) = -\/-cos(fc,x)e^''ii-^ii | "J^"^'' | (2.87) 

with eigenvalues 

/il,k = ftk, /i2,k = 7k- (2.^ 
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Using these eigenmodes, the different elements needed hnm (2.34) to compute the 
Casimir pressure are 



h 



5] 



kq 



5i 



kq 



/3k 



^^2 



Q!k-7k 
2 



K2 



ak-7k 



(ak-7k)^ 

/3k 
"k-7k 

^2 



(2.89) 



After simple algebra, the stress tensor on the plates is obtained as 

2 



(T..(0)) 



V 



E 



ri(Ki+K2/3^/(ak-7k)^ 2riK2,32/(ak-7k)2 



2au 



"k+7k 



(ri/3g/(«k-7k)^+r2)>^2 

27k 



(2.90) 



which for specific models (that is, specific values of a^, /3k, and 7k) could be computed 
and regularized to obtain the Casimir force on the plates. 

To show the kind of results that can be obtained we consider the simple reaction- 
diffusion two-field model Ok = + Dk^, /3k = A12, and 7k = A2 + Dk"^, with noise 
intensities Fi = F and F2 = 0, representing the system 



dt 
dt 



(2.91) 



Caution should be taken to avoid the case Ai = A2, for which the dynamic matrix is 
not diagonalizable and a Jordan block appears. The method developed in this article 
is not directly applicable, but the generalization is simple. Also, neither Ai or A2 can 
vanish, because there is no damping term to make the nonconservative noise vanish 
and the fields would perform an unbounded random walk. Furthermore, we assume 
that the stress only depends on 'ip2i i-e., T^x = ^^V^i- Therefore, any eventual Casimir 
force is produced by the fluctuations of the second field which are produced by the 
coupling with the first field. The stress on the plates is finally 



(T..(0)) 



2V 



J2 [(Ai + De){\2 + De){\i + A2 + 2Dk' 



As Ai 7^ A2, we can apply partial fraction decomposition to obtain 

1 1 2 



(T..(0)) 



TKXI2 



2DV{Xi - A2)2 



+ 



i.2 _|_ t,2 ^2 _|_ 1.2 1.2 _|_ 1.2 
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where kf = Xi/D, = A2/-D, and A;| = (Ai + X2)/2D. Then, we can perform each 
infinite sum as in the case of scalar white noise to obtain the Casimir force as 

(T m- -^"^^^ / (1 - e-^^^^-) (1 - e-^^^^-) 

^'^^^''^^ ATrDL^iX, - X2Y \ (1 - e-2^3L.)2 

It is interesting to note that, if Ai = and/or A2 = 0, this Casimir force diverges. 



2.7 Discussion and conclusions 

In this Chapter, we have developed a formalism to study Casimir forces in classical 
systems out of equilibrium based on the stochastic dynamical equations of the sys- 
tem under study. The equilibrium case is recovered as a particular limit where the 
fluctuation-dissipation theorem is valid. 

In particular, we study the interaction which appears between intrusions in a 
medium subject to any kind of noise. The study is restricted to additive noise and 
non quantum systems. The method only relies on the stochastic evolution equation of 
the fleld in the medium, and information about the interaction between the medium 
and the intrusions, as given by the boundary conditions of the flelds at the surface 
of the bodies. No assumptions are made regarding any characteristic of the noise, 
which could be internal of external, thermal or induced, white or colored, additive or 
multiplicative, and even non-Gaussian. 

This formalism reduces to the classical thermal Casimir effect when the medium 
is subjected to an additive Gaussian white noise with autocorrelation amplitude Q = 
kBT{C+C^) and its dynamics is described by a Hermitian operator, as shown in (2.36). 

Eq. (2.37) is the main result of the Chapter. It can be used to obtain the Casimir 
force (or fluctuation induced forces) for the steady state of any system, in Equilibrium 
or not equilibrium situations, described by a noise term which can be correlated or 
not in space, colored or white in time and for any geometry in principle. Furthermore, 
as shown in Sect. 2.6, it is possible study systems whose dynamics would be given 
by non Hermitian evolution equations, or even by dynamical operators explicitly time 
dependent [64] [56]. In principle, we can study time dependent Casimir forces with 
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this formalism. Eq. (2.37) can be also used as a numerical computation tool of Casimir 
forces in realist situations: Eigenvalues problem can be numerically solved for any ge- 
ometry and, after introduce the results of this problem in Eq. (2.37) and perform the 
numerical sum, we obtain a numerical result for the Casimir force. All this procedure 
must be done carefully in order to avoid spurious numerical divergences. In fact, an 
implementation of the problem can be easily done with a Finite Element Method im- 
plementation software. First of all, we would solve numerically the eigenvalue problem 
of the geometry over study with the boundary condition of each body and, to carry 
out the (divergent) series over eigenvalues, we add a regularization kernel to the series. 
The introduction of the regularizing kernel transform the considered divergent series 
into an asymptotic series. Therefore we have to carry out this sum with the specific 
techniques of the asymptotic series to obtain the correct result of the series [92]. We 
also used the formalism presented here to probe that Casimir forces nulls in equilibrium 
isotropic systems, and it is generally non zero for any other situation. 

We have obtained an exact formula for the Casimir force felt by a body (2.37), 
which shows how to derive the Casimir force for any geometrical configuration and 
noise. Equation (2.37) can be used to obtain non-equilibrium induced self forces over 
asymmetric bodies, as shown in [51]. It can also provide a numerical tool useful to 
evaluate Casimir forces for complicated geometries. 

Along this Chapter, we have used the formula to obtain the force between parallel 
plates in different media (in a reaction-diffusion model 2.4 and in liquid crystals 2.5) 
under the influence of different kinds of Gaussian noises, i.e., white noise to recover 
the thermal case already studied in the literature, and noises with nonzero spatial or 
temporal correlation lengths, where different forces appears. 

Finally we have shown an example of the evaluation of Casimir forces in a system 
with non-Hermitian evolution dynamics, which was an intractable problem until the 
development of the formalism presented herein. 
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2.8.1 Appendix A: Elizalde zeta function 

The computation of the Casimir forces makes use of the asymptotic expansion of the 
Ehzalde zeta function, which is defined as [93] 

ZA^,A,.)=J: ' (2.92) 

new ^ ' 

we are mainly interested in the case of diagonal matrix A. This function admits an 
asymptotic expansion valid for all complex s [94] 

Z^{s,A,uj) = ^—^ujP-^' (2.93) 



P 



I neZP 

where Ki,{z) is the inhomogeneous Bessel function of the second kind or Macdonald 
function. When p = 2, this is called the Chowla-Selberg formula. The demonstration 
of this formula is as follows. Using the known integral 

r dtf-^e~''' = (2.94) 
Jo 

Elizalde zeta function (2.92) can be written as 

Z,{s,a,,uj) = — / dtf-' J2 e-*-'e-*(""-^-), (2.95) 
We also need the Poisson theorem 

^ (5 (x - n) = e^"*'""'. (2.96) 



If we multiply on both sides of the equality by /(x) = e -^ nj^ q^^^^^ perform an 
integration over x G M^, we obtain 



J2 e-H""-^-) = Yl e-^-"-^-"-. (2.97) 
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After the substitution of (2.97) into (2.95), and separating the term m = 0, we obtain 



Zp{s,ai,uj) 



dtt'-'e 



r(s) Jo 

Finally, after carrying out the t-integral with 



(2.98) 



(2.99) 



and simplifying the result, we obtain Eq. (2.93). We are mainly interested in the case 
p = 1 with a > 0, for which 



v^r(s- i) , 

-UJ + 



47r^ 



1 [s]a ^-^ \auj/ 2 V a 

^ ' n=l 



(2.100) 



r(s)a 

When studying Casimir forces between plates with Dirichlet boundary conditions, the 
following series needs to be computed 



Fi(s,a,w) = ^ 



It is straightforward to obtain 



y L 



and therefore 



Fi(s, a, w) = Zi{s — 1, a, w) — co'^Zi{s, a, u) 



(2.101) 



(2.102) 



(2.103) 



2.8.2 Appendix B: Derivation of Eq. (2.33) 

In this Appendix we will obtain an expression for the cross-average of the mode am- 
plitudes (pnit) in the steady state given in Eq. (2.33) and derived from Eq. (2.31) 

hm (0„(t)C(^)) = lim e-^^"+^'-^' [ dr, [ dr2e^-^^+^'-^^ {Un)iU'r2)) • (2.104) 



00 J —00 



The problem is reduced to the evaluation of the cross-average of the noise amplitudes 

^n{t) and ^ra(t) 

(a(ri)C(^2)) = jdr, jdv2gl{vi)g,^{v2) (^(ri, rOf (ra, r^)) . (2.105) 



2.8 Appendix 
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Spatial and temporal noise correlations factorizes. We can use that to solve (2.105) 
obtaining 

(^n(n)C(^2)) = c(ri-r2) j dvi j rfrsfi-* (ri)/i(ri -r2)fi'm(r2) = c(ri -r2)/i„„, (2.106) 

where hnm is defined. Applying this result in (2.104), the cross-average of the mode 
amplitudes in the steady state is 

lim (0„(t)0:;,(t)) = hm f dn f dr2e-^-'^'-^'^-^-^'-^'^c{n-r2)hnm 



'OO •/ — oo 

OO POD 



dZi / dZ2e-^'-''-^^''c{z^ - Z2)Km, (2.107) 

Jo 

where the change of variables zi = t — ti and Z2 = t — T2 has been applied. After 
performing a second change of variables u = zi — Z2 and v = {zi + Z2)/2, and using the 
parity c{z) = c{—z), we obtain 



t— >oo 



00 roo 



lim {(f)r,it)(f)*^{t)) = / du / t;^;e-/^"(''+"/2)-M;n(«-«/2)c(^) 

J~oo J\u\/2 

/oo roo 
■00 >^|m|/2 

"OO g-(/in-/i^)M/2g-(Mn+M,*J|M|/2 
hnm. I duc(u) 



00 f^n + /^m 

/•oo /"OO 

- — / duc{u)e~^"'' + I duc{u)e~^*^'' 

f^n + l^m iJo Jo 

If we define c{p) as the Laplace transformation of c(t), we obtain Eq. (2.33) 

lim (UtWmit)) = hnm'-^^^^^^^^. (2.108) 
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"/ think I can safely say that nobody understands quantum mechanics. 

Richard P. Feynman 



3 

Stochastic quantization and Casimir forces. 



In this Chapter we show how the stochastic quantization method developed by Parisi 
and Wu can be used to obtain electromagnetic (EM) Casimir forces starting from the 
dynamical approach presented in Chapter 2 of this Thesis. Both quantum and thermal 
fluctuations are taken into account by a Langevin equation for the field. The method 
allows the Casimir force to be obtained directly, derived from the stress tensor instead 
of the free energy. It only requires the spectral decomposition of the Laplacian operator 
in the given geometry. As an application, we compute the Casimir force on the plates 
of a finite piston of arbitrary cross section. Fluctuations of the force are also directly 
obtained, and it is shown that, in the piston case, the variance of the force is twice the 
force squared. 

Fluctuation-induced (Casimir) forces [2] [14] are currently attracting renewed atten- 
tion, probably because of the accessibility of small systems at the nano- and micro-scale, 
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as it is in these small systems where Casimir forces are revealed. However, techniques 
to calculate forces for complicated geometries, beyond the usual ones with high de- 
grees of symmetry, are scarce. Recently, a powerful technique to calculate EM Casimir 
forces [21] [22] has been proposed. It is based on a multiscattering technique and has 
been successfully applied to many configurations, such as plates, cyhnders, spheres, 
wedges, etc. [95] 

Casimir forces have their origin in fluctuations of EM fields. With that idea in 
mind, a method that takes as its starting point the Langevin equation describing the 
evolution of systems subjected to fluctuations has recently been proposed [72]. This 
method was applied to thermal fluctuations, and one of its advantage is that it is 
applicable to systems in or out of equilibrium. However, its application to EM Casimir 
forces of quantum origin has not yet been developed. 

The goal of this Chapter is to apply the Langevin equation method to calculate 
forces of quantum origin. The method is based on the stochastic quantization method 
developed by Parisi and Wu [96] [97]. They construct a Langevin equation for a given 
field subjected to thermal fluctuations. However, this Chapter will show how quantum 
fluctuations arise naturally within that method. 

The plan of this Chapter is as follows. We start in Sect. 3.1 by presenting the 
Stochastic quantization procedure or Parisi-Wu formalism, which let us define the 
quantum fluctuations of the EM field. In Sect. 3.2 we dynamic formalism of Casimir 
effect is applied to the EM field subject to quantum-thermal fluctuations. As a result, a 
formula of the Casimir force over a given object is obtained. In Sect. 3.3, the developed 
formalism is applied to the faces of a perfect metal metal piston of arbitrary cross 
section. As a result, the Casimir force over each plate is obtained for all temperatures. 
The short and large distance limits are studied in the zero and high temperature limits. 
In Sect. 3.5, the variance of the Casimir force is defined and obtained for the system 
under study. We finish with some conclusions and three appendix. In Appendix A 
(Sect. 3.8.1), the intermediate calculations to obtain a formula for the Casimir force 
over a face of the arbitrary shape piston are presented. In Appendix B (Sect. 3.8.2), 
the formula of the variance of the Casimir force is obtained for linear surrounding field 
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in equilibrium. Finally, in Appendix C (Sect. 3.6), we present the equivalence between 
three formalism of Casimir effect: (1) the Partition function formalism; (2) the Stress- 
Tensor formalism and; (3) the Dynamical formalism presented in this Thesis. 
The contents of this Chapter is based on the work published in [98]. 



3.1 Parisi-Wu formalism 

The Casimir force is calculated in [72] via the stress tensor, T, which must be averaged 
over the probability distribution of the fields 0. Such a probability distribution p = P[0] 
is given by: 

P = ^e-^M/n^ (3.1) 

where S[(j)] is the action of the scalar or EM field with zero mass. Wick-rotated in the 
time variable [t = ir) [14] [99] [100], i.e., 

Sm = -- dr j dr<p* + V^j 0, (3.2) 

where 13 = l/ksT and Z is the partition function Z = J VcjyDcf)* e~^^^'^\ For the 
bosonic case, the field must obey periodic boundary conditions in time, that is, 
0(r + Ph,r) = 0(r, r). The stress tensor is normally a bilinear form in the field 
T = 'J'[0*, 0, r], the expression that defines the stress tensor operator f. 

Herein we use the formalism of Parisi and Wu [96] to evaluate the average of the 
stress tensor in an alternative way. The idea of Parisi and Wu consists in a formulation 
of quantum mechanics or quantum field theory in terms of a stochastic process. More 
precisely, a Langevin equation for the field is written as an evolution equation in an 
auxiliary time, which we will call pseudo-time, s. In this description, the field depends 
on the new pseudo-time variable: 0(t, r) — )■ 0(r, r;s). The Langevin equation takes 
the form [101] 

,/ \ '^'S'[0] , . 

a,0(r, r; s) = — + t]{t, r; s) 

^ +VA4> + Vir,r;s). (3.3) 



^2 
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The term 7]{t, r; s) is the source of fluctuations, given by a Gaussian white noise satis- 
fying the fluctuation-dissipation relation [81] 



0, 



(3.4) 



2kBT6{T - T')6{r - r')6{s - s'). 



The associated Fokker-Planck equation (also called advanced Kolmogorov Equa- 
tion) in pseudo-time s for the probability distribution p is 

6 



dsp - 



5<p 



5S[<P] ^ 1 5 



p = 0, 



(3.5) 



5(p f3i 

where we have used the properties of the Gaussian noise given in Eq. (3.4). The steady 
state in pseudo-time is reached in the limit s — )■ oo, where dsP = 0, then p obeys in the 
steady state 

p = 0. (3.6) 



6S[(t)] 1 6 
S(f) ^ ~pJ4> 



Eq. (3.7) is a linear variational equation in p, which admits the solution 

1 



(3.7) 



Then the solution of the Langevin equation in the stationary limit s — )■ oo reproduces 
the probability distribution given by Eq. (3.1) [99]. 

Having an expression for the stress tensor and a Langevin equation for the field, we 
can follow the procedure recently developed in [72] to obtain the Casimir force. The 
field (j) is written as (and a similar decomposition for the noise rj with coefficients rjnm)'- 



r; s) = ^ (pnm{s)gm{r)fn{v), (3.8) 

n,m 

where and Qm are the eigenfunctions: 
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VVri(r) = -Xlfn{r) -^g^dmir) = -a;^^„(r), (3.9) 

which are orthogonal under the scalar product in space or time. The above expres- 
sion indicates that /„(r) and encode the spatial configuration of the system, that is, 
the position of the bodies and their boundary conditions. In a similar fashion, Qmij) 
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and contain the (Wick-rotated) time dependence. As we are considering a bosonic 
field that obeys periodic boundary conditions in r, the eigenvalues are Um = 27:171/ /3Hc, 
m G Z, and the normalized eigenfunctions are gm{T) = ex]i{iuJmT) / y/Wi. Then, the 
coefficients (pnm satisfy the differential equation 

= - [A^ + ^l] <PnUs) + VnUs), (3.10) 

which can be integrated to give 

<Pnm{s) = f rfae(""+-^)('^-^)r/„™(a), (3.11) 
J —00 

where the noise coefficients satisfy, from Eq. (3.4) 

{Vnm{s)r]l,^,{s')) = 2kBT6{s - s')Snn'Smm'- (3.12) 

3.2 Casimir forces from the average stress tensor 

The formalism developed in [72] allows the average stress tensor to be calculated by 
substituting expression (3.8) into the stress tensor and taking the average over the 
fluctuations, rji^r, r; s), in the limit s — )■ oo. 



(T(r)) = hm {T[4>{s),4>{s),r]) = \im Yl {^nurnAs)K,,mM nU,^ Tn^^^]. 

s ^oo s ^ oo 

rii,mi 712 ,m2 

(3.13) 

The two functions (pnm are replaced by their values in terms of the fluctuations, given 
by Eq. (3.11), resulting into a double integral over two pseudo-times. Then, the average 
(T(r)) is carried out using Eq. (3.12), that eliminates one integral, with the final result: 



where Tnm(r) = T[/„, f^,r]. 

Furthermore, the sum over the temporal eigenvalues, Um, can be carried out to 
obtain the result 
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which can be related with the quantum fluctuation-dissipation theorem for the EM 
field [102]. Finally, to obtain the Casimir force over a certain body, the stress tensor 
must be integrated over the surface fl that defines the object 



Equation (3.15) is the main result of this Chapter. It gives an expression for the 
quantum Casimir force including the effects of a finite, non-vanishing temperature, in 
terms of the eigenvalues and eigenfunctions of the Laplacian operator. This expression 
has been easily obtained using the stochastic quantization approach to a quantum field, 
together with the formalism of the Langevin equation calculation of fluctuation-induced 
forces [72]. 

Let us note that, as it occurs in most of the calculations of the Casimir Force, the 
expression for (T(r)) in Eq. (3.15) is generally divergent at every point of space r if 
the sum over eigenvalues runs to infinity. However, the expression for the force, that 
is obtained integrating over the surface of the body, Eq. (3.16), is finite in the sense 
of divergent series, as numerical calculations show. It means that the integration over 
the body regularizes the divergences of the averaged stress tensor. If we assume that 
such regularization is carried out mode by mode, we can interchange the integration 
over the surface and the summation over eigenvalues, to obtain. 



which is a finite result. Therefore, the interchange of the integral and summation 
regularizes the Casimir force, avoiding the use of ultraviolet cutoffs. 

As we show herein, this provides a new method to calculate Casimir forces for 
a given geometry by diagonalizing the Laplace operator. So, this approach is suitable 
for numerical calculations of Casimir forces in complicated, realistic geometries. More- 
over, this method provides the force directly, not as a difference of the free energy with 
respect to a reference state, which is somehow difficult to establish. Also, it can be 
used as a starting point for a perturbative theory for, e.g., non-fiat geometries, rough 
surfaces or similar problems. 




(3.16) 




(3.17) 



3.3 Perfect metal pistons of arbitrary cross section 
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The expression for the Casimir force, Eq. (3.15), allows one to evaluate the quantum 
limit, that is, by setting the temperature equal to zero, where the second summand 
inside the brackets in Eq. (3.15) vanishes, i.e., 

lim(T(r)) = -y^. (3.18) 

n 

In the opposite, classical limit, when /i — )■ 0, a Taylor expansion of the square bracket 
in Eq. (3.15) gives 

\\MT{v)) = ksTY,'^- (3.19) 

n " 

This expression for the Casimir force has been used for classical systems [72], such as 
liquid crystals [46] or reaction-diffusion systems [50]. The two limits, quantum (3.18) 
and thermal (3.19), show that the driving force of the fluctuations has different origin. 
In the first case, the presence of the factor h indicates the quantum nature of the 
fluctuations, whereas in the second case, the factor ksT reveals its thermal origin. 



3.3 Perfect metal pistons of arbitrary cross section 

Let us consider a piston of arbitrary cross section, area A and perimeter P made of a 
perfectly conducting metal surface [103]. Two flat conducting plates of the same cross 
section of the piston are placed at a distance L apart along the x direction. The plates 
are perpendicular to the surface of the cylinder. We have represented the geometry 
under study for the particular case of a circular piston in Fig. 3.1 We calculate the 
Casimir force for these plates by applying the formalism developed in this Chapter. 
First, we solve the eigenvalue problem for the EM field and apply Eqs. (3.14) and (3.16) 
to obtain the force. For the EM field the normal component of the stress tensor reads 

T,, = El + Bl-]^¥?-]^B\ (3.20) 

which has to be averaged over the noise and then integrated over the surface of the 
plates, as shown in Eq. (3.16) and Fig. 3.1. These technical calculations are left to 
the Appendix A 3.8.1 of this Chapter. For perfectly conducting plates, the boundary 
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M ► 

L Loo oo 



Figure 3.1: Geometry of the perfect metal piston system, for the particular case of 
circular cross section. The the cavity and the piston are perfect metals. The piston is at a 
finite distance L of one of the faces of the cavity, and at a distance L^o of the other face. 
This distance tends to infinity in the calculations of the Chapter, but studies of finite Loo are 
straightforward. 

conditions are: E x n = and B ■ n = 0, where n is the normal vector at the surface. 
In this geometry, the EM field can be decomposed into transverse electric (TE) and 
transverse magnetic (TM) modes, which are discussed independently [104]. 

In Coulomb gauge (V ■ A = 0), the electric E and magnetic B fields are derived 
from the potential vector as E = —dfA and B = V x A. 

For the TM modes, the magnetic field is transversal to the x direction, and the 
vector potential A for TM modes can be written (in Coulomb gauge) as 

Atm = i-CVlD,d,CdyD,d^,Cd,D)e-"^\ 

Etm = iuj{-CVlD,d^CdyD,d^Cd,D)e-''^\ 

Btm = -w^(0, -Cd,D, CdyD)e-''^\ (3.21) 

Here u"^ = + X"^, where the fields C{x) and D{r±) satisfy 

d^C{x) = —klC{x), (Neumann BC on x = 0, L) 
ViD„(rx) = -A2D„(r^), (Dirichlet BC on S), (3.22) 

where S is the surface of the cylinder and r_L = {y, z). For the TE set, the electric field 
is transversal to x, so the vector potential is 

Ate = ^^^(0, -S d,N, S dyN)e-"^\ 
Ete = uj\0, -Sd,N,SdyN)e-'''\ 

Bte = iuj{SV\N, -d^S dyN, -d^S d,N)e-''^\ (3.23) 



3.3 Perfect metal pistons of arbitrary cross section 



55 



where the functions S{x) and A^(r_L) satisfy Eqs. (3.22) with the opposite boundary 
conditions: Dirichlet for S{x) and Neumann for A^(r_L). However, in this case, we must 
exclude the constant eigenfunction, with eigenvalue equal to zero, as it gives that 
Ate = and then Ete = Bte = 0. 

Substitution of the TE modes into the expression for the stress tensor, and integra- 
tion over one side of the plates gives, after a long but straightforward calculation left 
to the Appendix A 3.8.1, left to an expression of the Casimir force over one side of the 
piston as 



where = {n^Tr/LY, and are the Matsubara frequencies defined after Eq. (3.9). 
This expression is one of the main results of this Chapter, because we use it to derive 
the Casimir force from noise averages instead zeta regularizations or optical paths. For 
the TM modes, one obtains exactly the same expression, but A„ are the eigenvalues 
of the two-dimensional (2D) Laplacian with Neumann boundary conditions. We will 
denote the complete set of eigenvalues of the Laplacian with Neumann (excluding the 
zero eigenvalue) and Dirichlet boundary conditions by the index p. The expression 
above is the equivalent of Eq. (3.14) when the spectrum can be split into a longitudinal 
and transversal part, that is, = + A^. These series are divergent, but the net 
Casimir force, which is the difference between the force on the two sides of the plate, 
is finite. 

The sum over the variable in Eq. (3.24) can be carried out with the help of 
the Chowla-Selberg summation formula [93]. This formula extracts the divergent, L- 
independent part of the summation, which cancels when the integral in Eq. (3.24) is 
performed for both sides of the piston, resulting in 



Here, A = 27rkBT/hc is the inverse thermal wavelength. If we carry out the sum over 
n, the force results on the more compact form 




(3.24) 



oo 




(3.25) 



p mgZ ra=l 




(3.26) 
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In Ref. [105] it was obtained a formula for the free energy of the configuration considered 
here, that, after differentiation with respect to the distance between the plates, leads 
to the force above. We can proceed to evaluate the Casimir force at T = 0, that is, 
the purely EM case without thermal corrections. This case is obtained by noting that, 
when T — i- (equivalent to the limit A — in Eq. (3.25)), the sum over m can be 
replaced by an integral by the use of Abel-Plana formula [106]. Computing the integral 
for finite A and taking the limit A — )■ 0, the result is 

^^ = -^J2K W^L\) + K^{2nL\)] , (3.27) 

where Ka{x) is the modified Bessel function of order a. This expression gives the finite 
or regularized Casimir force between two plates at distance L. The precise geometry 
of the plates enters into the double set of eigenvalues of the Laplacian (with Neumann 
and Dirichlet boundary conditions) Ap. 



3.3.1 Short distance limit 

For short distances, however, the summation above can be calculated without explicitly 
knowing the eigenvalues of the Laplacian. Such eigenvalues must scale with the inverse 
of the typical size of the piston. So, for distances L much smaller than the section of 
the piston, we can replace the sum over the eigenvalues, Ap, by an integral, using the 
asymptotic expression for the density of states of the Laplacian in two dimensions for 
each set of Dirichlet or Neumann BC [103] [107]: 

p{k) = ^ke{k) + T^feik) + xS{k), (3.28) 

where 6{k) is the step function, A is the area of the piston, P its perimeter, rj = +1 
for Neumann BC and rj = —1 for Dirichlet BC and x depends on the curvature of the 
perimeter and on the number and angle of vertices as 

'^ = ^S£-7l+,^5:/-h.M0i. (3-29) 

I 1- -I j 

where is the angle of each piston's vertex and ^(7^) is the curvature of each smooth 
section of the perimeter. The resulting integrals can be performed to obtain the Casimir 
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force contribution of each mode as 
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(3.30) 
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where the contribution of the zero eigenvalue for the Neumann BC problem must be 
explicitly dropped because, as discussed above, it does not contribute to the Casimir 
force. The global Casimir force is the sum of the two polarizations contributions [103] 

^2 



-he 



" rA + ^(2x--l) 



240L4 



24L2 



(3.31) 



which reduces to the well-known result of the EM Casimir force for infinite parallel 
plates [2]. 



3.3.2 Large distance limit 

In the opposite limit, when L is much larger than the typical size of the plate, the 
argument of the Bessel functions in Eq. (3.27) is much larger than one. Because of the 
exponential behavior of the Bessel functions, only the smallest eigenvalue Ai contributes 
to the sum, with the result 

-T=o = -T^^7iAf e-^^^S (3.32) 



2V7iL 

result directly proportional to the degeneration gi of Ai. Here, a counterintuitive 
result is obtained. One would expect that the thin piston would tend to the known 
one-dimensional (ID) Casimir force, but instead an exponential decay of the force is 
found. The known ID result would be obtained if the zero eigenvalue were considered. 
However, this eigenvalue is excluded, as it leads to a vanishing eigenf unction. Therefore, 
the ID Casimir force cannot be obtained as the limit from three dimensions (3D) to 
ID. 

At intermediate distances, that is, L comparable to the size of the plates, one must 
solve the eigenvalue problem. To illustrate the intermediate behavior and the transition 
from Eq. (3.31) to (3.32), we study the case of a circular cylinder of radius R. In this 
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case, the eigenvalues are the zeros of the Bessel function Juij-) and its derivative, 
with degeneration Qn = 2 except in their lower zero (Ad = 2.40482 for the Dirichlet 
problem, and Ajv = for the Neumann case) The dependence with the radius of the 



a function of L/R when the force is multiplied by R^. Figure 3.2 shows the results 
for the Casimir force, Eq. (3.27), when = 1, = 10, = 100, = 1,000 
and A^ = 10,000 eigenvalues (taking into account their degeneration). We have also 
plotted the two limiting results: (i) the 3D Casimir force (3.31), valid for L <^ R 
with an algebraic behavior: FcR^ oc {L/R)~^ (solid line), (ii) the far distance limit, 
given by Eq. (3.32) with the smallest eigenvalue A?(i?) ~ (^^^)^ = 3.38996/i?^ with 
degeneracy gi = 2, that is valid for L ^ R (dotted line). The transition between both 
regimes is observed at L ^ i?. As expected, when few eigenvalues are summed, for 
instance A^ = 10, the resulting force is only valid in the limit of long distances. As the 
number of eigenvalues increase, the numerical result approaches the 3D Casimir force. 
For A^ = 1,000 eigenvalues, we have excellent results for L/R > 0.05. We remark, 
however, that the full curve can be obtained by only considering A^ ^ 100 eigenvalues 
for large distances and matching this numerical result with the asymptotic expression 
(3.31) for short distances, with a crossover distance L ^ 0.3i?. 

3.4 Thermal Casimir force 

In a similar fashion, we can calculate the thermal Casimir force when — )■ 0, or A — oo. 
Then, in Eq. (3.25), only the term with m = is different from zero. The sum over n 
can be done, with the result 



cylinder is Ap(i?) = A^ 



{R = l)/R^, so the Casimir force as expressed in Eq. (3.27) is 




(3.33) 



p 



For short {L<^R) and long (L^i?) distances, Eq. (3.33) reads 




Fh=o = -ksTgiXie 



-2L\i 



(3.34) 
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R 

Figure 3.2: Adimensional Casimir force on the plates of a cylindrical piston of circular 
cross section of radius i? as a function of the distance L between the plates. Different curves 
are obtained by summing eigenvalues of the Laplacian, using Eq. (3.27)), where N is 
indicated in the legend. Dashed line is the short distance limit and dotted line is the far 
distance result Eq. (3.32). For shorter distances we would need even more eigenvalues to 
converge to the analytical known limit given in Eq. (3.31). 

where gi is the degeneration of the smallest non zero eigenvalue Ai. Here the Area- 
correction term for short distances forces should be read carefully, because the sum 
over n and the integration over k do not commute and gives different results. If sum 
over n is done first, we obtain the result shown above, but if the integration over k is 
done first, the term proportional to (2% — 1) is zero. 



3.5 Fluctuation of Casimir forces 

The Casimir force has its origin in fluctuations, so it is a fluctuating quantity itself. 
The formalism developed in this Chapter allows the calculation of fluctuations of the 
force, defined as 

al = I I ([T(ri) ■ dSi][T(r2) • dS^]) - (3.35) 
Jn Jn 
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Calculation of ap requires the correlation of a product of four noises, 1], which, because 
of the Gaussian nature of rj, factorizes into three products of pairs of noises, as given 
by Eq. (3.4). 

For the piston geometry considered in this Chapter, the fluctuations of the force 
can be obtained by evaluating Eq. (3.35), technical details are left to Appendix B 3.8.2 
of this Chapter. In this case, and because of the geometry of the problem, each of the 
summands that appear in the four-point correlations of the noise gives F^. Therefore, 
the fluctuations of the force, for any temperature and cross section, are simply 



Similar fluctuations have been obtained for a purely thermal force [108] and measured 
in [68] by means of numerical simulation in a hydrodynamical system. The fact that the 
fluctuation of the force is as large as the force itself is a signature of fluctuation-induced 
forces. 

3.6 Equivalence between dynamical approach to Ca- 
simir effect, Stress Tensor formalism and Par- 
tition function formalism 

In this Section, we demonstrate the equivalence of the dynamical formalism of Casimir 
effect in the particular case of equilibrium systems with the Stress-Tensor formal- 
ism [109] of Casimir effect and with the partition function formalism, from which the 
multiscattering formalism is derived [21]. 

As seen in Chaper 2 of this Thesis, for equilibrium systems, the fluctuations of the 
surrounding fleld follows the next Langevin PDE given in Eq. (2.11): 



where C is the Onsager operator associated with transport coefficients, Q is the density 
functional of the energy functional of our problem given in Eq. (2.9), it is 



4 = 2F^. 



(3.36) 



dt(j) = -Cg(j) + ^{r,t) 



(3.37) 




(3.38) 
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and ^(r, t) is a zero mean Gaussian equilibrium noise, with autocorrelation (Eqs. (2.2) 
and (2.12)) 

(e(r, t)e(r', t')) = ksT [C + £+) 5{v - v')5{t - t'). (3.39) 

In this model, (f) is constrained because of the boundary conditions of each a intrusion, 
which in the more general linear case are of Robin type: 

(D«0(s) + iV„9,0(s)) |,e5. = 0, (3.40) 

where is the surface of each a intrusion, and and A^^^^ are constants which defines 
the kind of boundary condition of the surface Sa {Da = and Na 7^ for Neumann 
Boundary Conditions, A^^ = and Da 7^ for Dirichlet Boundary Conditions and 
Na 7^ and Da 7^ for general Robin Boundary Conditions). Then we can define the 
spatial eigenproblem 

CGUv) = /i„/„(r), {DaU^) + NadnU^)) |se5. = 0, (3.41) 

where CQ is a hermitian operator because our system is in equilibrium, then the right 
eigenvalue is gnir) = /*(r), and {f*\fm) = /^^ rfr/*(r)/„(r) = 6nm- 

It can be demonstrated that the partition function of our problem is (Eqs. (2.6) 
and (2.14)) 

Z = J V(Pe-^^^^\ (3.42) 

As -F[0] is a quadratic functional over 0, 2^ is a Gaussian functional integral, then the 
functional integral can be carried out to give [22] 

^ (3.43) 



As we have a system in equilibrium at a constant temperature, we can define the 
Helmholtz free energy of the system in terms of the partition function of the system as 

-l5F = \og{Z) = -Uog\g\. (3.44) 

Using this determinant in term of the eigenvalues of the linear operator and the 
properties of the logarithm, the Helmholtz free energy is [110] [111] 

= ^ E (An) , (3.45) 

n=l 
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where the functional determinant is well defined here because the spectral zeta regu- 
larization principle is used here. The quantities A„ are defined in the eigenproblem 



Gfnir) = A„/„(r) 



(D„/„,(s) + iV,9„/„(s))|,es. =0, 



(3.46) 



where we have assumed that operators CQ and Q share eigenfunctions (that is the case 
in the examples considered in this Thesis), and the relationship between A„ and /i„ is 
assumed to be /i„ = LnXn, where £/„ = These assumptions will be understood 

at the end of the Section. 

From Eq. (3.45), we define the Casimir force as 



'c 



2 ^ A 

^ n=l 



(3.47) 



The problem is to obtain an expression for the derivative of an eigenvalue in terms 
of an internal parameter. To do so, we need an integral representation of the eigen- 
value A„, that can be easily obtained from the eigenvalue problem. If we left-multiply 
the eigenvalue problem by f^ij) and carry out an integration over all the space, the 
eigenvalue can be written as 



A. 



dvf:(v)gUv)=2 / dv¥[U 



(3.48) 



Following [46], let us consider the variation of A„ under a displacement 51. If we have 
into account that a displacement 5C. of a surface implies not only a variation of the field 
in the bulk, but also a variation of the bulk itself, 



5\r 



2 / dr 
'n 



+ 2 (b dsh-6i¥. 

'an 



dfn ' dWU 

We integrate by parts the second summand of the volume integral, obtaining 



2 / dr 
'n 



d¥ ^ f d¥ 



6fn + 2 d) dsn 
'an 



-Sfn + Si¥ 



(3.49) 



(3.50) 



The integrand of the volume integral vanishes, because the equilibrium state is a mini- 
mum of the free energy, and the integrand is the Euler-Lagrange equations of the field 
derived from the energy functional of our problem given in Eq. (3.38). In addition to 
that, we can approximate at first order in 6i that 

Sfn{L) = fn{L + 5i) - UL) = V/„(L) ■ 5i+ O (di") , (3.51) 
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then 

d¥ 



5\n = 2 (t dsiifj, ■ 

an 



(3.52) 



We observe that the tensor operator under brackets is the stress tensor as defined in 
Eq. (2.17), then 

^ = 2/ dsh,.r''[U,f:,s]. (3.53) 
d£u Jan 

Using this resuh in Eq. (3.47), we arrive to our final result 



A^ = V^/ dsh^-r''[U,f:,s], (3.54) 



which is the Eq. (2.37) for equilibrium systems. 

Note that, when C ^ 1, the Casimir force is in this case, following Eq. (2.37): 



Atn + fJ'n Jan 



'c 

n=l 



Using the imposed assumption = LnXn, the Casimir force is in this case, following 
Eq. (2.37): 

Fc = -y;^/ dsh,-r^''[f^,f:,s], (3.56) 
n=i ^" 

which is the same force as the one given in in Eq. (3.54), but with the opposite criterium 
of sign because of the criteria of signs given in Eq (2.20). Note that if the operator Q is 
isotropic in space then its eigenvalues cannot have any information about the geometry 
of the system (A„ 7^ f{L)), and therefore the Casimir force vanishes in this special case, 
as explained in Sect. 2.3 of this Thesis. 

If we use the spectral representation of the Green function of the problem 

^(r,r') = V^^^^^, (3.57) 

Xn 

n 

we recover the Stress-Tensor formalism of equilibrium Casimir effect with the Fluctuation- 
Dissipation Theorem already applied [109]. Then, in this Section we have demonstrated 
the equivalence between the Stress-Tensor formalism, the partition function formalism 
(from which the multiscattering formalism is derived) and the dynamical approach to 
equilibrium Casimir effect. 
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3.7 Conclusions 

We have shown in this Chapter that stochastic quantization together with the Langevin 
formahsm provides a new approach to calculate Casimir forces in the quantum electro- 
dynamics (QED) case, including thermal effects. The starting point is the calculation 
of the Casimir force via the stress tensor, which is a function of the fluctuating fields. 
Parisi and Wu derived a Langevin equation to describe the dynamics of such fields, 
which can be integrated to give an expression for the force. The method presented 
herein is quite simple, and avoids some technical complication of other approaches. 
Moreover, it calculates the force directly, instead of the free energy. Another advantage 
is that it provides a numerical method to calculate the Casimir force in complicated ge- 
ometries, such as those of interest for microelectromechanical systems (MEMS) devices. 
The method only requires spectral decomposition of the Laplacian operator in the given 
geometry, and summation of the eigenvalues and the value of the eigenfunction along 
the boundary of the object, as shown in Eqs. (3.15) and (3.16). Quantum (T — )■ 0) and 
classical {h — > 0) limits are recovered by Eqs. (3.18) and (3.19) respectively. 

In addition, the equivalence of the dynamical formalism in equilibrium cases with 
the partition function formalism and with the stress-tensor formalism have been ob- 
tained and, thanks to the use of the dynamical formalism applied here, we have been 
able to define a variance of the Casimir force and to give a general expression of it. 
This formula have let us to obtain the variance of the arbitrary section perfect metal 
piston at any temperature. 

3.8 Appendices 

3.8.1 Appendix A: Derivation of averaged stress tensor over 
the piston surface 

In this appendix we will obtain the integral over one side of the piston of the average 
over the noise of the stress energy tensor given in Eq. (3.24) as a sum over and A„. 
Because of the system geometry, it is evident that y and z components are zero, then 
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the Casimir force (3.16) over one side of the piston is defined as 

Fc,. = ^ (T..) ds = J^ds (eI + Bl- ^E^ - ^B^^ , (3.58) 

where V is the surface of one side of the piston. The total Casimir force is the sum of 
the contributions of both sides of the piston. Then we need the two point correlation 
function of A, which is obtained from Eq. (3.14) as 

1 1 ctj nfrla^^^fr') 



/3 ^ + A2 II 

where aj^n(r) is the i component of the eigenvector Aa„(r) = — A^a„(r). Then, for the 



electric and magnetic fields we have 

Pi::!^?n + K II ^ 



{Em')) ^ 1 ^,,l,/it^;f\ (3.60) 



{Bi{v)B^{v )) = ^ 2^ ^2 ,x2 112 ' (3-61) 

P + II II 



nm 



where e„(r) = — 2a;a„(r) and b„(r) = V x a„(r). These expressions are the Quantum 
Fluctuation-Dissipation theorem for the EM field [102]. 

From (3.22) and its equivalent result for TM modes, it is immediate to obtain that 
the functions C(x) and S[x) of TM and TE modes respectively, are S[x) = sin(A;a;x) 

and C{x) = cos{kxx), with = {rixir/Ly. Then, for Z{x) = S{x) = C(x), we obtain 

L I 

dxZ^{x) = -. (3.62) 
2 

In addition to that, if gn is solution of the eigenproblem V\gn = — A^^ over the 
surface of one side of the piston V, then a product with g*^ and an integration over V 
gives 

Ai,„ = I ^^^9nf , (3.63) 

II 9n II Jv 

where the norm || gn \\ is defined as 

II 9n f = / ds\gn\^. (3.64) 
Jv 

We separate the study in polarizations because the geometry of the system does not 
mix them. 
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TE modes 

The norm of Ate is defined as 



II Ate ir= I dx I cIsIAteI'^. (3.65) 
Jo Jv 

Then, having into account Eq. (3.23), it is easy to obtain 

II Ate ir= (3-66) 

We obtain (-Ete,x)^ ds = because we are working with TE modes, and 

(-Bte.o;)^ ds = jlp (Ete)^ ds = because these integrals are proportional to sm{kxx) 
over the plate. The contribution of the magnetic field is 

ds{B,,) =-^J:^. ||,^||. ■ (3.67) 

r III II II II 

We study this case in detail as an example of the calculations involved in this Appendix. 
The rest of integrals are obtained on a similar way. Using Eq. (3.23), we obtain 



2 

V JV 



(isb„(s)b* (s) = uj I ds 



S^WInY + (d^Sf (dyNf + (d^Sf (d^Nf , (3.68) 



where S{x) = sin(A;^x), then over the plate we have S{x) = and 
dxS{x) = limx^o kx cosik^x) = k^. Then 

/ dsK{s)K{s) = uj^kl [ ds [{dyNf + id.Nf] = u^kl [ ds iV±Nf . (3.69) 
J-p Jv Jv 

Using Eq. (3.63), the integral is carried out to 

/ dsh^[s)K{s) = u'klXl (3.70) 
Jv 

Finally, using Eqs. (3.66) and (3.70) in Eq. (3.67), we obtain the contribution to the 
Casimir force of the TE polarization as 



' nm " 
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TM modes 

The norm of Atm is, using the definition of Eq. (3.65) and Eq. (3.21): 

II Atm r= (3.72) 

Obviously, (-Btm.x)^ = because we are working with TM modes. The rest of 
integrals over one side of the piston are obtained as 



V 



where we have used that C{x) = cos^k^x), then over the plate C(0) = 1 and dxC(0) = 
lim^^Q —kx sin(A;a.x) = 0. Then, the contribution to the Casimir force of the TM 
polarization, having into account all non zero contributions is 

nm '"■ " 

Eqs. (3.71) and (3.76) are used to obtain the Casimir force over a piston shown in 
Eq. (3.25). 

3.8.2 Appendix B: Derivation of the fluctuations of the force. 

In this appendix we are going to obtain a formula for the variance of the Casimir force 
from Eq. (3.35) and derive the result given in Eq. (3.36) for pistons. 

In the model presented here, the dynamics of the fiuctuations of the medium is 
given by the SDE (Eqs. (2.1) and (3.3)) 



2 



where //(r, r, s) is a zero mean Gaussian noise with autocorrelation 

(?7(r, r, s)rj{T\ r', s')) = V5 {r — r') 5 {v — v') 5 {s — s) . (3.78) 
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Therefore, Eq. (3.77) generalizes the SDEs given in Eq. (3.3) when Ai = —A and in 
Eq. (2.1) when we consider the pseudo-time as the real time s = t and we drop the 
dependence on r. 

We have assumed that each object is subject to an stochastic force F, which is a 
stochastic variable defined as 



■ Ujds, (3.79) 
and whose mean is the so defined Casimir force (Eq. (2.19)) 



(F^) = jf(r-') -n^ds. (3.80) 



al,, = / (fT<> ■ Ujds It^ ■ rikdr) - ( It^ ■ n^ds) . (3.81 



The average is performed over the stress tensor because we are assuming that the 
surface of the bodies does not fluctuate. Then the variance of the Casimir force is 
defined as the variance of the stochastic force F as 

2 

7 

In the case of scalar of EM fields, T*-' is bilinear in the field, then the integrand of the 
Casimir force is proportional to the two point correlation function. We also need the 
four points correlation function in order to obtain the variance. In this derivation, we 
will assume that the spatial operator Ai is an hermitian operator, then its eigenvalues 
fin are real and left and right eigenfunctions are equal, then we have the spatial and 
temporal eigenproblems 

1 ^2 

Mfn{r) = /i„/„(r), -^^f^dmir) = -Umgm{r), (3.82) 

which are a generalization of Eq. (3.9) as explained above. This eigenvalue expansion 
let us factorize the general solution of Eq. (3.77) in the same way as in Eq. (3.8) as 

0(r, r; s) = ^ 0„„(s)5(m(r)/„(r), (3.83) 

n,m 

the same factorization must be applied to the noise. We can apply this factorization 
to Eq. (3.77) to obtain 

ds4>nm = + /^n] 0nm + Vnm{s) , (3.84) 
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whose general solution is 

(t^n^is) = f dae^^"-+^-y^-'\n^{a). (3.85) 

J — oo 

The four points correlation function in the steady state is defined as 

(0(ri,ri)0(r2,r2)0(r3,r3)0(r4,r4)) = lim J] 5Z 5Z 5Z [W'^n^^rnAs)) ^ 



s— >oo 

ni,mi n2,ni2 n3,m3 n4,m4 \fc=l 



4 

X 

k=l 



WfnAT^k)gn,ATk), (3.86) 



where the average of the product of four modes over the noise is obtained by the use 
of Eq. (3.85) as 



\fc=l 



lim / ]^0nfc,mfc(s) ) = lim / dai / da2 / da^ / da^x 



X e ('^"^1 ^''"1 ) (""l e ('^^2 +^""2 ) (<^2 -s) g (w^g +^J„3 ) ((T3 - s) g (w^_^ ) ((74 - s) ^ 

X (r?nimi(<7l)^?n2m2(0-2)?7n3m3(0-3)^n4m4(0-4)) • (3.87) 

Eq. (3.87) depends on the average over the noise of the product of four factorization 
terms of the noise itself. The four point correlation of the modes of the noise is obtained 
by the use of the definition given in Eq. (3.85) as 

(^mmi {<yi)rin2m2 ('72)^?n3m3 {<y:i)Vrnnn (^^4)) = 

/f f nhp rhfi phfi phfi 

dr2 / dvs / dr^ / dri / dT2 / dr^ / dr^x 
J J Jo Jo Jo Jo 

X fn, (ri)/n2 (r2)/n3 (r3)/n4 {^4)9^^ {ri)g^^ {T2)gm, {n)gm^ (^4) X 

X (??(ri, Ti, ai)r7(r2, rg, a2)'n{Y^, rg, o-3)?7(r4, r^, a^)) . (3.88) 

Having into account that ,^(r, t) is a Gaussian noise, the four points correlation function 
of the noise is given as 

(r7(l)r7(2)r/(3)r/(4)) = (r/(l)r/(2))(r7(3)77(4)) + (r7(l)r7(3))(r/(2)r/(4)) + (r/(l)r/(4))(r7(2)r/(3)), 

(3.89) 

where we have used the compact notation rj{i) = ri{ri,Ti,ai). Then, we use the auto- 
correlation noise given in Eq. (3.78) to obtain 

(r/(l)r/(2)77(3)77(4)) = r^~6{l-2)~6{3-4) + r^~6{l-3yS{2-4) + r^~Sil-4)S{2-3), (3.90) 
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where we have used the compact notation 6{i — j) = 6{ri — rj)6{Ti — Tj)6{ai — crj). 
With the result given by Eq. (3.90), we can carry out the integration over the space of 
Eq. (3.88) by the use of the spatial Dirac delta to obtain 



iVnimi (O-O^nama {(^2)Vn3m3 {.(y-Mnirm {cT^)) = 

Sm2,m3S{ai - a4:)6{a2 - (Xs)- (3.91) 



With the last result given in Eq. (3.91), we can carry out the integrals over the pseudo- 
time s of Eq. (3.87) obtaining the last intermediate step of the calculation as 



lim ( TT , 



Kk=l 



^ni ,n2 ^mi ,m2 ^n3 ,714 ^m3 ,1x14 



+ r 



[('^mi + /^ni) + (^ma + /^na)] [(^^^3 + /ins) + {^li4 + /^"4)] 
^rii ,713 ^mi ,m3 ^n2 ,ri4 ^m2 ,m4 



[(^mi + /ini) + (W^3 + /ina)] [(^^^ + /^na) + {(^^4 + f^n4)] 
^ni ,n4 ^mi ,1714 ^n2 ,713 ^m2 ,013 



(3.92) 



Then the four points correlation function defined in Eq. (3.86) results in 



(0(ri, ri)0(r2, r2)0(r3, r3)0(r4, r4)) = 

/ni(ri)/ni(r2) /n3(r3)/n3(r4) / N / N / N / N 

^ ~r~2 > V~r~2 > ^9mAn)9mA'r2)9m3{r2j9m3{U) 

nimi nams V '"l ' V '"-3 ' "■■>/ 

I ^ /ni(l'l)/ni(l'3) /n2(r2)/n2(r4) / n / n / n / n 

+ T 2^ 2^ 6,2 -I ^,2 gm,(ri)^„.,(r3)^n.,(r2)(/^,(r4) 



Here, we assume that ti = T2 = T3 = = r, then the dependence on gmkij) is 
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simplified to 

(0(ri, r)0(r2, r)0(r3, r)0(r4, r)) 



EE 



/ni(ri)/ni(r2) /n3(r3)/n3(r4) 



^ £2 ^ ^ /«i(ri)/„i(r3) /n2(r2)/„2(r4) 

mmi n2m2 

We are interested in the special case with ri = r2 = r and rs = r4 = r', then 



/ni(ri)/ni(r4) /n2(r2)/n2(r3) 



(3.94) 



(0(r,r)0(r,r)0(r',r)</.(r',r)) 



4 ™ ni;^^ Ki + /^"i) K2 + 



Finally we obtain from Eqs. (3.95) and (3.81) the variance as 



0", 



(3.95) 



(3.96) 



EE 



V /T|/„(r), /„.(r)| ■ <ir / T|/„.(r'), /„(r')| • 

+ /^no ) Jv Jr' 



n,i^,n2m2 Rl + /^"O R2 + Z^-^) 

where the second term in Eq. (3.81) has been simplified with the first term of Eq. (3.95). 
It is an integral formula of the the variance of the Casimir force over a given body, one 
of the main results of this Thesis. It is possible to carry out the sums over mi and m2 
using = with m G Z and T = 2kBT, obtaining 



2k bT he 



1 + 



Therefore, the variance can be written as 



Or 



he 



1 + 



e 



he 



(3.97) 



1 + 



^^EE 



jnUr), /„(r)| ■ dv j T|/„(r'), /»(r')l • 



X 



X 



(3.98) 



and the low (quantum) and high (classic) temperature limits of this formula are 
straightforward to obtain as 

^ E E inuv). uv)] ■ dv I nuv'), uv')] . dv'. 



lim (Tp 



(3.99) 
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— /r[/„(r), Ur)] ■ dr (f r[/^(r'), /n(r')] ■ dr' . 

nfJ'm Jr Jr' 

(3.100) 

In the special case Ai = —A, we have /i„ = A^. As a particular case, we obtain the 
variance of the Casimir force of a perfect metal piston of arbitrary section. In this 
particular case, /„(r) = /„(r') and they are constant over each side of the piston. In 
addition to that, 'jr[/„(r), /m(i")] is a bilinear form for the EM field, then we can swap 
n and m indices because they are dummy indices to obtain 



he 

^ ^ 2A„ 



e "B 



he 

2 Arr) 



1 + 



X 



X j(r[/„(r),/„(r)] ■ dr£r[/^(r'),/„(r')] ■ dr' , 
This result admits the simplification 



(3.101) 



E 



he 



) fc oT 



jr[fn{v)Jn{v)]-dv\ . (3.102) 



Using the definition of Casimir force given in Eq. (3.80), we obtain that, for any 
temperature and for pistons of arbitrary section, the variance of the Casimir force is 
twice the Casimir force itself 



cr. 



2Fi 



'F. --c,., (3.103) 
which is the result presented in Eq. (3.36), valid for all perfect metal pistons with 
arbitrary section and for any given temperature. The same result for the variance is 
obtained for classical systems by the application of the same procedure as shown above 
to Eq. (3.100) for any given spatial (linear) operator Ai whose eigenvalues are 

As in the case of evaluation of Casimir forces, we do not need any cutoff to obtain 
a finite result for the variance of the Casimir force. 

The interchange of the integral and summation regularizes the variance Casimir 
force on the same fashion as the Casimir force was regularizated too, avoiding also in 
this case the use of ultraviolet cutoffs. Other regularizations, that in some cases may 
lead to non-universal forces or fluctuations are, for instance, the subtraction of the 
vacuum stress tensor [112] or by averaging the stress tensor over a finite area or a finite 
time [113] [114]. 
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Therefore, we have found an universal form for the fluctuations, as opposite to other 
authors [108] [115]. The difference has its origin in that we do not compute the stress 
fluctuations or the fluctuations of the force on each side of a plate. Rather, we first 
compute the total fluctuating force on the body (which is finite) and then we compute 
its variance. 
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'In mathematics you don't understand things. You just get used to them." 

John Von Neumann 



4 

Extension of Langevin formalism to higher 
temporal derivatives and its application to 

Casimir effect. 



In this Chapter we present another appUcation of our model of evaluation of fluctuation 
induced forces to a generalization of Langevin dynamics. 

Generalizations of Stochastical Differential Equations is an active field in last years. 
For example, the Langevin Difference Equation has an evident interest in numerical 
simulations. The replacement of Gaussian noises by Levi noises is a fundamental tool in 
the study of turbulence and generalizations of Fluctuation-Dissipation theorem. Also 
the study of multiplicative noises have interest in noise-induced phase transitions. 

The generalization presented here has a different nature. We will obtain the two 
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point correlation function of a field whose temporal dynamic is described by a general 
differential functional, i.e. by a linear functional on temporal derivatives of an integer 
order greater or equal to one (that is the case considered before). The case of fractional 
derivatives is left to a future work [101]. 

Usual Langevin equations are a fundamental tool in the study of equilibrium and 
non equilibrium Thermodynamics, but there exist cases when Langevin equations can- 
not describe the relaxation dynamics of our systems. The case of the electromagnetic 
field subject to stochastic sources is a simple first example, note that it is a different 
problem from the case considered in Chapter 3, where a first order temporal derivative 
describes the evolution of the field in a virtual time, while in this second order 

temporal derivative (characteristic of the wave equation) describes the evolution of the 
filed, and its sources are descibed by a white noise. 

We will begin with a description of the dynamics in Sect. 4.1, followed by the 
derivation of the Green function of the field in Sect. 4.2. In Sect. 4.3, we will obtain 
the two point correlation function in the steady state and compare it with the Green 
function of the spatial part of the dynamics. In Sect. 4.4, as an example of the generality 
of the proposed model of evaluation of Casimir forces, we will obtain the Casimir force 
between two parallel plates immersed in these fluctuating media. We finish this Chapter 
with a Discussion in Sect. 4.5. 

The contents of this Chapter have not been published anywhere. 

4.1 Langevin equation with higher temporal deriva- 
tives 

Let us consider the generalization of the linear stochastic differential equation (2.1) for 
the field 0(r, t), 

F{dt)<P{v,t) = -A^0(r,t) +e(r,t), (4.1) 

where is a spatial linear operator and F{dt) is an entire function in dt- The term 
C,{r,t) is a Gaussian noise. Here we assume that C,{r,t) is a zero mean Gaussian white 
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noise characterized by the autocorrelation 

(e(r,t)r(r',t'))=r5(r-r')5(t-0- (4-2) 

Generahzation to colored Gaussian noises can be performed following the same reason- 
ing that in Chapter 2 adapted to the formalism presented here. 

When F{x) = x, we recover the usual Langevin dynamics given in Eq. (2.1), and 
the two point correlation function of the field in the steady state is given by 

hm (0(r,t)0*(r',t)) = ^ ^MM^, (4.3) 

t-s>oo Z ^—^ fin 

k 

where /n(r) are the solutions of the spatial left and right eigenproblem 

A^/„(r) = ^„/„(r), (4.4) 
M+gn{r) = fi:gn{r). (4.5) 

4.2 Green function of the Generalized Langevin Equa- 
tion 

As an intermediate step of the derivation of the two point correlation function of the 
field in the steady state, we need the Green function of Eq. (4.22) without noise. 
This Green function obeys 

[F{dt) +M]G (r, r', t, t') = 5 {r - r') 5(t - t'). (4.6) 

After a symmetric Fourier transformation on time t {f is not transformed) and fac- 
torizing the spatial part of the dynamics by the use of the expansion of eigenfunctions 
given by Eq. (4.4), we obtain that the Fourier transform on time t of the expansion 
terms of the Green function obeys 

iujt' 

[F{-tu)+fin]Gnir,r',u,t') = /„(r)^:(r'), (4.7) 

V 27r 

where the expansion of spatial Dirac delta on eigenfunctions has been used. Then the 
complete Green function can be written as the sum of fractions 

G(r r' CO t') - V ^"^^^^-^^^^ (A 8) 
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As a crucial step, we use the Mittag-Leffler theorem [106]. So that these fractions can 
be decomposed as a sum of simple fractions in obtaining 



G(r,r',w,t') = ^= V 



E 



m=l '"■ 



/n(r)^:(r'). (4.9) 



Here: 



• {w^lm^^ is the countable family of solutions of the dispersion relation of the 
problem F{—iu) = fin for each spatial mode, i.e., they are poles of G in the 
variable u. 

• {Rm}^^ are the correspondent residue for each pole cu^ of G. It is important to 
observe that, if F{—iu) has units of then the Green function has units s". 
As a consequence, these residue must have units of s"~^ because of dimensional 
coherence. 

Then we have obtained the Green function in the temporal frequency space. We 
have to invert the Fourier transform in order to obtain the Green function in the correct 
representation in time as 

G(r, r', t, t') = Y,Y. [^:^e---(*-*')^(sign (Me [<]) {t - t'))] /„(r)^:(r'), (4.10) 

neZ m=l 

where the tilde means that modes with real part equal to zero contribute with a 1/2 
weight. This is the main result of the Section. We observe that these Green functions 
do not explode at infinite time because modes are splitted in three families depending 
of the sign of their real part: 

• If Me [cu^] > 0, then this mode is relevant for all t > t'. 

• If Me [cu^] < 0, then this mode is relevant for all t < t' . 

• If Me [wj^] = 0, then the inverse Fourier transform gives an undamped term 

\e{t-t'). 
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When F{x) is an entire function on x with a zero of order > 1, the Mittag-Leffler 
theorem [106] says that the fraction that we obtained in Eq. (4.8) can be expanded to 

e i e v--\ It 



V , .""^'^ , , (4.11) 



n=l 

whose inverse Fourier transform results in 



-iuj + uon 



{t)=e--^'-''^e{Re[iON]){t-t'))J2RN,n \ . (4.12) 

n=l 



In the important particular case of the usual Langevin equation with first order 
temporal derivative, the known Green function can be derived as a special case of 
Eq. (4.10) as 

G(r,r',t,tO = 5^e-"(*-*')/.(r)^:(r0. (4.13) 

4.3 Correlation function in the steady state 

Once the Green function of the generalized Lanvegin equation has been derived, we 
are able to obtain the correlation function of the field in the steady state. This is the 
key quantity to obtain, because the average of any measurable quantity in the steady 
of the field can be written in terms of this correlation function. 

We can write the general solution of the generalized Lanvegin equation (4.22) in 
terms of the Green function of the field (4.10) as 

/oo /» 
dt' / dr'G{r,r',t,t')^{r',t'). (4.14) 

The correlation function in the steady state is defined as 

(0(r)0*(r')) = lim lim (0(r, t)0*(r', t')) • (4-15) 

t'>t 

This average must be obtained from 

r,t)</.*(r',t'))= (4.16) 

p poo p poo \ 

' dr" / dt" Giv, r", t, t")^iv", t") / dv'" / dt"' ^(r'", ^^^^(r', v'", t', t'") ) . 

J —oo J J —oo I 
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By the use of Eq. (4.2), this integral can be reduced to 

(0(r,t)0*(r',t')) = dr" j dt" G{r,r" ,t,t")G*{r' ,r" ,t' ,t"). (4.17) 

Using Eq. (4.10) and carrying out the integrations in time and space (where the or- 
thogonahty of /„(r) has been used), we obtain that the two point correlation function 
is 

M ' M ' 

{<P{r,t)<P*{r',t')) = rY,Y.Y. RXfn{r)f:{r')x (4.18) 

raeZ m=l m'=l 

X j dt' [e-"-(*-*')e-"r'(*-*')^(sign(Me[<])(t-tO)^(sign(Me[<,])(t-t')) 



Finally, by the use of the definition of the correlation function in the steady state given 
in Eq. (4.15), we obtain after a straightforward calculus 

^'^ on on* 



where the tilde means that modes with real part equal to zero contribute with a 1/2 
weight. Modes with Me [wj^] < do not contribute to the final result because we have 
assumed causality when we imposed the condition t' > t. When Im [w^] ^ 0, having 
into account that the sum is symmetrical in {m,m'}, the result is reduced to 



(0(r)0*(r')) = (4.20) 

AI ' Ad ' 

Re [R'^K^,] {ul^' + u!^/) + Im [K^K^,] {uj^" - uj^/') 



neZ uj" ,>0 uj">0 \ rn ■ mi ■ v m m / 

where we have splitted the modes in real and imaginary parts, uf^ = uj"^ + iw^"- 
Equation (4.20) is the main result of this Chapter. When modes with Me [w^] = 
appear, we have to take into account that they contribute with a 1/2 weight, and 
terms with Me [cjj^] = Me [w^/] = do not contribute at all, because limj;_5.o ;j2^y = 0. 

In the important particular case of the usual Langevin equation with first order 
temporal derivative, the correlation function in the steady state can be derived as a 
special case of Eq. (4.19) as 

(0(r)0*(rO) = r 5^ hMl^_ (4.21) 
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Another interesting case is the correlation function of a field subject to the wave equa- 
tion in the presence of white noise, in this case we have 



C2 



[r,t)-A0(r,t)=e(r,t). (4.22) 



In the temporal frequency domain, the inverse of the wave operator admits the factor- 
ization 

-r2 ^2 r 1 11 

(4.23) 



uj"^ — c^k"^ 2c|k„| 



oj + c|k„| oj — c|k.„| 
then = c|k„| are the positive modes of the dispersion relation, and i?^ = ~2lF7 
its correspondent residue, then the correlation function in the steady state is given by 

(0(r)0-(rO) = ^E ^"^|ki^'^^ ' (4-24) 

where /c„ > is imposed because causality. Similar results can be obtained for even 
greater orders of the temporal derivative. It is easy to verify that, if F{x) = (^f~^ and 
we include a mass term m'^(j) the left part of Eq.(4.22), then the correlation function 
in the steady state has the form 

(0(r)0*(r')) = ^ E (4-25) 
C^(P)^ (k^ + m2)'-? 

where C (p) is a non monotonous increasing function in p, too complicated to be shown 
here. For p = {!,■ ■ ■ , 6} it is easy to obtain C{p) = {2, 8, 18, 32, 50, 36}. 



4.4 Casimir force between infinite parallel plates 

As an application of the developed formalism, we are going to calculate the Casimir 
force between infinite parallel plates, located at x = and x = L, immersed in a 
medium 0(r,t), whose fluctuations in the steady state are described by the generalized 
Langevin equation 

(M.^ =K2A</)-Ki0 + e(r,t), (4.26) 

where ^(r,t) is a white Gaussian noise of intensity F and p is a positive integer. We 
assume that our plates vanish the fluctuations of over their surfaces, then is subject 
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to Dirichlet boundary conditions over the surface of each plate. In addition to that, 
we assume that the free energy is the one of a nematic hquid crystals in presence of 
an external magnetic field given in Eq. (2.59). Then the stress tensor and the used 
eigenf unctions are the ones given in Sect. 2.5, but the average of the stress tensor over 
the surface of a plate is given in terms of the two point correlation function in the 
steady state of this stochastical process, given in Eq. (4.25). 

The stress tensor on the surface of a plate is given by (see Eqs. (2.65) and (2.66)) 

^^^^^0)^ = - E (4-27) 

^^(P) „ (^2 + ^2 + ^2 + ^2)2-- 

Applying the limit L^, — )■ oo, we obtain 

2 



— r /"°° I if' I 



Using polar coordinates and regularizing the resulting expression using Eqs. (2.103) 
and (2.100) with s = 2 — ^, we obtain an infinite series on modified Bessel functions, 
which cannot be carried up to a closed analytical result in general. In any case, it is 
possible to obtain a closed expression for the limits of large and small adimensional 
variable x = k^L. At long distances compared with the correlation length, that is 
L ^ k^^, the force decays as 

Fc/A = ^ M . (4.29) 

In the opposite limit, when the plates are at a distance much smaller than the corre- 
lation length, or L <ti kQ^ , the force is 

rsinf^)rf2) cfi + -) 

If p = 1, we recover the results given (2.68) and (2.69) respectively. For the interesting 
case p = 2, we obtain, at long distances compared with the correlation length L ^ k^^, 
the force 
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and in the opposite limit L <^ kg^, the force is 



^■cM=jg^^. (4.32) 



4.5 Discussion 

In this Chapter we have derived the Casimir force between parallel plated immersed 
in a medium whose fluctuations arc described by a Langevin equation with a complex 
temporal evolution differential term. The result obtained here is an example of the 
kind of calculus for Casimir forces that we can perform within the dynamic formalism 
of Casimir forces developed in this Thesis. 

Obviously, the Langevin equation showed here describes out of equilibrium systems, 
because the fluctuations in equilibrium require a first order temporal derivative to 
describe the temporal evolution of such fluctuations. 

We have also found that the two point correlation function in the steady state is 
an essentially different object of the Green function of such field. These two functions 
are related both each other only for Hermitian systems whose temporal evolution is 
described by a first order temporal derivative. 

The shape of the obtained Green functions is interesting. When an eigenvalue 
could produce a divergence (in the infinite future or past), the Green function cancels 
its contribution in the problematic section of time either future or past. Therefore, 
spatial operators A4 with negative eigenvalues can be treated in this formalism. 

It is important to remark that in this Chapter we do not claim that some kind of 
steady state is reached for the systems studied here. We have studied what would be 
the Casimir force between plates if such steady state exists and it is described by a 
generalized Langevin equation in temporal derivatives. 



Extension of Langevin formalism to higher temporal derivatives and 

84 ITS APPLICATION TO CASIMIR EFFECT. 



MULTISCATTERING FORMALISM 

OF Casimir effect 
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"Nature's great book is written in mathematical language." 

Galileo Galilei 



5 

Introduction to the multiscattering 
formalism of Casimir effect 



In this Chapter we present a short introduction to the multiscattering formahsm of 
Casimir effect. This formahsm was first introduced by Bulgac et. aL in [116]. It was 
generahzed to the electromagnetic (EM) Casimir effect by Emig et. al. in [21] and [22]. 
This formalism is general, and let us to obtain the Casimir energy of a system of 
generalized dielectrics of arbitrary shape for any given temperature. In particular, 
Casimir energy between perfect metal parallel plates [2] and the more general Lifshitz 
formula [14] are derived from this formalism. London energy [3], and retarded limit of 
Van der Waals energy given by Casimir and Polder in [1], and the more general case of 
Feinberg and Sucher potential [117] are covered too. The Casimir energy of a system 
inside a container can also be obtained [118]. 
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The contents of this Chapter is not an original work of the author of this Thesis. It 
is a theoretical introduction of the multiscattering formalism of Casimir effect, needed 
to perform the (original) results obtained in following Chapters of this Thesis. This 
introduction is based on [21] [22] [119]. 



5.1 Partition function description of the problem 

The system under study is described as follows. There is a thermal bath at a given 
temperature T (the case T = is included), within which there are stationary 
dielectrics. The EM action is 

1 



A^{x'), (5.1) 



and the boundary conditions of the EM fields are the continuity across the surface of 
the dielectric of the transversal components of E and H fields 

Ei„ X n = Eout X n, Hi„ x n = H^^^ x fi. (5.2) 

Having into account that these boundary conditions are linear in electric and magnetic 
fields, and that these fields have a linear relationship with the potential vector A'^, then 
the partition function of the constrained A'^ field can be written symbolically as 
r ^ 

Z= VA'^{x^)l[6[C^A^ix^)] eM/'^-X(-^)[^9|-A]A-(x-)^ ^^^3^ 

^ a=l 

The functional Dirac delta for each dielectric a can be written, by means a functional 
Fourier transformation, as a functional integral over the induction currents j'j^{x^) 
inside the volume and surface of each dielectric a. The origin of j'^{x^) are the imposed 
boundary conditions over the surface of each dielectric. 

5 [C„A^(a;^)] = /"l?j-(a;^)eU'^-'ir(-'M^{-'). (5.4) 



Then the partition function of our problem is 

N 



Z = J VA^'ix^) f[ J Vj'^{x^)e^^-f'^''''^'^^'-'''^l^^'-^]^'^'-'''^^^^^^ 

(5.5) 



a=l 
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In order to carry out the functional integrals that appear in Quantum Field Theory 
(QFT), we must apply a Wick rotation of them, which also let us impose the tempera- 
ture of the thermal bath explicitly. This rotation also connect the QFT with Statistical 
Physics by means of the Matsubara formalism [100]. Then we apply the Wick rotation 
by a change of variable to an adimensional imaginary time x° = ict. Then the EM 
action (5.1) Wick rotated is transformed to 

J^S = -^l i'^'j ■l-Al [95 + A] = (5,6) 

Note that the same Wick transformation was applied in Chap. 3 of this Thesis, where 
the Parisi-Wu formalism was used to derive the Casimir force in the EM field. Because 
of Matsubara formalism, bosonic fields (as EM or scalar field) are periodic in the 
interval G [0, phc\, then we can Fourier expand over this interval to obtain 

A'^(x°, x) = ^ A^(x)e*^"^°. (5.7) 

Periodicity condition over x° imposes k„ = with n G Z. In addition to that, if 

n 7^ 0, then y4^„(x) = yl^*(x), then we can factorize the functional differential to 

oo 

VA^{x',^) = I)<(x) n^^^^(x)2^Ar(x), (5.8) 

n=l 

while the euclidean action is expanded to the infinity series 

oo / „ N „ 

SE = PhcJ2[ C?xAr(x) [kI + A] A^,„(X) + J2 ^^X^jXn(^)^n(^) + 
n=0 V"^ a=l 

(5.9) 

where c.c. means the complex conjugate of the written expression. Then the partition 
function is factorized in Matsubara frequencies 



oo 



Z = l[2n, (5.10) 



71=0 

where each partition function has the general form 

N 



Z„^o= / I)A^:(x)PAr(x)n [ Pj-^(x)Pj- w 



X 



X g-§[/'^^^>i;.MW(4+A)Aj;(x)+ELi/««xin,;(^M"W+'=-'=-]. (5.11) 
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Here some caution is needed with the zero Matsubara frequency, because in this case 
we have just a real field without complex conjugate contribution. This will result in a 
1/2 weight in its contribution to the Casimir eff'ect. We can carry out these Gaussian 
integrals in Ai^(x.) to obtain 

+ A| 

where ^o(x, y, ^n) is the dyadic Green function, which takes into account the continuity 
condition of the induction currents. It can be written for imaginary frequency as 

C;o(x, x', k) = Go(x, x', - ® V'Golx, x', n), (5.13) 

where Gnfx, x', k) = -j-^-r^^^^ is the Wick rotated Green function for the scalar field. 
Then Eq. (5.11) is reduced to 

N 
a=l ^ 

(5.14) 



5.2 Free energy of the system 

To solve Eq. (5.14), we apply a multipolar expansion of the dyadic Green function in 
a given base with label i as 

Qoi^, x', K„) = Y,Ci [Mr*(x, /€„) ® Mr"(x', Kn) - Nr*(x, k^) ® Nr"(x', , 

(5.15) 

where we are assuming that x > x'. The quantities M£(x, k„) = V x [^^(x, k„)x] 
and N^(x, /€„) = -V x M^(x, the vector multipoles at imaginary frequency, 

which correspond to the two polarizations of the EM field. TM polarization is given 
by M^(x, K„), while TE polarization is given by N£(x, /t„). ^^(x, k„) is the scalar 
multipole in the same coordinate system. We can find this Dyadic Green functions and 
their multipolar expansions for different coordinate systems in [22] [95] [120] [121]. 
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And we apply a change of variable from induction currents to external multipoles 
moments. These external multipoles moments for an object a are defined in terms of 
the induction currents as 

QM,ii^n) = J dXaj a(Xa)^/lf, (Xq,K„), QE/i^n) = J 'j^XqJq,(Xq)N^ (xQ,,/t„). 

(5.16) 

We will also need the internal multipole moments, defined as 

4'M,i{^'n,) = J dXaj a(xQi)^/li (xQ,/t„), = J dXaja{^a)^i {^a^f^nj- 

(5.17) 

Then, following the formalism developed in [21] [22] [122] [119], the action of the problem 
given in Eq. (5.14) as 

'^J2J2^''^ = ^J2Y1 / / [j*a,^MQo{-^a, K„)j^(X/3) + C.C.] , 

a=l 13=1 a=l /3=1 ^ ^ ^ ^ 

(5.18) 

is transformed in terms of external multipole moments as 

^ N N N N 

^ E E ^^^/^ = T ^ ^ q:m-^Q/., (5.19) 

^ a=l 13=1 ^ a=l 13=1 

with M matrix defined as 

The M"'' matrix contains all the information of our system. Diagonal terms contain the 
inverse of scattering T matrix of each object, which take into account the geometry and 
dielectric properties of the considered object. Non diagonal terms contain translation 
U matrices, which take into account the relative positions and orientations between 
objects. If these objects are immersed inside a container, the formalism changes a bit 
[123] [22] [124]. 

After applying the change of variable to external multipole to the partition function, 
Eq. (5.14) is transformed to 

N oo 



1 V 'JKJ p 

P=E,M a=l 1=1 
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This is a Gaussian integral, which can be carried out to give 



For the zero Matsubara frequency, we have 



because A'^^q{x.) is real 



Helmholtz free energy is defined in terms of the partition function as 



oo ^ oo 

Z = \[Zn J-=--^log(Z„), (5.24) 



n=0 n=0 



where we have factorized Z in Matsubara frequencies given in Eq. (5.10). Using 
Eqs. (5.22) and (5.23), we obtain 

-. CO ' 

^T = ^Y. log|M°^M|, (5.25) 

^n=0 

where the prime indicates that the zero Matsubara frequency has a 1/2 weight. Casimir 
energy is the excess of free energy because of the relative distance between objects. As 
a consequence, we regularize Eq. (5.25) eliminating the contribution to the free energy 
which does not depend on the distance. If we subtract to Eq. (5.25) the energy when 
the objects are infinitely far away form each other, the Casimir energy is defined as 



T 



_ oo ' oo ' 

log|M"/^(/s:„)| - -sY. log|Mf (k„)| , (5.26) 



n=0 ' n=0 



with the definition M^(k„) = lim|x^^|-s.oo M°'^(k„). This is a diagonal matrix whose 



because non-diagonal components 



determinant is equal to |M^(k„)| = n^=i 
are zero (lim|x„^|^oo U"^(/tn) = 0). It is easy to see that Helmholtz free energy can be 
written as 

^ oo ' 

•^T = 7T$^log|l-N(/€„)|, (5.27) 



where N"^ ( characteristic matrix of the system that goes to zero when oo. 

The (1 - W^) matrix can be written as 

(1 - N"^) = + (1 - T"U"''. (5.28) 
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By the use of Abel-Plana formula [106], this Casimir energy can be written as 

he ^ ^ , , he f°° dK 1 , / |l-N(m)| \ 
JF^=— / d/tlogl-NK / log ^ . 5.29 

In the zero temperature limit (T — 0), we obtain the celebrated formula for the 
quantum Casimir energy 

he 

J-o = — / rf^log|l-N(«:)|, (5.30) 

-^71" Jo 

while in the classical limit {h 0), equivalent to the high temperature limit (T oo), 
the Casimir energy is given by the zeroth Matsubara frequency contribution 

J-,; = ^log|l-N(0)|. (5.31) 

In the following chapters, we will use the formulas showed above and variations of 
them to study different properties of EM Casimir, as the appareance of negative en- 
tropy because the Casimir effect between matel spheres, defferent geometries as the 
case of non-parallel cylinders or even different limits, as the Pairwise Summation Ap- 
proximation. 
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"Three's a crowd." 

Anonymous 



6 

Three-body Casimir effects and 
non-monotonic forces 



In this Chapter, we study the not pairwise additive behavior of Casimir interactions. 
This non pairwise additive property leads to collective effects that we study for a 
pair of objects near a conducting wall. We employ a scattering approach [119] to 
compute the interaction in terms of fluctuating multipoles. The wall can lead to a 
non-monotonic force between the objects. For two atoms with anisotropic electric and 
magnetic dipole polarizabilities we demonstrate that this non-monotonic effect results 
from a competition between two and three body interactions. By including higher 
order multipoles we obtain the force between two macroscopic metallic spheres for a 
wide range of sphere separations and distances to the wall. 
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Three-body Casimir effects and non-monotonic forces 



A hallmark property of dispersion forces is their non-additivity which clearly dis- 
tinguishes them from electromagnetic forces between charged particles [125]. Work on 
the interactions between multiple objects is limited mostly to atoms or small particles 
which are described well in dipole approximation [126]. This approximation cannot 
be used for macroscopic objects at separations that are comparable to their size since 
higher order multipole fluctuations have to be included [127] [21]. In such situations, 
also other common "additive" methods such as proximity or two-body-interaction ap- 
proximations fail. Three-body effects for macroscopic bodies have been studied in quasi 
two-dimensional (2D) geometries that are composed of parallel perfect metal cylinders 
of quadratic [109] or circular [128] [129] cross section and parallel sidewalls. For this 
setup non-monotonic forces have been found and interpreted as resulting from a com- 
petition between electric and magnetic polarizations which are decoupled for quasi 2D 
geometries of perfect metal structures. In this Chapter we investigate collective 3-body 
effects between compact objects, including anisotropic polarizabilities, and a wall in 
three dimensions using the multiscattering approach [21] [119] [22]. This allows us to 
observe the influence of polarization coupling and anisotropy on non-monotonic effects. 

We consider the retarded Casimir interaction between a pair of atoms with anisotropic 
electric and magnetic polarizabilities near a conducting wall, see Fig. 6.1. We identify a 
competition between 2- and 3-body effects and prove that this leads to a non-monotonic 
dependence of the force between the atoms on the wall separation H for each of the four 
possible polarizations of fluctuations (electric/magnetic and parallel/perpendicular to 
the wall) separately. For isotropic polarizabilities we find that only the force component 
due to electric fluctuations is non-monotonic in H. 

For atoms, magnetic effects are almost always rather small in the retarded limit. 
Contrary to this, for conducting macroscopic objects contributions from electric and 
magnetic multipole fluctuations are comparable. To study the effect of higher-order 
multipoles, we consider also two perfect metal spheres near a wall, see Fig. 6.1. Based 
on consistent analytical results for large separations and numerical computations at 
smaller distances we flnd a non-monotonic dependence of the force between the spheres 
on H. Unlike for atoms, this effect occurs at sufficiently large sphere separations only. 
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Figure 6.1: Geometry of the two-sphere/atom and sidewall system. Shown are also the 
mirror images (grey) and two- and three-body contributions (sohd and dashed curly lines, 
respectively) . 

We follow this plan for the Chapter: In Sect. 6.1 we present the formalism used to 
describe the system of 2 compact objects (asymmetric atoms or perfect metal spheres) 
in presence of a perfect metal plate. It is a variation of the multiscattering for- 
malism [21] [22] presented in [119], where the multiscattering formalism and image 
method are used together. In Sect. 6.2 the appearance of nonmnotonic forces between 
anisotropic atoms because the presence of a perfect metal plate is studied. In Sect. 6.3 
we study the appearance of nonmonotonic forces between perfect metal spheres in the 
presence of a perfect metal plate, analytically and numerically, and how this effect 
is modified with the distance between the spheres and with the plate. Finally, the 
relevance of the effect presented here is discussed in Sect. 6.4. 

The contents of this Chapter is based on the work published in [130]. 

6.1 General approach 

As presented in Ref. [119] the Casimir energy of two bodies in the presence of a perfectly 
conducting sidewall can be obtained using the scattering approach by employing the 
method of images, which introduces fluctuating currents on the mirror bodies. The 
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Casimir energy of the original system is then given by the energy of the original and 
the image objects and it can be expressed as an integral over imaginary wave number, 



which is given by the T-matrix that relates the regular and scattered electromagnetic 
(EM) fields for each body, and by the U-matrices that describe the interaction between 
the multipoles of object a and object /3, U"'^, and between the multipoles of object a 
and the image of object /3, U^'"^. The T-matrices depend only on the properties of 
the individual bodies such as polarizability, size, and shape. The U-matrices depend 
only on the distance vector between the objects and decay exponentially with distance 
and wave number k. The matrix Moo accounts for the subtraction of the object's 
self-energies and hence follows from M by taking the limit of infinite separations, i.e., 
by setting all U-matrices to zero. For a multipole expansion the matrix elements are 
computed in a vector spherical basis for the EM field with partial wave numbers / > 1, 
m = [119]. 

In the following we study the force F = —dS/dL between the two objects at 
separation L and hence eliminate the contributions to the energy that depend only on 
the sidewall separation H, see Fig. 6.1. We expand the determinant of Eq. (6.1) as 

|M| |M|;^^ = |l + TU^| |l + TU^| X 

X |l - (1 + TUO"^T(U^i +U^'2^)(1 + TU^)~^T(U^2 + U^'^2)| . (6.3) 

The first two determinants on Eq. (6.3) yield together twice the interaction energy 
between a single object and the sidewall, since U''^ = U^'^^ = U^'^^ describes the multi- 
pole coupling between one object and its image and hence depends only on H. Hence, 
we consider only the energy Soo that corresponds to the last determinant of Eq. (6.3) 
and provides the potential energy of the two objects in the presence of the sidewall so 
that F = —dSo^/dL. In the absence of the sidewall, H — )■ oo, the matrices U'^'"'' all 




(6.1) 



with the matrix 




(6.2) 
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vanish and £oo_ simplifies to the energy between two spheres [21]. For an interpretation 
in terms of multiple scatterings, it is instructive to use the relation log 1^41 = Trlog(74) 
and to Taylor expand the logarithm and the inverse matrices, 

= -- rd^y-Tr ( Sr=o(-ir(™0"T(U2^ + U^'^^)x Y 
- 27rio ''^P 1^ x^^^o(-l)"'(TU0"'T(Ui2 + U^'i2) j 

The trace acts on an alternating product of T- and U-matrices which describe scattering 
and free propagation of EM fluctuations, respectively. Multiple scatterings between an 
object and its image (TU'^) are followed by a propagation to the other object-image 
pair, either to the object (U^^) or its image (U^'^^), between which again multiple 
scatterings occur before the fluctuations are scattered back to the initial object or its 
image (U^^ or U^'^^) and the process repeats. This expansion is useful for small objects 
or large separations. 



6.2 Two atoms 

First, we consider the case of two identical, ground state atoms near a wall, see Fig. 6.1. 
The separation between the atoms is L and the separation of each of them from the wall 
is H. In dipole approximation, the retarded limit of the interaction is described by the 
static electric (a^, q;||) and magnetic {(3^, f3\i) dipole polarizabilities of the atoms which 
can be different in the directions perpendicular [z) and parallel (||) to the wall. The 
T-matrix of the atoms is diagonal and has finite elements only for the dipole channel 
(partial waves with / = 1), given by 

Tfo = ^«./^^ Tit = ^«||/t' (6.5) 

for electric and 

T^^ = l^^', T,Z=Ip^^^' (6.6) 

for magnetic polarization with m = ±1. For atoms, the polarizability is much smaller 
than L^, and hence it is sufficient to compute the interaction to second order in the 
polarizabilities. This amounts to neglecting all terms other than p = 1 and n = n' = 
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in Eq. (6.4). The resulting energy £^oo is then compared to the well-known Casimir- 
Polder (CP) interaction energy between two atoms (without the wall), 

he 



[33«J +13al -Ua\\P, + 33/3J +13(3^ -14/3||a^] 



(6.7) 



which corresponds to the sequence TU^^TU^^ in Eq. (6.4). The total interaction energy 



IS 



(6. 



he 



87rD'< 



with D = y/L"^ + 4iJ2. The 2-body energy £2,\iD, L) comes from the sequence TU^'^'^^TU^' 
in Eq. (6.4) and hence is the usual CP interaction between one atom and the image of 
the other atom (see Fig. 6.1). The change in the relative orientation of the atoms with 
I = LjD leads to the modified CP potential 

26aJ +20^2 -14£2(4aJ - Qaya^ + ^al)^^3t{a\\ - 
+26/3|f +20/32 -14£2(4/3|| - 9/3|j/3, + 5/32) + 63£4(/3|| - 
-14(a|,/3|,(l-£2) +e'^a\i(3,) - 14(/3||a|, (l-f^) +fP\^a,) 

"(6.9) 

The 3-body energy £^{D,L) corresponds to the matrix products TU^-'^TU^'^^ and 
TU^'^^TU^^ in Eq. (6.4) and hence describes the collective interaction between the two 
atoms and one image atom. It is given by 

(3^6 + i5£5 ^ 28£^ + 20£3 + 6£2 - 5£ - 1) (^aj - /^j) 
- {3e + 15£5 + 2Ae - 10£2 - 5£ - 1) {al - f3f) 

(6.10) 

For isotropic electric polarizable atoms this result agrees with that of Ref. [126]. It 
is instructive to consider the two limits H <^ L and H ^ L. For H <^ L one has 
D ^ L and the 2-body potentials are identical, £2^\{L,L) = £2^\{L). The 3-body 
energy becomes 



12 



£3iD,L) 



Ahe 



1 



£3(1, L) 



he 



47rL 



-[-33aJ + 13a^ + Ua\\l3, + 33/3J - 13/3^^ - 14/3||a^] . (6.11) 



The total energy £"00 is now twice the energy of Eq. (6.7) plus the energy of Eq. (6.11), 
and hence 

he 



lim £oo(L, H) 



SttL^ 



[I3(2a,)2 + 33(2/3,1)2 _ i4(2/3||)(2a,)] . (6.12) 
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Then, £^oo becomes the CP potential of Eq. (6.7) with the replacements — > 2a^, 
a\\ — )■ 0, — )■ 0, /3|| — )■ 2/3||, as can be observed in Fig. 6.3. The 2-body and 3-body 
contributions add constructively or destructively, depending on the relative orientation 
of a dipole and its image which together form a dipole of zero or twice the original 
strength (see Figs. 6.2 and 6.3). For H ^ L the leading correction to the CP potential 
of Eq. (6.8) comes from the 3-body energy which in this limit becomes (up to order 

The signs of the polarizabilities in the leading term ~ can be understood from the 
relative orientation of the dipole of one atom and the image dipole of the other atom, 
see Fig. 6.2. If these two electric dipoles are almost perpendicular to their distance 
vector they contribute attractively to the potential between the two original atoms, 
but if these electric dipoles are almost parallel to their distance vector they yield a 
repulsive contribution. Contrary to the electric behavior, if these two magnetic dipoles 
are almost perpendicular to their distance vector they contribute repulsively to the 
potential between the two original atoms, but if these magnetic dipoles are almost 
parallel to their distance vector they yield a attractive contribution. For isotropic 
polarizabilities, the leading term of Eq. (6.13) vanishes and the electric part ~ of 
the 3-body energy is always repulsive, while the magnetic contribution ~ of the 
3-body energy is always attractive. 

The above results show how the force between the two particles varies with H. If the 
two particles have only either or polarizability, their attractive force is reduced 
when they approach the wall from large H due to the repulsive 3-body interaction, as 
shown in Fig. 6.3. At close proximity to the wall, the fluctuations of the dipole and 
its image add up to yield a force between the particles that is enhanced by a factor of 
4 compared to the force for if — t- oo. Corresponding arguments show that the force 
between particles with either ay or polarizability is enhanced at large H and reduced 
to zero for — 7- 0. This proves that the force between particles which both have either 
of the four polarizabilities is always non-monotonic. The situation can be different if 



Ss{H,L) 



he 
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cancel 
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Figure 6.2: Typical orientations of electric (E) and magnetic (M) dipoles and image 
dipoles for H/L — )• and H/L — )• cxd. 



more than one polarizability is finite, especially for isotropic particles. In the latter 
case all contributions (electric, magnetic, mixed) are enhanced for if — and only the 
electric term is reduced at large H so that only the electric part gives a non-monotonic 
force (Fig. 6.3). In general, the monotonicity property depends on the relative strength 
and anisotropy of the electric and magnetic polarizabilities. 




Figure 6.3: Contributions to the energy between 2 atoms described by Eqs. (6.5) and 
(6.6), splitted in contributions of electric polarizabilities (left figure) and magnetic polariz- 
abilities (right figure). In blue, the perpendicular to the plate (z) contribution, in red the 
parallel to the plate (||) contribution and, in black, the energy for isotropic polarizabilities. 
The limits for H <^ L (dotted curve), H ^ L (dashed curve) and the limit H ^ oo (dot- 
dashed curve) are also plotted. As the limits for H <^ L and H ^ L show, all figures plotted 
here are nonmonotonous with H, except the case of isotropic magnetic polarizabilities (black 
curve of the right plot). 
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6.3 Two macroscopic metallic spheres 



6.3.1 Analytical results 

Secondly, we study two macroscopic perfect metallic spheres of radius R for the same 
geometry as in the case of atoms where the lengths L and H are measured now from 
the centers of the spheres, see Fig. 6.1. The T-matrix is diagonal and the elements 



T, 



M 
Im 



2K, 



(6.14) 



T, 



Im 



(6.15) 



,7r Ii+i/2if^R) + 2KRIl^,^^iKR) 
' ' 2 Ki+i/2{kR) + 2kRK[^^/^{kR) 

are given in terms of the modified Bessel functions Ip, K^. First, we expand the energy 

in powers of R by using Eq. (6.4) which implies that we expand the T- matrices for 

small frequencies but use the exact expressions for the U-matrices. For R <ti L, H and 

arbitrary H/L the result for the force can be written as 

i+2 



he 



IT 



j=6 



(6.16) 



The functions fj can be computed exactly. We have obtained them up to j = 11 and 
the first three are (with s = \/l + Ah"^) 



h{h) 

f7{h) 
fs{h) 



s"^(18 + 312/1^ + 2052/i^ + 6048/i^ + 5719/i^) + 18 - 12/^^ + lOOl/i*" 



160/1 



12 



5-11(6210 + U055Ah^ + 1315364/;,^ + 6500242/i^ 



-11/ 



17830560/1*^ + 25611168/^1° +15000675/;, 



,12^ 



-6210 - 3934/^2 + 7QAh^ - 78/i^ + 71523/i 



12 



(6.17) 



The coefficient fj of R^ vanishes since a multipole of order / contributes to the T-matrix 
at order /2^'+i so that beyond the two-dipole term ~ R^ the next term comes from a 
dipole (/ = 1) and a quadrupole (/ = 2), yielding /§. For H ^ L one has fei^h) = 
-1001/16+3/(4/i6)+C(/i-8), /g(/i) = -71523/160+39/(80/i6)+C(/i-S) so that the wall 
induces weak repulsive corrections. For H L, fei^h) = —791/8 + 6741/i^/8 + 0{h^), 
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/8(/i) = -60939/80 + 582879/iV80 + C»(/i^) so that the force amphtude decreases when 
the spheres are moved a small distance away from the wall. This proves the existence 
of a minimum in the force amplitude as a function of H/R for fixed, sufficiently small 
i?/L, see Fig. 6.4. We note that all fj{h) are finite for — > oo but some diverge for 
/i — )■ 0, e.g., /g ~ /ii ~ /i^^, making them important for small H. 




0.5 1.0 1^ — H 1.5 2.0 05 1.0 1-^ — !l 1.5 20 



Figure 6.4: Functions fei^h) (left figure) and fsih) (right figure) (Eq. (6.17)) as a function 
of h = j^. The black curve is the exact function. The series expansion for H <^ L (blue 
curve), H ^ L (red curve) are also plotted, and the limits /i — t- (dashed blue curve), and 
/i — >• oo (dashed red curve) are also represented to show the nonmonotonicity of f&{h) and 
h{h). 



6.3.2 Numerical results 

To obtain the interaction at smaller separations or larger radius, we have computed 
the energy £oo and force F = —dSoo/dL between the spheres numerically. For the 
energy, we have computed the last determinant of Eq. (6.3) and the integral over k of 
Eq. (6.1) numerically. The force is obtained by polynomial interpolation of the data 
for the energy. The matrices are truncated at a sufficiently large number of partial 
waves (with a maximum truncation order l^ax = 17 for the smallest separation) so 
that the relative accuracy of the values for Soo^ is ~ 10~^. The data for H/R = 1 
are obtained by extrapolation in Imax- The results are shown in Figs. 6.5 and 6.6. In 
order to show the effect of the wall, the figures display the energy and force normalized 
to the results for two spheres without a wall. Fig. 6.5 shows the energy and force 
as a function of the (inverse) separation between the spheres for different fixed wall 
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distances. The proximity of the wall generally increases the interaction energy and 
the force between the two spheres. The effect is more pronounced, the further the 
two spheres are separated. For sufficiently large H/R, the energy and force ratios are 
non-monotonic in L and can be slightly smaller than they would be in the absence of 
the wall. Fig. 6.6 shows the force between the two spheres as a function of the wall 
distance for fixed L. When the spheres approach the wall, the force first decreases 
slightly if R/L < 0.3 and then increases strongly under a further reduction of H. For 
R/L > 0.3 the force increases monotonically as the spheres approach the wall. This 
agrees with the prediction of the large distance expansion. The expansion of Eq. (6.16) 
with j = 10 terms is also shown in Fig. 6.6 for R/L < 0.2. Its validity is limited to 
large L/R and not too small H/R; it fails completely for R/L > 0.2 and hence is not 
shown in this range. 




Figure 6.5: Numerical results for the potential energy (dashed curves) and force (solid 
curves) between two spheres as function of R/L for different sidewall separations H/R. Both 
force and energy are normalized to their values in the absence of the sidewall. Inset: Magni- 
fication of behavior for small R/L. 
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{H -R)/R 



Figure 6.6: Numerical results for the force (dots) between two spheres as function of the 
sidewall separation H/R for different sphere separations R/L. Shown are also the analytical 
results of Eq. (6.16), including terms up to j = 10 for R/L < 0.2 (solid curves). Inset: 
Magnification of the non-monotonicity. 



6.4 Discussion 



In this Chapter we have demonstrated that nonmonotonicities of Casimir force appears 
not only between cyUnders [109] [128] [129], but also between compact objects. In par- 
ticular, nonmonotonicities of Casimir force appears between two atoms and between 
two perfect metal spheres, in presence of a perfect metal plate. For perfect metal 
spheres, there exists a critical distance from which this nonmonotinicity of the force 
disappears, although the variation (enhancement in this case) of the force between the 
spheres continues because the perfect metal plate contribution is not pairwise additive. 
As shown by Eq. (6.3), Casimir energy between three objects is not pairwise additive, 
but the effect of the inclusion of a third object in the system is not a simple enhance- 
ment of the force between objects (although in general this seems to be the case). We 
have shown here that, depending of the relative distance between the three objects, 
the effect of the non-pairwise behavior can be either an increase or a decrease of the 
force between objects, when it is compared with the isolated case. 

The results presented here for isotropic atoms are potentially relevant to the inter- 
action between trapped Bose- Einstein condensates and a surface [131] at close surface 
separations. The results for macroscopic spheres could be important for the design of 
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nano-mechanical devices where small components operate in close vicinity to metal- 
lic boundaries. Generally, the wall-induced enhancement of the interaction can make 
the experimental observation of Casimir forces between small particles more feasible. 
The dependence on the anisotropy of polarizabilities applies not only to atoms but to 
general polarizable objects and suggests interesting effects for objects of non-spherical 
shape, e.g. spheroids. 
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'Science is built upon approaches tliat gradually approach the truth. 

Isaac Asimov - Nightfall 



7 

Pairwise Summation Approximation of 

Casimir energy 



In this Chapter, we obtain the Pairwise Summation Approximation (PSA) of the 
Casimir energy from first principles in the soft dielectric and soft diamagnetic limit, 
this analysis let us find that the PSA is an asymptotic approximation of the Casimir 
energy valid for large distances between the objects. We also obtain the PSA for the 
electromagnetic (EM) coupling part of the Casimir energy, so we are able to complete 
the PSA limit for the first time for the complete electromagnetic field. 

109 
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7.1 Introduction 

Since 1948, when Casimir introduced the energy that got his name [2], calculation 
formulas have been looked for. Many analytical and numerical methods have been 
proposed, such as the zeta function technique, the heat kernel method, semiclassical 
methods or Green function (local) methods just to mention a few of them [14]. However, 
exact results have been obtained only for some simple geometries. 

The asymptotic analysis of Casimir energies has also a long history. In fact, Casimir 
himself and Polder in year 1948 gave the first asymptotic formula for the Casimir en- 
ergy between two electrically neutral bodies in terms of their electric induced dipoles 
[1]. Some time later, a generalization of that formula, known as the Pairwise Summa- 
tion Approximation (PSA) was derived [15] for electric media. The main assumption 
is a linear superposition of the Casimir-Polder interactions between the induced po- 
larizabilities of each element of volume body. Then, the PSA energy is expressed as 
an integral over the two object's volumes and it is proportional to the objects polariz- 
abilities. The formula has been recently reobtained by R. Golestanian in [26] and by 
K. A. Milton et. al. in [132] in the soft dielectric limit. 

Besides, a new asymptotic method for calculating Casimir energies in term of the 
induced multipoles of the interacting bodies has been proposed in [21] and [116]. This 
formula provides a procedure for the calculation of the Casimir energy between N 
arbitrary shaped compact objects [130]. 

In this Chapter, the asymptotic calculation of the electromagnetic Casimir energy 
and a presentation of a systematic asymptotic expansion procedure for the integral 
dipole formula on higher orders is presented. We will reobtain the classical results of 
Casimir and Polder [1], and of Feinberg-Sucher in [117] in the rarefied or soft dielectric 
and soft diamagnetic limit. 

The result we obtain here is a generalization of Milton et. al. one [132] but assuming 
that the bodies are soft diamagnetic as well as soft dielectric. For this purpose, we 
will use the multi-scattering expansion of the Casimir energy formalism given in [21]. 
To our knowledge, it is the first time this formalism is used to derive the PSA. We 
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obtain that the PSA is the first order of a perturbation expansion in the difference 
between the electric and magnetic permeabihty constants of the objects respect the 
electric and magnetic permeability constants of the medium were the objects are placed. 
The interest of this result is that we obtain this pairwise summation formula for the 
complete electromagnetic Casimir energy which is derived as an asymptotic limit of 
an exact and free of divergences formula. That means that now we can establish the 
range of validity for that approximation. In fact, we will justify why this formula is 
valid in the far distance objects limit. This derivation also gives us a perturbative 
procedure for corrections of this approximation and posterior expansions to more than 
two objects or finite temperature cases. 

We will follow this plan for the Chapter: Using the soft dielectric and soft diamag- 
netic approximation, in Sect. 7.2 we will obtain the PSA of Casimir energy for the 
zero temperature case starting from the exact Casimir energy formula given in [21]. 
We will also study the far distance limit to reobtain the asymptotic Casimir energies 
given in [1] and [117]. In Sect. 7.3 we will obtain PSA Casimir energy formulas for 
any temperature, studying the high and low temperature limits. In Sect. 7.4 we will 
study the PSA Casimir energy for the three bodies system and for the general N bodies 
system. We will obtain a kind of superposition principle of Casimir energy in the PSA 
approximation. Finally, in Sect. 7.5 we will study the first perturbation energy term 
to the PSA approximation and we will discuss about the nature of the PSA limit of 
the Casimir energy. 

The contents of this Chapter is based on the work published in [28], [133] and [134]. 

7.2 Diluted limit at zero temperature 

Our objective is the calculation of the complete electromagnetic Casimir energy be- 
tween two bodies in the soft dielectric limit. For this purpose, we will use the Casimir 
energy formula between two compact bodies at temperature T = 0, given in [21] as 




(7.1) 
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Where k is the frequency, 1 is the identity matrix, N is the matrix N = T1U12T2U21. 
Here, is the T scattering matrix of the a - body under the electromagnetic field, and 
Va/B is the propagation matrix of the electromagnetic field from a - body to /3 - body. 
We will use the position representation instead the multipole representation used in 
[21], so we can identify Vap = Goa/3, where is the free dyadic Green function. 

Taking into account that log 1^4 1 = Trlog(A) and that log(l — x) = —J2'^=i 
transform Eq. (7.1) into 

Equation (7.2) is an asymptotic expansion of Eq. (7.1) in the Pmax < 00 case, which 
means that our calculus will be valid in the large bodies distance limit. The T operator 
is related with the potential V (we will see what V is later) by the Lippmann - Schwinger 
equation, that can be written as [135] 

T = (1 - VGo)-^ V. (7.3) 

Applying a Born expansion to Eq. (7.3), we can obtain an approximation of the T 
operator in the soft dielectric and soft diamagnetic limit as 

00 

T = J2iyGo)"V ~ V. (7.4) 

n=0 

This approximation is more valid for weaker V, so here is where the soft dielectric and 
soft diamagnetic limit is applied. 

Now we will study the lowest expansion order. In the lowest expansion order {p = 1 
and T = V^) we get the asymptotic approximation of the Casimir energy between two 
bodies as 

he 

E = -i^j^ rfA;Tr(ViGoi2Wo2i). (7.5) 

We can separate the magnetic and electric part of each body potential (which is a 
diagonal operator in positions space, because it is local) and each field contribution in 
the free dyadic Green function. 



7.2 Diluted limit at zero temperature 



113 



To define the potential Vi, we use the generalized constitutive relations of materials 
with magnetoelectric coupling 



D = eE + aH, 
B = /3E + /xH. 



(7.6) 



In general, aij = Pij, because a linear magnetoelectric coupling is allowed only in 
materials odd under time reversal or odd under inversion {aij — ?■ —aij). The magneto- 
electric polarizability aij is zero otherwise. The potential of each body can be defined 
as the difference of energy of the EM field because the existence of this body. Having 
into account that the energy of the EM is defined as 



E 



rfx'^ (E ■ D + H • B) 



(7.7) 



where D = cqE and H = yUoB in the vacuum, the EM energy in presence of N 
generalized dielectrics is 



1 r ^ 

S=- dx^ (eoE^ + ^xoH^) + J2 ^Si 



{71 



We use the generalized constitutive relations given in Eq. (7.6) to obtain the excess of 
energy because the existence of each dielectric as 



AS,; 



dx^ (E, H) 



(e* - eo) 


ai 


A 


{l^i - A^o) 




(7.9) 



then the potential is defined as 



Vi 









ai \ 


V A 


f^i J 



Xi [r) , 



(7.10) 



where q = ej — eo, /ij = /ij — yUo in the vacuum and Xi (r) is the characteristic function 
of the i-body volume (1 inside the body and in the rest of the space). The definition 
presented here let us define the potential Vi when the medium is not the vacuum 
and when the objects are not isotropic (think for example in uniaxial and biaxial 
crystals [136]). Note that we represent the electromagnetic properties of the objects 
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by their potential energy instead the complement representation by their boundary 
conditions [18]. It will let us to obtain the PSA as an integral over the volume of the 
bodies as required by the PSA. The contributions of the free dyadic Green function 
are obtained in the Appendix in Eq. (7.113) - (7.116). 

As seen in Eq. (7.113) - (7.116) of the Appendix, operators are defined over three 
different linear spaces: 1) An EH - space, whose components are the electric and the 
magnetic field; 2) over the space coordinates, because we are working with a vector 
and a pseudovector; and 3) over positions. We must solve the trace over these three 
spaces: EH - space, vector coordinate space and position space. First we solve the 
trace in the EH - space for the important case of non-magnetoelectric coupled objects 
and we obtain: 

fc /.oo / \/-EnEE\rEr'EE i \/-E r^EH-frH r^HE \ 

E = -— dkTi ' ' ' ' (7.11) 

2vr In \ i-\/HryHE-[/EryEH , i/H ryHHi/H ryHH J 

-'U \ (-^012^^2 ^021 + ^1 '-^012 ^2 ^^021 ) 

We identify each term with this obvious notation: 

-E' = Eee + Eeh + E^e + Ehh. (7.12) 



7.2.1 Purely electric and purely magnetic Casimir energy 

The purely electric case has been solved by K. A. Milton et. al. in [27]. 
Using Eq. (7.113) and Eq. (7.12), we have to solve: 

Eee = -^f dkTr {VfG^gV.^G^i) . (7.13) 

Replacing each potential by its value and assuming isotropy, it is followed that they 
are proportional to the identity matrix. Then we can drop them for the coordinates 
trace, but not for the spatial positions trace: 

Eee = -ei62^ dk J dv, j rfr2XiX2Tr [G^^^G^i) . (7.14) 
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Using Eq. (7.113) of the Appendix and R = |r — r'| we have, in matrix form: 



G^^{R, k) = -{3 + 3kR + eR^) 



A-kR^ 
I 



^ Ri R1R2 R1R3 ^ 



R2R1 R2 R2R: 



+ {l + kR + k'^R^) 



-kR 



v 



-R3-R1 R3R2 R3 J 

1 ^ 

1 
1 



(7.15) 



So the trace in spatial coordinates is easily solved as 



Tr {G^iG^i) = (6 + l2kR + IQk^R^ + Ak^R^ + 2k^R^) 



,-2fc_R 



(47r)2i?6- 



(7.16) 



We can simplify our calculus using the dimensionless variable u = kR to factorize 
spatial and frequency contributions of Eq. (7.14) 



Eee = — Ci^ 



(7.17) 

Finally, using Jg°° duu'^e~"'^ = -s^, the following result is obtained: 



E 



EE 



-23eie2- 



hc 



dvidv2 



(47r)3 j^j^ |ri -ral^' 



(7.18) 



We obtain the same result as in [27] because the dyadic Green function used in [27] is 
the pure electric part of the Green function matrix we use here. 

Having into account that the pure magnetic part of the Green function matrix is 
equal to the pure electric part G^^fi^R^k) = GQj^{R,k), the calculus of the purely 
magnetic part of the Casimir energy is similar to the electric one. Consequently we 
obtain the following result for the purely magnetic part of the Casimir energy: 



Ehh = — 23/ii/i2 



he 



dridr2 



(47r)3 72 |ri -r2r 



(7.19) 



To our knowledge, this is the first place where this result is obtained. 
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7.2.2 Coupled electromagnetic Casimir energy 



Here we are going to calculate the contribution of the electromagnetic coupling part of 
the Casimir energy. Using Eq. (7.114) and Eq. (7.12), we have to solve: 

he 



E 



EH 



27r 



(7.20) 



(7.21) 



Replacing each potential by its value and assuming isotropy, it is followed that they 
are proportional to the identity matrix. Then we can drop them for the coordinates 
trace, but not for the spatial positions trace: 

Eeh = -^1/^2^ Jo I I ^''2XiX2Tr (G^if G^f ) . 

Using Eq. (7.114) of the Appendix and R = |r — r'| we have, in matrix form: 

R3 

-R3 

y i?2 —Rl 

SO the trace in spatial coordinates is easily solved: 



Ggf{R, k) = -k-^e,,AR = Go{R, k)^[k + 



k 
R 



-R.\ 



Ri 




(7.22) 



TT{Gg",Ggf)=-2eGl{R,k)(k 



R 



(7.23) 



We can simplify our calculus using the dimensionless variable u = kR to factorize 
spatial and frequency contributions of the formula 

Eeh = ei/i2^ J dr, J dv,^ due~'^ {u' + 2u' + u') . 

Finally, using duu^e~"''^ = -Jttt, the final result is obtained as 



(7.24) 



Eeh = 7eiyU2- 



hc 



dvi / dr2 



X1X2 



(7.25) 



(477)3 J ^ J - J17 

There is an antisymmetry between G^^ and G^^ shown in Eq. (7.114) and Eq. (7.115) 
of the Appendix, that is GQ^{R,k) = — G^^(-R, fc). Therefore, we can obtain the 
coupling between the magnetic part of the first object and the electric part of the 
second one in a similar way. Then, we can obtain the coupling between the magnetic 
part of the first object and the electric part of the second one: 

he f f dYidY2 



E 



he 



7e2/ii 



(4vr)3 JJ^ \v1-Y2V' 



(7.26) 
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These two cross contributions to the Casimir energy in the soft dielectric and diamag- 
netic hmits are the main resuh of this Chapter. It is interesting to note the sign change 
in this part of the Casimir energy with respect to the other contributions. Then, these 
terms can invert the typical attractive nature of the Casimir energy to repulsive in this 
limit for objects with very different electromagnetic nature. Finally, thanks to these 
new two terms, we obtain the complete electromagnetic Casimir energy between two 
objects in the soft dielectric and diamagnetic limit for the first time. Then the global 
asymptotic electromagnetic Casimir energy is finally obtained: 

B = ^0 / / (7.27) 

where 7 = 23eie2 — 7ei/i2 — 7e2/ii + 23/ii/i2 is the complete multiplicative constant, 
taking into account all the electromagnetic effects in this limit, nor just electric effects, 
where 7 would be 23eie2. 

7.2.3 Asymptotic Casimir energy 

We can also obtain from Eq. (7.27) (always in the soft electric and magnetic limit) the 
first asymptotic energy order for two objects. These formulas are the first order of a 
multipolar expansion of the integrand in Eq. (7.27) in the coordinate system of each 
object. That means that we have to assume that the distance between the objects is 
much greater than their characteristic lengths Ra, that is Ra <^ R- In this limit we 
can separate the problem into two scales and we can replace |ri — r2| = R, where R is 
assumed to be a constant. Then the integral can be easily solved in this limit to: 

where Va is the volume of the a - object. As it is demonstrated in the Appendix B 7.10, 
in the soft limit order we can approximate the electric and magnetic polarizabilities as 
= and = jl^. In [137] it was derived for an sphere. By using the same 
method with the approximation of that the effective field over the dielectric is equal 
to the induced field in the soft limit, we arrive at ~ arbitrary shaped 

object. Equation (7.27) simplifies into 

E ~ (23Qf af + 23af af - 7af af - 7af af ) . (7.29) 
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This is the soft response hmit of the Feinberg and Sucher potential [117]. It coincide 
as well with the soft response limit of the asymptotic Casimir energy obtained in [21] 
for spheres. The presented scheme looks as we would use a local two point potential 
(where we substitute the polarizabilities of the bodies by their local susceptibilities) as 
in the present case with the potential given in Eq. (7.27), and integrate over the volume 
of each body to obtain the PSA of the Casimir energy. If we also study the asymptotic 
distance limit of the PSA, then we reobtain the Feinberg and Sucher potential in the 
soft response limit, now proportional to the polarizabilities of the objects instead the 
susceptibilities. 



7.3 Diluted limit at any temperature T 

In this Section we will calculate the Casimir energy in the diluted limit at any finite 
temperature. Then we will focus on different approximations to low an high tempera- 
ture limits, recovering the zero temperature case showed before and giving new formulas 
of the PSA for any temperature. We begin this study with the Casimir energy formula 
(7.1) for any temperature: 

oo ' 

E = A;bT^ log|l-N(fc„)|, (7.30) 

n=0 

with Matsubara frequencies kn = 2Tx^^n = An. The tilde means that the n = case is 
weighted by a 1/2 factor. As usual, we apply the Born approximation to the T matrix 
scattering obtaining at first order = Va (where a labels the body in interaction). We 
apply again that log |y4| = Trlog(y4) and log(l — x) = — J2'^=i ^ 5 then we transform 
Eq. (7.30) into: 

00 ^ 00 ' 

E = -fcBT^-^Tr(NP(An)). (7.31) 

p=l P n=0 

At first order in p we obtain: 

00 ' 00 ' 

E = -kBTj2 Tr (N(An)) = -^bT^ Tr (T1U12T2U21) . (7.32) 

n=0 n=0 

As we are working in the positions representation space instead in the multipolar 
representation space, we have to represents the operators Tq and in the positions 
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space. For that issue, we have into account that Uq/3 matrices represent the field 
propagation between two points. Then they are represented by the free vacuum Green 
function of the interaction field. Using the Born approximation, we obtain the next 
formula for the diluted approximation of the Casimir energy: 

oo ' 

Ec^Et = -ksTY, Tr ( l^i Go 1 21^2^021) • (7.33) 

n=0 

As we did with (7.5), we can separate the magnetic and electric part of each body 
potential and each field contribution in the free dyadic Green function. We trace over 
the EH space obtaining: 

As before, we identify each term with this obvious notation: 

Et = Eee + Eeh + Ehe + Ehh- (7.35) 

7.3.1 Purely electric and purely magnetic energy 

Using Eq. (7.113) and Eq. (7.116) of the Appendix, we can solve the trace over the 
coordinates of Eee and, similarly, of Ehh- The matricial form of the purely electric 
part of the dyadic Green function and of the purely magnetic part of the dyadic Green 
function are: 

G^if (i?, k) = -RiRj (3 + 3kR + en^) + 5ij [l + kR + k^R^) |^ 

= Ggf{R,k). (7.36) 
Where R = — r^. The trace in spatial coordinates is easily solved obtaining: 

oo ' „ „ —2kR 

Eee = -kBTj2 ^1^2 / dvi / dr2-^^-^ [Q + l2kR + lOk^R^ + Ak^R^ + 2k^R^] . 

n=0 ^ ' 

(7.37) 

Replacing k by An we obtain: 

F _ 7 2 [ f dridr2 ^' 2ARn^^ 



n=0 



X 



[6 + 12ARn + lOA^R'^n'^ + AA^R^n^ + 2A^R%^] . (7.38) 
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Having into account that J2'^=o e = jz^? is easily deduced: 

oo oo 



(7.39) 



n=0 n=0 

Denoting A = -RA, we can carry out the sum obtaining: 



E 



EE 



ksT f f dvidv2 1 

(4^"'"' Ji '2 (e2A- 1)5 



X 



[e2A (_9 _ + lOA^ - + 2A^) + e^^ (6 + 36A - lOA^ - 12A3 + 22A^) 
+e^^ (6 - 36A - lOA^ + 12A3 + 22A^) + e^^ (-9 + 12A + lOA^ + + 2A^) 
+3e^°^ + 3]. (7.40) 

We obtain a similar result for the purely magnetic Casimir energy replacing €162 by 

/il/i2. 



7.3.2 Coupled magnetic - electric Casimir energy terms 

We perform the same calculations using (7.114) and (7.115) of the Appendix as for the 
purely electric case. The trace over spatial coordinates is already done in (7.23), so 
we obtain the formula for the coupling terms of the Casimir energy just replacing the 
integral in (7.24) by a sum: 

dridr2' 



e-'"'' {R'k^ + ^R'e + R'e) . (7.41) 



n=0 



Here fc„ = 27r^^n = An are the Matsubara frequencies. After solving this sum and 
denoting A = RA, we obtain: 

2kBT^ ^ r [ dYidY2 A^ 



E 



EH 



X 



'1 J2 



r^X 

-1)5 

e4^(-l -6A + 11A2) 



e'^ (1 - 2A + A 
+6^^ (-1 + 6A + 11A2) + e^^ (1 + 2A + A^ 



(7.42) 



Replacing 61/12 by /iie2, we obtain the formula for Ehe- These formulas for the Casimir 
energy in the diluted limit are valid for any temperature, but they are too much 
complicate for analytical analysis at any temperature. Therefore, we study the limits 
at high and low temperatures. 
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7.3.3 Low and zero temperature limit 

Here we will recover the zero temperature limit of the diluted limit. Then we will make 
a perturbative analysis valid for low temperatures. We can make it easily by expanding 
the Taylor series in A = RK of the integrand and studying just the first orders deleting 
the rest ones. By using (7.40) and (7.42), we reobtain (7.27) for the zero temperature 
case. We need to take the fifth order series term of the Taylor expansion of Eq. (7.40) 
and Eq. (7.42) to get the next non zero perturbation term of the Casimir energy as 



because Aii? = A2E = A^E = A^E = 0, where AnE is the n-th order correction term 
to the low temperature expansion. 

7.3.4 High temperature and classical limit 

In this Section we will obtain the high temperature limit of the Casimir energy in the 
diluted limit, whose first term will be the classical limit of the Casimir energy. To 
obtain this limit, instead of solving the sum in Eq. (7.37) and Eq. (7.41), we will just 
keep the first term of the sum. Then the classical limit of the Casimir energy in the 
diluted limit reads the first sum term: 




(7.43) 




(7.44) 




(7.45) 



pel 



'EH 



Ei 



HE 



(7.46) 



0. 



And the first perturbation to that limit 



is the next sum term: 





(7.47) 



6162 



yUlAi2 
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= = ^ [ [ ^e-^^^ (R'A' + 2R'A' + R'A') . (7.48) 

ei/i2 /iie2 (47r)2 JJ2 ^ ! y > 

Where AiE'^^ is the first correction to the classic hmit of the Casimir energy E'^K With 
these results we see that we loose the coupling between electric and magnetic energy- 
terms in the classical limit. 

Finally we must remark that, if we take the asymptotic distance approximation as 
in Sect. 7.2.3, we reobtain the results given in [115] in the soft response limit. 



7.4 PSA for three bodies system 



In this Section we are going to calculate the Casimir energy between three bodies in the 
Pairwise Summation approximation. In this asymptotic limit we will obtain that the 
energy of the system will be the addition of the PSA energy of each pair of objects. This 
linear behavior of the Casimir energy was expected in that approximation although we 
know that it is in general false [130]. We begin this study from the Casimir Energy 
formula for three objects given in [122] [130]: 



E 



he 
2^ 



dk log 



(7.49) 



The M matrix (whose coefficients are non commutative matrices) is: 





-U12 


-Ui3 


-U21 


T2' 


-U23 


-U31 


-U32 





v 

So, using the logarithm product rule, the Casimir energy between three bodies is 



-E3 



he 
2^ 



dk\og\l -N12I + 



he 



dk\og\l -N13I + 



he 
2^ 



27r 

Here the N matrix is Nq/j = Ta\iai3^i3^pa, and: 

M = (1 - N13)-' (T3U32 + N31) (1 - N12)-' (T2U23 + N21) . 



dk\og\l - M| . (7.50) 



(7.51) 
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Expanding the inverses we obtain the double series 

oo oo 

R=Y, (^13)"' (T3U32 + N31) (^12)"' (T2U23 + N21) . (7.52) 

ni=0 n2=0 

The lowest order of the Born series of the M matrix will come from the first order 
expansion series making rii = n2 = 0, that is 

M ~ (T3U32 + N31) (T2U23 + N21) . (7.53) 

In addition to that, we just consider the sum term with the minimum number of T 
matrices products, because that will be the lowest order expansion in susceptibilities 
of the M matrix. That means that we reduce the highly non linear M matrix to 

M - (T3U32T2U23) = N23. (7.54) 

Replacing Eq. (7.54) in Eq. (7.50), we obtain the next PSA of the Casimir energy 

between three objects: 

he /"°° he C°° he C°° 

E^ = — dk\og\l-ni2\ + — rffclog|l -N13I + — / c/A;log|l -N23I . 



2ti Jq 2tt Jq 27T jQ 

(7.55) 

Therefore we obtain that the PSA approximation of the Casimir energy between three 
objects is the sum of the PSA energy of each pair of objects. In other words, we rederive 
a kind of superposition law of energies in the diluted PSA limit as expected, because 
the usual presentation of the PSA approximation is the assumption that we have a 
superposition behavior in that asymptotic limit. Here we have justified the validity of 
that approximation. 

If we take the asymptotic distance limit to Eq. (7.55), we will obtain a superposition 
of two bodies PSA energies in the asymptotic limit. The nonlinearity of the Casimir 
energy must be given by higher orders expansion terms even for systems with three 
bodies. The same superposition behavior is expected for the general N body case, and 
it will be proven in the next Section. 

7.4.1 PSA for general N bodies system 

In this Section we are going to generalize the PSA energy of three bodies given by 
Eq. (7.55). For this purpose we are going to use an iterative procedure which will give 
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us the new terms to include to the PSA of the Casimir energy of — 1 bodies when we 
include another new object to our system. Let us represent the M matrix of Eq. (7.49) 
as the sum of its diagonal and its non-diagonal parts as [122]: 



or symbolically as 



+U, 



(7.56) 



(7.57) 



then it is easy to find that the inverse of M is the next perturbative series 

oo 

M-^ = T^(-UT)". 



(7.5J 



n=0 



On the other hand, the M matrix of the N objects system is related with the 
of the — 1 objects system by block matrices in the next way 



matrix 



rrr—l 



Where 7 index goes from 1 to — 1. With this result we calculate |Mjv| obtaining 



l,7l 



(7.59) 



In the PSA approximation is valid the substitution M^^ ~ T, where T is a diagonal 
matrix whose — 1 diagonal elements are with 7 index defined as before. Then, 
using this approximation we can approximate the M^r determinant as: 



Hat 



_i||TAr|-i|l-T^U-,.^_iTU^_.i 



,71- 



(7.60) 



Where we multiply by blocks the matrix of the last determinant obtaining 



NiAT ... n \ 



N 



2N 












(7.61) 



7.4 PSA FOR THREE BODIES SYSTEM 



125 



Finally, we obtain the desired result: 

N-l 

\Mn\ = |Miv-i||Tjvr^ n 1^ ~ ^^^1- C^-^^) 

7=1 

The other needed matrix is 

\M^,n\ = \M^,n-i\\Tn\~^- (7.63) 
And with the initial conditions: 

|Mi| = |Moo,i| = \Ti\-\ (7.64) 
we obtain these determinants in a closed form as 

N N l-l 

\^n\ = n n n ii - ^^-i c^-^s) 

k=l 1=2 m=l 

and 

N 

iMoo,ivi = n I'^'^r'- C^-se) 

k=l 

Taking the logarithms of Eq. (7.49), the Casimir energy for the N bodies system is 
approximated as 

oo N l-l 

En = ^ dkJ2J2^og\l-Nim\ (7.67) 

"^0 1=2 m=l 

in the PSA limit. Eq. (7.67) is the generalization of Eq. (7.55) for the body case, 
and it show us that in the PSA limit, we always obtain a superposition principle of the 
two body PSA Casimir energy despite the fact that Casimir energy is not a nonadittive 
interaction. We will obtain the non-linear effects in the next orders of the expansion. 

This result is qualitatively different of the usual PSA procedure of integrating each 
point of a point potential over the volume of each object. The reason is simple, 
a point potential is proportional to polarizabilities [138], but in the soft limit 
approximation, the lower allowed term is proportional to two polarizabilities. Therefore 
we should study the (A^ — 1) expansion term to obtain the first one proportional to A^ 
polarizabilities. So in a consistent calculation of PSA energies for three or more objects, 
the asymptotic approximation of the PSA energy in the soft material approximation 
is different of the A^ point potential function. In the next Section we will find where 
the contribution of the three point potential is relevant in the diluted limit. 
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7.5 Second order expansion of diluted limit 

In this Section we are going to study the second order expansion of the PSA in the 
diluted hmit. This is a comphcated long calculus so, instead the complete electromag- 
netic case, we will restrict ourselves to the purely electric case. That means that we 
will ignore the magnetic properties of the bodies. It could be interesting to show the 
complete study of that series term, but it is a long calculation to show here. However, 
it possesses theoretical utility, as it shows the lack of the superposition principle in 
the Casimir energy calculations. The next results have been done with the help of 
Mathematica [139]. 

Again we start from the Casimir energy between two compact objects as given by 
Eq. (7.2). We studied the first order expansion of the energy in Sect. 7.2. Now we 
are interested in the next order expansion term in polarizabilities, so we maintain the 
study of the fist term of Eq. (7.2) taking just the case p = 1. In addition to that, we 
take the second order approximation of the Born series of the T matrix: 



Applying the linearity of the trace and taking our attention just to the second order 
expansion in polarizabilities (that is, to the third order term in polarizabilities), we 
obtain the following result for the second order of the diluted limit result: 



E2 = -— dkTi {VrGoriViUi2V2U2r) - — / dkTi (\/if/i2'V2'G'o2'2V^2f/2i) • 



We will center our study just on the first integral of (7.69), because the analysis of 
these two integrals is the same. First we make the trace over the EH space, which 
is automatic here because we have canceled the magnetic properties of the bodies. In 
other words, we are not studying the coupling between electric and magnetic induced 
dipoles. After that we trace over the space coordinates, so we need the matricial form 
of the electric part of the dyadic Green function given in Eq. (7.113) of the Appendix, 
where R = r — r'. We will need to define the function: 



n[a, (3Uk) = Rip (1 + kR^p + k'^Rlp) + (3 + 'ikR^p + k'^Rl^) R^f,. (7.70) 



oo 




(7.68) 



n=0 





(7.69) 
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Where k is the frequency, a and j3 labels the bodies coordinates and i shows the vector 
component of Ilai3 used. So we can solve the trace over the space coordinates obtaining 
the following result for the first term of (7.69): 

+n[i,2]yik)n[2,i']yik)n[i',i]yik) 

+n[i, 2Uk)n[2, i'Uk)n[i', lUk)]. (7.71) 

Here the space integrations are performed just over the body volumes because it is 
there where dielectric and diamagnetic potentials are defined. It is possible to perform 
the k integral, because we have an exponential multiplied by a polynomial integrated 
between zero and infinity. It is not shown here because it contains 216 sum terms. 

For any finite temperature we also obtain a long result, but in the classical limit we 
can at least obtain a tractable formula. This limit consists of keeping only the k = 
mode. In that case we obtain the result 

; keT f f f 616261/ 



1287r3 Ji' RiiRli'Rl'i 

X [7^[l, 2],.(o)7^[2, i%(o)7e[i', i],(o) 
+7^[l,2],(o)7^[2,l'],(o)7^[l^l],(o) 
+7^[l,2],(o)7^[2,l%(o)7^[l',l],(o)]. (7.72) 



Where TZ[a, /3]i(0) = (i?^^ + 3i?^^j), so we can write: 

Et = [drJdrJ drr '''''' 



X [[Rl^ + 3i?i2^.) (i?2i/ + 3i?2i'x) [Ryi + '^Ri'ix) 

+ (-^12 + 3-R?23/) [RIv + '^Rlvy) {Rl'l + 3-Rl'lj^) 

+ (-R12 + ^Ri2z) {RIv + ^Rli'z) {Rl'l + ^Ri'iz)]- (7.73) 

When performing this integral, a term proportional to R^,^ and another one to R^,\ 
appear. These terms can be problematic because i?i/i = belongs to the integration 
interval for all the points of the 1 body. But these singularities can be removed with the 
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appropriate regularization procedure. Note that if we make 1' = 3 in Eq. (7.69), this 
integral would also be one of the contributions of the second term expansion of a three 
objects system. In that case we have not got any problem in the volumes integration 
for this term and this expansion term looks quite similar to the local analog of the three 
points potential given in [138]. In fact, in the three (or more) bodies system that term 
is the first term that breaks the superposition behavior founded in 7.4. Replacing the 
term by the operator —A in Eqs. (7.113) and (7.116) of the Appendix, and replacing 
the k term by the operator \/—A in Eqs. (7.115) and (7.114) of the Appendix, which is 
a change valid on shell because Go{R, k) = ^^^^ we recover the formalism of Power and 
Thirunamachandran [138] in the soft dielectric and soft diamagnetic limit. This scheme 
let us take into account not only electric phenomena as in [138], but also magnetic and 
electro-magnetic coupling effects so as an extension to finite temperature cases. Taking 
this argument into account, we can understand the structure of the perturbation terms 
of the PSA limit of the electromagnetic Casimir energy between bodies. The Casimir 
energy in the PSA limit has the structure of a series of infinity terms whose ra-order 
term is the sum of integrals over n + 2 bodies (which can be repeated or not) of the 
{n + 2) points local em potential. Each series term come from the expansions made in 
Eq. (7.2) and in Eq. (7.4) and is proportional to the product of n + 2 permittivities. 
If we take the asymptotic distance limit of these integrals, we will obtain a series of n 
points EM potential now with polarizabilities instead susceptibilities as in [138], but 
in the soft response limit. When we use perturbations of the tree term of the Born 
expansion, it appears divergent terms in the integrand as seen in Eq. (7.73). These 
singularities can be removed with the appropriate regularization procedure [140]. 



7.6 PSA for objects with general magneto— electric 
coupling 

In this Section, we extend the study of PSA to materials with a general magneto- 
electric linear coupling with the electromagnetic field, we use the potential Vi given in 
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Eq. (7.10) to obtain the PSA Casimir energy of frequency independent dielectric with 
magnetoelectric couphngs at zero temperature as 

where 

7o = 23eie2 - 7ei/i2 - Tp-ih + 23/ii/i2 + 7aia2 + 23ai/32 + 23/3ia2 + 7/3i/32- (7.75) 

This result extend the Feinberg and Sucher potential [117] to objects with magneto- 
electric couplings, but only in the diluted limit. In the high temperature limit, the 
PSA of Casimir energy is 

-keT f f dridr2 , . 

(47r)^ J^J2 \ri-r2\^ 

where 

7d = 3eie2 + 3/ii/i2 + 3ai/32 + 3/3ia2- (7.77) 

It is also possible to obtain results valid at any fixed temperature from Eq. (7.30), 
the results are Eq. (7.40), but with (q;i/32 + /3ia2) instead €162, and Eq. (7.42), but with 
(— aia2 — instead ei/i2. Then the correction terms to high and low temperature 

cases are given by Eqs. (7.43), (7.44), (7.47) and (7.48), being careful with the sign for 
terms proportional to aia2 and /3i/32- 

As we can see, if magnetoelectric couplings have different signs and are strong 
enough to compensate the usual electric and magnetic coupling of dielectrics, a repul- 
sive Casimir energy between objects can be achieved in the diluted limit. We will study 
the appearance of this repulsion in a more physical system of Topological Insulators in 
the next Section. 



7.7 PSA for Topological Insulators 

Recently, in [141], A. Grushin and A. Cortijo demonstrated the existence of an equilib- 
rium distance between topological insulators with opposite topological polarizabilities 
sign between parallel plates. This result was extended for all temperatures in [133], 
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leading to a reduction of the equilibrium distance with the temperature until the dis- 
appearance of the equilibrium distance in the high temperature regime. Regions of 
Casimir repulsion for all distances were also reported for some values of w = ^ and |^| 
(where w is the ratio between the oscillator strength Ue and the resonant frequency 
of the material, while 6 is the magnetoelectric topological susceptibility, their definition 
is given below). In this part of the article, we apply the formalism developed in the 
previous Section to TI. 

Following [141] and [133], the electromagnetic response of three-dimensional topo- 
logical insulators, which determines the reflection coefficients, is governed by the La- 
grangian 

£ = £o + = E ■ D + B ■ H. (7.78) 

The first term of Eq. (7.78) reads Cq = ^ (eE — /i^^B), and it is the usual electro- 
magnetic term, which gives the usual constitutive relations of dielectrics without magneto- 
electric couplings. The second term of Eq. (7.78) is a nontrivial axionic [142] or topolog- 
ical magnetoelectric term which gives rise to the topological magnetoelectric [143] [144] 
term in the Lagrangian £e = If E ■ B, where a is the fine structure constant (a = |^) 
and 6 is the magnetoelectric topological susceptibility. Therefore, the constitutive re- 
lations of TI are easily obtained from the Lagrangian of TI given in Eq. (7.78) as 

D = eE + a (^-^ B, 

H = /i-^B - a (^^^ E. (7.79) 

Time-reversal symmetry indicates that ^ = 0, '7r(mod 2it) being 6 = it the case for TI 
and ^ = the case for trivial insulators. When the boundary of the considered object 
is included, the action corresponding to the Lagrangian Eq. (7.78), Sq + Sq, is a fair 
description of the TI only when a time-reversal breaking perturbation is induced on the 
surface to gap the surface states [143]. In this work we consider that the time-reversal 
perturbation is a magnetic coating of small thickness / compared with the rest scales 
of the problem (therefore we will neglect its contribution to Casimir effect) which gaps 
the surface states. In the described situation, 9 is quantized in odd integer values of vr 
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such that 

e = {2n + l)n, (7.80) 

where n G Z, determined by the nature of the magnetic coating, but independent of 
the absolute value of the magnetization of the coating. Positive or negative values of 
6 are related to different signs of the magnetization on the surface [145], which we 
consider is perpendicular to the surface of the body. Being a topological contribution, 
6 is defined in the bulk as a constant whenever the bulk Brillouin zone is defined [144]. 
Following the discussion in [133], we assume that the frequency dependent dielectric 
function e{u) is described by an oscillator model of the form: 

e{in) = eo + V , ^ ^ , (7.81) 

In this model, w/j j is the resonant frequency of the ith oscillator, while Wg.i accounts 
for the oscillator strength. The damping parameter 7/j j satisfies 7ij^j <^ ^R,i^ playing 
therefore a secondary role. As there are few experimental results of e(m) for TI, we will 
assume that there is just one resonance and do not take into account the contribution 
of the damping constant, therefore the dielectric function considered here has the form: 

e{iK) = eo + -—^. (7.82) 
1 + 

where w = In addition to that, we assume fi = fiQ. 

In order to apply the formalism developed in the last Section, we need the fields D 
and B function of E and H, then 

2^ 




D = 

B = fia (^^^ E + fiH, (7.83) 

2 

where a is the fine structure constant (a = j-). Note that the PSA can be widely 
applied to TI, because it is expected that the magnetoelectric topological susceptibil- 
ity 6 was small (and the magnetoelectric coupling proportional to the fine structure 
constant), and the interesting behavior of TI because of the Casimir effect appears also 
for the soft dielectric response of the TI. 



132 



Pairwise Summation Approximation of Casimir energy 



The case of multiple resonances or any other more exact model of e{iK) [141] and 
magnetoelectric coupling (which could have into account the contribution of the surface 
of the TI to the magnetoelectric response) can be easily performed within the formalism 
presented here. 

Because of the frequency dependence of e(m) for TI, it is not possible to obtain the 
integrand of the PSA energy at zero temperature in a closed form in terms of simple 
functions. Then the PSA Casimir energy at zero temperature for topological insulators 
is given by the integral 

he f f dridr2 



PSA 



(47r)3 7i J2 |ri - rsl 



f7o, 



(7.84) 



where 70 now depends on the distance between points of the TI and on the dielectric 
properties of both TI as 

7o = wlwlxfi{x) + 60aia2 + 2'ia\ci\ + [a^wl + ci\wf) xf2{x), (7.85) 



where aj = a — , w = x 
/2(x) the next functions: 



c 



|ri — is the adimensional distance and fi{x) and 



/i(^) 



6xGi 3 



0.4,1 



sets I +5GI 



X 



2 ' ^' 2 



3r^3,i 



+ 2x^01 



X 



(7.86) 



( 



^xGl 



3,1 



X 



^xG\ 



3,1 





i, 1 



3,1 



^3,1 



+ 3Gt3 
) + x'Gl 



3,1 



X 



X 



^3,1 
,3 



where G^" 



/ 



(7.87) 



are the Meijer G functions. Eqs. (7.86) and (7.87) and their 
limits at short and large x are plotted in Fig. 7.1. All coefficients of 70 are positive 
{xj\{x), xf2{x), 23 and 60, because /i(x) and /2(x) are smooth positive functions in 
the real positive axis, as seen in Fig. 7.1) for any value of 61 and 62- 
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2 4 6 8 X 10 

Figure 7.1: Representation of given in Eq. (7.86) (red color) and f2ix) given in 

Eq. (7.87) (black color). The blue curve is their common limit at large x, the green curve is 
the limit of f2{x) at small x and the yellow curve is the limit at small x of /i(x). 



At short distances we can approximate fi{x) ~ 37r and f2{x) ~ Gvr, then 70 can be 
approximated by 

70 = eOaiOa + 23a?a^. (7.88) 
If sign(^^i) = — sign(02), this energy becomes positive (i.e. repulsive) for all 

m<-'^,. (7.89) 

For diluted materials, this condition is beyond the positive energy condition \6i62\ < 



^Vei(0)e2(0) < [l-iS], then we can consider it is always fulfilled for diluted 

dielectrics. 

At large distances we can approximate fi{x) ~ ^ and f2{x) ~ ^, then 70 can be 
approximated by 

7o = wlwl + dial + d\w\ + dlwl + — 0.10.2- (7.90) 
If sign(6'i) = — sign(6'2) and wi = W2 = w, the energy becomes positive (i.e. repulsive) 
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for all 



2 23/2 -2 n ~A 240 _ _ 46_„ 2 
^ ^ 46 I - "2 + Y "1 + "2 - "^"i"2 - ^«l«2 



(7.91) 



As a consequence, there are several different regimes for the system in the quan- 
tum limit. When sign(6'i) = sign(6'2), the Casimir energy is enlarged because of the 
contribution of topological charge 6, but when sign(^i) = — sign(^2); different regimes 
appear. 

When sign(^^i) = — sign(^^2) and for low enough absolute values of 61 and 62, the 
condition given in Eq. (7.89) is fulfilled while the condition given in Eq. (7.91) is not, 
then we have a repulsive Casimir energy at short distances and an attractive Casimir 
energy at large distances, then there must exist an stable equilibrium distance in this 
case. 

But if the condition given in Eq. (7.91) is also fulfilled, the magnitude of the posi- 
tive arguments of Eq. (7.85) is not large enough to compensate the repulsion because 
topological charges 6, then we obtain repulsion for all distances. 

In the high temperature limit, another result is obtained. In this case the PSA 
gives the Casimir energy as 



where, contrary to the zero temperature case, "fd does not depend on the distance 
between points, and it is given by 



Depending on the values of w and \6\, we obtain either repulsion for all distances or 
attraction for all distances. The condition to be fulfilled in order to obtain repulsion 
in the classical limit when wi = W2 = w is 



In Fig. 7.2, we represent the behavior of the PSA energy as a function of 6 and w 
when Wi = W2 and 61 = 6 = —62 in the quantum and in the classical limit, having into 
account that there is a forbidden region because of the positive energy condition [146]. 




(7.92) 



QQ 

7d = W1W2 + aiW2 + (^2^!! + «i«2 + 2aia2- 



(7.93) 




(7.94) 
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In both cases we find different regions of parameters of repulsion for all distances, but 
in the quantum limit there is a region of existence of an equilibrium distance, while in 
the classical limit a region of attraction for all distances appears. 




TT 



Figure 7.2: Attraction versus Repulsion in the classical hmit (T — t- oo) and in the 
quantum hmit (T — t- 0) as a function of w and \6\. For the Classical limit, the green region 
shows repulsion for all distances and the rest of phase space shows attraction for all distances. 
For the quantum limit, the red and green regions show repulsion for all distances and the white 
region shows parameters for which an equilibrium distance appears. Several curves of constant 
adimensional equilibrium distance deq = ^|ri — are also plotted. The black region is a 
forbidden region for the parameters due to the positive energy condition |a| < \/ e(0) [146]. 

It is expected that at finite temperatures the equilibrium distance found in the 
quantum limit would be reduced when temperature increases, and Fig. 7.2 suggests 
that behavior, but we have let this study as a future work. 
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7.7.1 Casimir energy between spheres and plates 

A claim would be done to the method performed here. It is a method just valid in the 
diluted limit. Then the results presented here are general results in the geometry of 
the bodies, but not in the dielectric response of materials. 

In this Section we will use the PSA formulas presented in the previous Chapter 
to obtain the Casimir energy to specific systems. In particular, we study the sphere- 
plate system because of its experimental relevance and the two infinite plates system 
to compare the PSA results with the exact ones. Both systems were studied at zero 
temperature. 

We were not able to obtain analytical results, so a numerical method has been 
implemented. 

The two infinite parallel TI plates have been already studied in [141]. Our interest 
here is to compare the results obtained by the PSA with the exact ones. We impose 
6*1 = —62 = TT in order to obtain an equilibrium distance, and different w from w = 0.2 
to w = 0.6. The obtained results of PSA and the exact ones are compared in Fig. 7.3, 
where we find that PSA tends to overestimate the magnitude of the Casimir energy 
the more the larger w. This fact reflects the nature of the approximation made in 
PSA, where we must assume the diluted limit is valid. On the other side, an excellent 
approximation of the equilibrium distance between plates is obtained, then the PSA 
gives a good qualitative result. 

We have also used the PSA to study the sphere-plate system because of its ex- 
perimental relevance in Casimir effect experiments. In this case, we also impose 
61 = —^2 = TT and w = 0.45 for both plate and sphere in order to obtain an equi- 
librium distance at T = (see Fig. 7.2). We vary the dimensionless radius of the 
sphere Rg = from Rs ^ until = 1 to observe its effect in the equilibrium 

distance. In Fig. 7.4, the PSA energy per unit of volume of the sphere is plotted as a 
function of the distance between the nearest points of the plate and the sphere. As a 
result, an equilibrium distance has been obtained for all the studied sphere radii. This 
equilibrium distance reduces when the radius increases until reaching a constant value. 
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Figure 7.3: Comparison between exact Casimir energies in units of £"0 = 
curves) and PSA (points) at zero temperature for TI infinity plates 
with 9i = —O2 = TT and different e(0) as a function of the adimensional distance d [141]. PSA 
tends to overestimate the absolute value of the energy, but captures the equilibrium distance 
without appreciable error. The overestimation of the energy increases with e(0), as expected 
because of the nature of the approximation. 



The appearance of an equilibrium distance is easy to understand because of the nature 
of the integrand in PSA calculations at zero temperature. As discussed in Eq. (7.85), 
from a given distance, the points of the sphere nearer to the plate tend to increase the 
Casimir energy (giving a repulsive contribution), while the rest of points of the sphere 
tend to reduce this energy (giving an attractive contribution). As at short distances 
Eq. (7.85) diverges as ^ and ^ diverges, it is evident that there is always a distance 
where repulsion compensates attraction, reaching the system an equilibrium distance. 



138 



Pairwise Summation Approximation of Casimir energy 



10^ 



100 



0.01 








\ 


-m-Rs 


= 10 


-3 


-•-Rs 


= 10" 


-2 


—m—Rs 


= 10 


-1 


[^Rs 


= 1 


J 




0.002 



0.004 



h 



0.006 



0.008 



0.010 



Figure 7.4: Dimensionless PSA Casimir energies in the sphere-plate Tl system per unit 



of volume of the sphere e 



with Q\ 



vr and w = 0.45 as a function of the 



dimensionless distance between the plate and the contact point of the sphere h = d — Rg, 
where d is the dimensionless distance from the surface of the plate to the center of the sphere. 
Positive energies are the whole curves, while negative energies are the dashed curves. The 
equilibrium distance for each sphere radius is represented by a point. These points are joined 
by a soft line to have an idea of the dependence of the equilibrium distance with the radius 
of the sphere. 



7.8 Conclusions 

In this Chapter, we have calculated the full electromagnetic Casimir energy in the di- 
luted limit between two bodies using the formalism given in [21] reobtaining the energy 
given in [132] in the pure electric case. Although this formalism fails in the perfect 
metal case, it is valid for soft dielectric and diamagnetic bodies and for bodies with 
magnetoelectric coupling. This formalism has been generalized to finite temperatures 
cases (giving exact results at finite temperatures) and extended for an arbitrary number 
of bodies. We have shown that the Casimir energy, which is a nonadittive interaction, 
has a superposition behavior in the first expansion of the PSA and how this superpo- 
sition behavior is broken when the integration of different point potentials appears 
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in the perturbation series to contribute in the whole energy as the next perturbation 
terms of the PSA. 

This formahsm let us reobtain known results of the literature as particular cases. 
For example, high and low temperature correction terms were obtained in [115] from a 
functional formalism, here we have obtained the exact result at any finite temperature. 
In addition to that, we obtain a generalization of the N point potentials formalism 
of Power and Thirunamachandran [138] valid for diamagnetic objects and at finite 
temperatures. 

The extension of this formalism to anisotropic materials is straightforward, although 
we have not done it here, it is shown how to do it for the case where it would be needed. 

We have also used the results derived here to study the Casimir energy for a two 
infinity TI plates, to compare the PSA with the exact results given in [141]. As a result, 
we found that PSA gives better quantitative results for lower w, as expected, because 
PSA works better in the diluted limit. On the other side, we obtain an excellent 
approximation of the equilibrium distance between plates, therefore we obtain good 
qualitative results for all 9 and w studied. We have also used the PSA to study the 
Casimir energy between an infinite TI plate and a TI sphere because of the experimental 
relevance of this geometry. 

The results presented in this article are valid for objects of arbitrary shape, but in 
the diluted limit. Then we have a systematic procedure to calculate Casimir energies 
valid for bodies with soft dielectric, soft diamagnetic and soft magnetoelectric response. 

7.9 Appendix A: Obtention of the matricial Green 
function 

Here we have calculated the Casimir energy using a matricial Green function for the 
electromagnetic fields instead of the dyadic Green function used in [21]. In this Ap- 
pendix we will see that both formalisms are equivalent just using as field sources the 
induced polarization and magnetization vectors instead of the induced currents into 
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the bodies. The induced currents into the bodies are [104]: 



j = -ikP + V X M. 



(7.95) 



We introduce this hnear change of variable in the partition function, so there is not 
any new relevant term in the action of the problem. This action will be transformed 
in the next way: 

S= [ dxjGoj, (7.96) 



Where 



S 



dxA;2p^oP + / dx.V x M^qV x M 



dxikPGoV X M - / rfxV x M^o^^P 



6ij - p ViVj 



Go, 



(7.97) 



(7.98) 



is the Green dyadic function. We integrate by parts each action term obtaining: 



S 



EE 



cixPjG^^Pj, 



(7.99) 



S 



EH 



dxikPGoV X M, 



/ 



d'xikPi 



Seh = J d^Pi [ikeiap'^f'] GoM" + J dxP, 



The second term is zero: 



(7.100) 



GoM". (7.101) 



ej„^VjM"V'^Go = M ■ V X VGq = M ■ = 0. 



(7.102) 



7.9 Appendix A: Obtention of the matricial Green function 



141 



So we finally get: 



S 



EH 



(7.103) 



She is similar to Seh, but with a changed sign because of the complex conjugation of 
the induced current into the action: 



She = - / rfxV x MGoikP, 



And the last term: 



Shh = / dxV X M^oV X M, 



S 



HH 



I 



c^xe*"^VaM^ 



G'oe,a6M"V'^ (7.108) 



(7.104) 



She = J (ixM^Gj^P^' = J dxM^ [-Cg^] P'- (7.105) 



(7.106) 



GoejabVM', (7.107) 



(7.109) 

Using e'^^^Sfejab = (^"^ - (^f and VfeV^'G'o = k'^Go, we get the next result: 



= [5:^2 - V„V'"] Go. 



(7.110) 



The other term requires even a more tedious work, but it is easy to obtain that 



j^e^^^ejab^^ViV^V'Go = 0, 



(7.111) 
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because this differential operator is zero. So Shh is 

Shh = j d-KM^ [e'^'^V p5.,fi,e,,,V') M 



= j rfxMiCj^Mj-. (7.112) 

After a Wick rotation, we obtain the used form of the matricial Green function, which 
components are, using Go{R, k) = and R = |r — r'|: 

G^,f{R,k) = [fc%,+V,V;.]Go(i?,fc), (7.113) 

G^^{R,k) = -ke,,kVkGoiR,k), (7.114) 

Ggf{R,k) = ke,,kVkGo{R,k), (7.115) 

G^^f{R,k) = [fc%,+V,V;.]G'o(i?,A;). (7.116) 

Which are the results used in Eq. (7.5). It is also possible to obtain the same result 
from fluctuation - dissipation theorem, as made in [147]. 



7.10 Appendix B: Dielectric polarizability in the 
diluted limit 

In this Appendix we are going to demonstrate that, in the diluted limit, the dielectric 
polarizability is given by 

V 

«^ = e-^. (7.117) 
47r 

Following [137], the macroscopic induced polarization vector P of a given dielectric is 
defined as 

P = xeE, (7.118) 

where E is the external applied electric field and xe is the electric susceptibility of the 
dielectric, defined as e = 1 + At[Xe- Then xe = On the other side, the microscopic 
induced polarization vector p of a given dielectric is defined as 

p = a^E', (7.119) 
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where E' is the local effective electric field. This local effective electric field is defined 
in terms of the local electric field E and the macroscopic polarization field P, but in 
the diluted limit we can assume that the contribution of P is negligible, then, just in 
the diluted limit, 

E' = E. (7.120) 

Because of the definition P = ^ for homogeneous bodies, we can combine Eqs. (7.117) 
and (7.118), 

p = a^E' = VP = XeVE. (7.121) 
As we are in the diluted limit, using Eq. (7.120), we obtain 

a"" = XeV = e^, (7.122) 

which is the result we wanted to demonstrate. Exactly the same reasoning can be 
applied to obtain the diluted limit of the magnetic polarizability as = jlj-. 
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"In this house, we obey the laws of Thermodynamics!" 

Homer Simpson 



8 

Casimir energy and entropy in the 
sphere sphere Geometry 



In this Chapter, we calculate the Casimir energy and entropy for two spheres described 
by the perfect metal model, plasma model, and Drude model in the large separation 
limit. We obtain nonmonotonic behavior of the Helmholtz free energy with separation 
and temperature for the perfect metal and plasma models, leading to parameter ranges 
with negative entropy, and also nonmonotonic behavior of the entropy with temperature 
and the separation between the spheres. This nonmonotonic behavior has not been 
found for Drude model. The appearance of this anomalous behavior of the entropy is 
discussed as well as its thermodynamic consequences. 

In 1948, Casimir predicted the attraction between perfect metal parallel plates [2] 
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and between neutral polarizable atoms [1] due to quantum fluctuations of the electro- 
magnetic field. Some years later, Schwinger extended this formalism to dielectric plates 
at finite temperature [148]. Recently, a multiscattering formalism of the Casimir effect 
for the electromagnetic field has been presented [21] [22] (see also [149] and [150]). 

The Casimir effect has some peculiarities. In particular, it is a non-pairwise inter- 
action; the Casimir thermal force (the excess of Casimir energy at a given temperature 
compared with the Casimir energy in the zero temperature case) between two isolating 
bodies is not necessarily monotonic in their separation, as seen in the sphere-plate 
and cylinder-plate cases [151]. In addition, for some geometries, intervals of negative 
entropy appear, as in the case of two parallel plates described by the Drude model 
[152] or, as recently shown, in the interaction between a perfect metal plate and sphere 
[153]. 

In this Chapter we study the Casimir effect between two spheres in the large sep- 
aration approximation using different models for the electric permeability: (1) the 
perfect metal model, (2) the plasma model, and (3) the Drude model. As a result, 
we find negative entropies in certain ranges of temperature and separation between 
the spheres for the perfect metal model, and for low penetration length for the plasma 
model. In addition, we find nonmonotonic behavior of the entropy with the separation 
while the force is attractive for all separations, making it appear as though the natural 
evolution of the system tends to increase the entropy in certain ranges of temperature 
and separation. For long plasma length and for the Drude model we do not find this 
anomalous behavior of the entropy. We also discuss the thermodynamical meaning and 
consequences of negative entropies in Casimir effect. 

The remainder of the Chapter is arranged as follows: In Sect. 8.1, we describe the 
multiscattering model used herein to obtain the Casimir energies and entropies for the 
two spheres. In Sect. 8.2, we obtain the large separation limit of the Casimir energy 
between perfect metal spheres, as well as the entropy and force. We also obtain en- 
tropies at smaller separations numerically. We study the plasma model system in Sect. 
8.3 and the Drude model in Sect. 8.4. We discuss the thermodynamic consequences of 
these results in Sect. 8.5. Finally, we discuss the results obtained in the Conclusions. 
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The contents of this Chapter is based on the work pubhshed in [154]. 



8.1 Multiscattering formalism of the Casimir en- 
ergy 

To calculate the Casimir energy, entropy, and forces between two spheres, we employ 
the multiscattering formalism for the electromagnetic field [21] [22]. This formalism 
relates the Casimir interaction between objects with the scattering of the field from 
each object. The Casimir contribution to the Helmholtz free energy at any temperature 
T is given by 

oo ' 

E = kBTY,^og\l-NM\, (8.1) 

n=0 

where k„ = ^ are the Matsubara frequencies and Ay = 2i!kBT thermal wave- 

length. The prime indicates that the zero Matsubara frequency contribution has height 
of 1/2. All the information regarding the system is described by the N matrix. For a 
system of two objects, this matrix is N = T1U12T2U21. Tj is the T scattering matrix of 
the i*^ object, which accounts for all the geometrical information and electromagnetic 
properties of the object. Ujj is the translation matrix of electromagnetic waves from 
object i to object j, which accounts for all information regarding the relative positions 
between the objects of the system. 

For a sphere of radius R with electric and magnetic permeabilities e and /x, the T 
matrix is diagonal in {imP, i'm'P') space, with elements given by 
^im/'m' = -^u'^mm'^pp'Tf^, whcrc the are defined as [22] 

keiKR)dR{Rte{nKR))-ndR{Rke{KR))ieinKRy ^ ' ^ 



rj,E _ k{i^R)dR{Riii{nnR))-edR{Riii{nR))ii{nnR) 

ke{KR)dR{Rie{nKR))-edR{Rh{KR))ie{nKRy ^ ' ^ 



where ii{x) = y^/^.,.! (x), ke{x) = \ ^f^K ^j^iXx) ^ and n = ^/Jle. Expressions for the 
matrices can be found in [22] and [124]. 
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8.2 Perfect metal spheres 

In the perfect metal limit, we apply e — )■ oo for any fx to Eqs. (8.2) and (8.3). In this 
case, the T matrix elements are independent of /i, taking the well-known universal form 

- -'3«'<3_'2i^,^.(«:i?)' ^^-^^ 



8.2.1 Casimir energy in the large separation limit 

To obtain the large separation limit of the Casimir energy, we need the dominant part 
of the T matrix in this limit. We define the adimensional frequency q hj n = q/d. The 
main contribution in the large separation limit comes from the lowest-order expansion 
term of the T matrix elements in 1/d. As it is known that, at small k, the T elements 
scale as k^^"*"^ [21], the dominant contribution comes from the dipolar polarizabilities 
part of the T matrix, taking the form 

3 \ d ' 3 V d 



1 qRY ^EE 2/?i?, 



By the use of the universal relationship between determinants and traces, log|y4| = 
Trlog(yl), and applying a Taylor expansion in terms of ^ of Eq. (8.1), we obtain 

oo ' 

E = -kBTj2 Tr (N(A5;^n)) . (8.7) 

n=0 

Using the translation matrices in a spherical vector multipole basis [22] and the large 
separation approximation of the T matrix, the trace of the N matrix is obtained by 
straightforward calculus. Here we denote with a sub-index T the results valid for all 
temperatures, with a sub-index the results in the quantum limit (T — )■ 0), and with 
a sub-index cl the results in the classical limit {h — t- 0), which is equivalent to the 
high T limit. Carrying out the sum over Matsubara frequency, we obtain the Casimir 
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contribution to the Helmholtz free energy as 



hcR^ ze- 



362 



27rrf7 2 (e2- - 1)^ 
(2 (15 - 29^2 + 99^^) cosh(2) + 15 cosh(5z) 

+ (-45 + 58^^ + 18z^) cosh(3z) (8.8) 

+2Az (6^2 - 5 + (5 + 3z^) cosh(2;2)) sinh(z)) , 

where z = d/Xx- We also define the adimensional Casimir energy as Ead{z) = ^^Et- 
From this result, the quantum (T — t- 0) and classical {h 0) limits with their first 
corrections are easily obtained as 

UShcR^ SksTTi'RWkBTY 
^' = -^6^^ Tld—[-h^) ' ^^-^^ 

Ed = -^^^^ - ^bT^ (15 + 30z + 29z^ + ISz"^ + 9z^) e'^^ (8.10) 

where z = d/Xx- Note that in Eq. (8.9), the first correction to zero temperature case 
is proportional to T^, contrary to the plate-sphere case [153] [151] and to the cylinder- 
plate case [151], where a result proportional to T'^ were obtained for both systems. The 
result presented here is new for perfect metal spheres, but expected, because it has been 
obtained also for compact objects in the diluted limit [28] [11 5]. As corrections have 
the same negative sign of the Casimir energy in both limits, they describe an increase 
of the magnitude of the Casimir energy in high and low temperature limits. It is no 
longer the case if we study the energy for all temperatures. In fact, for some ranges 
of separation and temperature, thermal photons tend to reduce the Casimir energy 
between the spheres, as shown in Fig. 8.1. This effect is captured in the contribution 
of the next term to the quantum limit, this term is AEq = +228M^^^ (¥)^ 
tends to reduce the magnitude of the energy. Fig. 8.1 is an indicator of the appearance 
of negative entropy in this system because of the negative slope of the energy curve. 
It is clear in Fig. 8.1 that for some temperatures and distances, the Casimir energy 
between spheres is lower than in the zero temperature case. But it is less evident the 
validity Eq. (8.9). In fact, there is a tiny increase of E/Eq at low ^ until reaching 
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Figure 8.1: Casimir energy between perfect metal spheres as a function of ^ compared 
with the energy at zero temperature. The dotted curve is the quantum limit, the dashed 
curve is the classical limit, and the solid curve is the asymptotic finite-temperature Casimir 
energy. Note that these curves are independent of the radius of the spheres but are only valid 
in the large separation limit. 



a local maximum at ^ !^ 1.0388 oi E/Eq ^ 1 + W It is not visible in Fig. 8.1 
because the difference of scales. The local minimum is reached at ~ 3.21733, where 
E/Eq 1 — 3.46 X 10~^. This minimum is reached, in the important cases of room 
temperature (T = 300K) and at the boiling point of N2 (T = 77K) at dT=3oo = 3.89/im 
and dT=77 = 15.21yum respectively. 



8.2.2 Casimir entropy in the large separation limit 

In the canonical ensemble, the entropy is defined as S* = —drE, where with an abuse of 
the notation we denote the Helmholtz free energy by E. In the large separation limit, 
the Helmholtz free energy depends on the adimensional variable z = = 27r^^^, so 
we can write the entropy as 

dzdE ^dksdE 
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Figure 8.2: Casimir entropy between perfect metal spheres as a function of compared 
with the classical limit. The dashed curve is the classical limit, and the solid curve is the 
asymptotic finite-temperature Casimir entropy. Note the regions of negative entropy and 
negative slope of the entropy with the parameter ^ 



and define the adimensional entropy as Sad{,z) = -^^^S = —dzEad{,z). From this result, 
the quantum (T — )■ 0) and classical {h — )■ 0) limits are easily obtained as 



5 



1 T?^ T?^ 

S,i = + (15 + 30;^ + 27z^ + Uz^ + 9z^ - 18z'>) e'^^ (8.13) 

where z = j^. So, the entropy is a growing function with temperature in both limits, 
but this is not the case for all temperatures (the next low temperature expansion term 
of Eq. (8.12) is ASq = - ^^^^^f kedR^ ^^^^^ already indicates that the entropy 

could change its growing behavior with the temperature), as we can observe in Fig. 
8.2, where a region of negative entropy and another region of negative slope of the 
entropy are observed. 

Because of the limit at low temperature of the entropy, we know that the entropy 
is positive for low T (not seen in Fig. 8.2 because it is small compared with S/Sd, but 
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it can be observed in Fig. 8.3), so there are three points where S = 0, including the 
origin. These two new points where the entropy nulls correspond to the local maximum 
and minimum observed in the Casimir energy in Fig. 8.1. 

Negative entropy of the Casimir effect has already been obtained between Drude 
parallel plates in [152] and in [155], and more recently between a perfect metal plate 
and sphere in [153] and between a Drude sphere and plate in [156] and [157]. This is 
the first time, to the best of our knowledge, that negative entropy appears between 
spheres because of the Casimir effect. 

These results are only valid when the separation between the spheres is large com- 
pared with their radius, regardless of the radius of each one. Therefore, it is possible 
that the interval of negative entropy would disappear when the separation between the 
spheres reduces, regardless of the temperature of the system. A numerical exploration 
of entropies at smaller separation has been performed to verify this issue. 

8.2.3 Numerical study at smaller separations 

As noted in the previous Subsection, asymptotic results are no longer valid when the 
separation between the spheres becomes comparable to their radius. For this reason, 
a numerical study of entropy was performed for these cases. We computed Eq. (8.1) 
numerically for all temperatures from T = until reaching the classical limit for fixed 
ratio between the radius R and separation, r = ^. 

For spheres, the T matrices are diagonal but infinite matrices, so a cutoff in {i, [') 
space is needed to obtain finite-dimensional matrices. In addition, another cutoff in 
Matsubara frequency is needed to obtain a finite series. 

The proposed method is an asymptotic approximation, at small separations more 
and more modes are needed to obtain convergent results. This means that there exists 
a minimum separation between the spheres below which we are not able to take into 
account enough multipoles to obtain good results. We use multipoles up to £ < 15, 
which means that we are restricted to r^ax = 7r~ — 0.45, when contact is reached at 
'^contact = ^ ^ = ^ = 0-5, where the energy diverges. 

In the small separation limit, the proximity force approximation (PFA) is known to 



8.2 Perfect metal spheres 



153 




Figure 8.3: Log-log plot of the absolute value of the entropy divided by its classical limit 
for perfect metal spheres at constant r as a function of Starting from the asymptotic 
solid curve (r — )• 0), we increase r in steps of 0.1 up to r = 0.4 (dotted curve). Dashed curve 
is the case r = 0.41, where the interval of negative entropy has disappeared. The curve for 
r = 0.45 is also shown 



be a good approximation to the Casimir energy [21]. It is also known that perfect metal 
plates do not experience negative entropies, so we do not expect to observe negative 
entropies between spheres in the small separation limit. 

In Fig. 8.3, the entropy of the system of two perfect metal spheres is plotted as 
a function oi z = for constant r = ^. The large separation result is shown too. 
We choose a log-log representation of the absolute value of the entropy divided by 
its corresponding classical limit. Therefore, zeros are observed as log divergences, and 
we can also observe the cases of negative entropy. Starting in the large separation 
regime, we observe an interval of negative entropy. As we increase r (reducing the 
separation between the spheres), the region of negative entropy tends to reduce until 
it disappears between r = 0.40 and r = 0.41. Power-law decay of the entropy (8.12) 
at low temperatures is observed as a linear decay of the curve at low and constant 
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behavior in the high-temperature hmit (8.13) is also reached in the simulation. 

8.2.4 Casimir force in the large separation limit 

In the previous Subsection we demonstrated that intervals of negative entropy appear 
due to the Casimir effect between perfect metal spheres. In these cases, for any given 
temperature, a minimum of entropy exists for a given, nonzero separation. Naively, 
this would imply a violation of the third law of Thermodynamics if the Casimir force 
is not zero at the minimum of the entropy. However, this is not the will be 

explained in Sect. 8.5. In addition to that, we find that the force is always attractive, 
independent of the increase or decrease of the entropy, which also would naively imply 
a violation of the second law, because the system can be enforced to perform a process 
in which the entropy is reduced instead increased. However, this is not the case, as 
will be also explained in Sect. 8.5. However, the appearance of negative entropy does 
have an effect on the force. The asymptotic Casimir force F = —ddE can be written 
in terms of the adimensional Casimir energy as 

Fad{z) = ——F = lEad{z) - zd,Ead{z), (8.14) 

where z = j^. In Fig. 8.4, the adimensional asymptotic force between the perfect 
metal spheres compared with the zero-temperature force is plotted as a function of j^. 
Here, for constant temperature, we observe nonmonotonic behavior of the force with 
the adimensional parameter 

Nonmonotonic behavior of the adimensional Casimir force implies nonmonotonic 
force behavior with temperature, but not with separation, because of the extra depen- 
dence of the force on the separation in Eq. (8.14). In fact, it is easy to verify that the 
force behaves monotonically with separation, and the nonmonotonicity of the entropy 
with separation implies nonmonotonic behavior of the force with temperature, because 

dF ^ ^dS 

dT dTdd dd' ^ ' 

so the appearance of negative slopes of the entropy with separation implies nonmono- 
tonicity of the Casimir force with temperature, despite the attractive force for all 
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Figure 8.4: Asymptotic force divided by its zero-temperature limit as a function of j^. 
The dotted curve is the zero-temperature result. The dashed curve is the classical limit, and 
the solid curve is the result at finite temperatures. The nonmonotonic behavior of the force 
with temperature results from the negative slope of the solid curve at constant separation 



separations and temperatures. 

As observed in Fig. 8.4, the force between the spheres is always attractive, but it 
is not monotonic with temperature; asymptotically, for any given temperature, there 
exists a range of separations for which the force decreases with temperature. This is 
the first time, to the best of our knowledge, that nonmonotonic behavior of the Casimir 
force with temperature has been described between compact objects. Nonmonotonic- 
ities of the force between a plate and cylinder and between a plate and sphere were 
already obtained in [151], but in that case the nonmonotonicity already appears for the 
scalar field. Nonmonotonicity does not appear between spheres for the scalar field; this 
is a characteristic effect of the electromagnetic field, because cross-polarization terms 
of the Casimir energy are essential for this nonmonotonicity to appear. 
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8.3 Plasma model 

In this Section, we assume that the electric susceptibihty of both spheres is described 
by the plasma model, i.e., 

e(zcK) = l + \y^ (8.16) 

and /i = 1. To obtain the large separation limit of the Casimir energy, we need the 
dominant part of the T matrix in this limit. The main contribution comes from the 
dipolar polarizabilities part of T matrix, taking the form [153] [156] 

TK™' = (3 + 2/^ - 3?/coth {y)) , (8.17) 

Tff,^, = -('^V, (8.18) 



-lm,lm' 3 1^ ^ 

with y = 27r^. The coefficients of the magnetic polarizability depend on the plasma 
frequency, tending to the perfect metal result as Ap — ?■ and to zero in the transparent 
limit Ap — i- oo [156]. 

8.3.1 Casimir energy in the large separation limit 

Once we have the asymptotic T matrix, the Casimir energy can be obtained by a 
straightforward but long calculation. It is possible to obtain analytical results for 
finite temperatures, but they are too long to show here. The zero-temperature Casimir 
energy is 

hcR'^ 

(207 + 222y2 + U3y^ - AlAy coth(2/) + 207?/^ cot\i{yf - 222y^ coth(?/)) , 

and the high-temperature limit is given by 

Ea = -kBT— M +— (3+i/2_3ycoth(?/)) J , (8.20) 

with y = 2txy^. The perfect metal limit is reached for Ap — )■ 0, as expected. When the 

plasma wavelength Xp <^ R, only the electric sector of the T matrix contributes, so in 

this case the energy in the zero- and high-temperature limits is given by 

23hcR^ ^ , ^3R^ , , 
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8.3.2 Casimir entropy in the large separation limit 

The entropy of the system is obtained in the same way as for the perfect metal case, 
but for the plasma model we have two different regimes. When Xp <^ R, we reach 
the perfect metal limit, obtaining negative entropy and nonmonotonic behavior of the 
entropy with separation and temperature as before. However, when Xp ^ R, the 
spheres are transparent to the magnetic field and the problem reduces to a scalar 
problem. In this case, the entropy is always positive. In Fig. 8.5, the adimensional 
asymptotic entropy, compared with its classical limit, is plotted as a function of ^ for 
these two limits. In Fig. 8.6, we present the points (^, ^) corresponding to zeros 
of the entropy (continuous curve), its temperature derivative (dashed curve), or its 
separation derivative (dotted curve). As discussed in the case of perfect metal plates, 
the anomalous behavior of the entropy found here would naively imply a violation 
of second and third laws of Thermodynamics, because the Casimir force is always 
attractive, irrespective of the slope of the entropy. However, this is not the case, as 
we will discuss in Sect. 8.5. As noted above, in the perfect metal limit we obtain a 
region of negative entropy, while for low Xp this region disappears. The point where 
negative entropies appear depends on the magnetic susceptibility fi. For /i = 1, negative 
entropies appear around Ap ~ 2R. 



8.4 Drude model 

In this Section we assume that the electric susceptibility of both spheres is described 
by the Drude model, i.e., 

<^c«:) = 1 + (8.22) 

and yU = 1. To obtain the large separation limit of the Casimir energy, we need the 
dominant part of the T matrix in this limit. The main contribution comes from the 
dipolar polarizabilities part of the T matrix, taking the form [153] [156] 



^ ( ^ 1 , Ti„,w-3 ( ^ ) • (8.23) 
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Figure 8.5: Entropy between plasma spheres as a function of ^ divided by the classical 
limit for the two limit cases. The solid curve is the entropy in the perfect metal limit A^j <C i?, 
where a region of negative entropy and nonmonotonicities of entropy with separation and 
temperature are present; the dashed curve is the entropy in the transparent limit \ji ^ it!, 
where the usual monotonic behavior of entropy is shown 



Now, the response of the material changes dramatically from the case of the plasma 
model. In principle, one would expect to obtain the T matrix for the plasma model 
as cr — > oo, but it is obvious that this is not the case. Not only do we not recover 
the plasma model in this limit, but also the T matrix diverges. The reason is that 
the plasma model is some kind of singular limit of Drude model, with a qualitative 
different behavior. Now the dominant contribution is given by the electric sector of the 
T matrix, whereas the magnetic part does not contribute for asymptotic energies. 



8.4.1 Casimir energy and entropy in the large separation limit 

Once we have the asymptotic T matrix, obtaining the Casimir energy is a straight- 
forward calculation. Carrying out the sum over Matsubara frequency, we obtain the 
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the temperature derivative of the entropy (dashed curve), and the separation derivative of 
the entropy (dotted curve) 



Casimir contribution to the Helmholtz free energy as 

Et = F X 

(e2- - 1)^ 

((6 - lOz^ + 22z^) cosh(2) + (-36^ + 12z^) smh{z) 

+ {I2z + 4^3) sinh(3z) + 3 cosh(5z) + (-9 + lOz^ + 2z^) cosh(3z)) , (8.24) 

where z = d/Xx- The zero- and high-temperature Casimir energy hmits are given by 

23hcR^ ^ , ^3i?6 , , 

E„ = -^. E,^-k,T^. (8.25) 

These asymptotic Casimir energies are universal because they do not depend on any 

material property. The lack of magnetic polarizability contributions is characteristic 

of the Drude model, leading to an effective scalar Dirichlet problem. The entropy is 

obtained as in the perfect metal case. As only the electric polarization contributes, we 

expect the usual behavior of the entropy as a monotonic function of as shown in 

Fig. 8.7. The zero- and high-temperature Casimir entropy limits are given by 

^ rOAkBRWrrkBTV imGkBdRW nksTy 
^°-^15^l^j 157^ l^j ' ^^-^^^ 
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Figure 8.7: Entropy as a function of the adimensional parameter ^ divided by its 
classical limit for the large separation limit of the system of Drude spheres. Entropy is 
positive and monotonic in temperature and separation 



r.6 

Sci = 3kB^. (8.27) 

In this case, the correction term to the quantum Umit tends to reduce the entropy, but 
it is not strong enough to change the sign of it, as seen in Fig. 8.7. 



8.5 Thermodynamical consequences 

In this Section we discuss the thermodynamical consequences of the obtained results. 
In this article we have obtained and compared the large separation limit of the Casimir 
energy and entropy for two spheres using three different electric susceptibility models. 
For perfect metal spheres, at any nonzero fixed temperature, we observe an interval 
of separations for which entropy is negative, while at zero temperature the entropy is 
always zero. 

In addition, the Casimir force is attractive for all separations and temperatures. 
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So, we could naively think that we have possible violations of the third and second 
laws of Thermodynamics, due to the existence of processes where the entropy of the 
system tends to increase and to the negative entropy intervals at finite temperature 
and distances, respectively. 

According to the Krein formula [25], we know that the Helmholtz free energy of 
the electromagnetic field has three independent additive contributions (see Eq. (8.36) 
of the Appendix): one is from the thermal bath, being proportional to the volume of 
the space [107]. Another is the sum of contributions of objects immersed in the bath 
considered as isolated objects, each contribution being also a function of the volume 
and surface of each object [107]. The third is of geometrical nature, which we could 
call the Casimir part of the Helmholtz free energy and that we actually calculate in 
Eq. (8.1). This term depends on the electromagnetic nature of the objects in the 
system, but it has a geometric nature, because it depends on the relative separations 
and orientations between the objects, being zero iff the relative separations between 
all the considered objects become infinite. 

Considering the whole system, the nonmonotonicity of the Casimir entropy with 
temperature is compensated by the contribution of the vacuum, because one scales 
with the global volume [107] (Eq. (8.37)) and the other with the separation between 
the spheres (Eq. (8.13)). This is not the case for the nonmonotonicity with separation. 
At constant temperature, the entropies of the thermal bath and of each object remain 
constant, but the Casimir force is always attractive while the entropy can increase 
or decrease with separation. Therefore, the potential violation of the second law of 
Thermodynamics still requires an explanation. 

The Krein formula also states that internal sources of entropy inside the spheres 
cannot appear in order to compensate these regions of anomalous behavior of the 
entropy, unless these internal sources are independent of the separation between the 
spheres. 

The second law of Thermodynamics states that global entropy must increase for any 
process, but only in closed systems. As we are working in the canonical ensemble, we 
are implicitly assuming that there exists an external reservoir which keeps our system 
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at constant temperature, so the system is not isolated and entropy can increase or 
decrease without violation of the second law. In the canonical ensemble, the condition 
equivalent to the second law is that the global Helmholtz free energy must decrease 
for any process, and this is true for the studied system. Therefore, the appearance of 
nonmonotonic entropy behavior in the canonical ensemble just implies nonmonotonic 
behavior of the force with temperature, as seen in Eq. (8.15). 

8.6 Conclusions 

In this Chapter, we have obtained and compared the large separation limit of the 
Casimir energy and entropy for two spheres using three different electric susceptibility 
models. For perfect metal spheres, at any nonzero fixed temperature, we observe an 
interval of separations for which entropy is negative, while at zero temperature the 
entropy is always zero. 

We can trace the origin of this anomalous effect to the functional form of the Casimir 
energy in this large separation regime. It consists of four components, two of which are 
of attractive nature, equivalent to the sum of two scalar problems, one for the electric 
polarization and the other for the magnetic polarization. The other two components 
come from the cross-coupling between the electric and magnetic polarizations of the 
two spheres, and tend to reduce the Casimir energy between the spheres (Eq. 12 
of [21] and Eq. 6 of [28]). These cross-coupling polarization terms are responsible 
for the impossibility of factorization of the electromagnetic Casimir energy into two 
equivalent scalar problems in general, and appear because the translation matrix Uq/? 
is not diagonal in polarization space. In this article we show that these cross terms 
not only reduce the Casimir energy, but also for some separations and temperatures, 
their contribution to the entropy is greater than the contributions of direct coupling 
between the electric and magnetic polarizabilities, resulting in an interval of negative 
entropy and nonmonotonic behavior of the entropy with separation and temperature. 

In addition, the Casimir force is attractive for all separations and temperatures. 
So, we could naively think that we have possible violations of the third and second 
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laws of Thermodynamics, due to the existence of processes where the entropy of the 
system tends to increase and to the negative entropy intervals at finite temperature 
and distances, respectively. 

In Sect. 8.5, having into account the complete thermodynamical system, that 
the system is described by the canonical ensemble, and with the help of the Krein 
formula and Weyl formula, we have demonstrated that there are not violations of 
second and third laws respectively. Therefore, the appearance of nonmonotonic entropy 
behavior in the canonical ensemble just implies nonmonotonic behavior of the force with 
temperature, as seen in Eq. (8.15). 

The interval of negative entropy for the spheres appears because of the cross- 
coupling between the polarizations of the electromagnetic field, which leads us to con- 
clude that it is a characteristic of the electromagnetic field and does not have an analog 
in the Casimir effect due to scalar fields. 

Applying the PFA to this problem, we obtain that, for perfect metal spheres near 
contact, the entropy is always positive, so we performed a numerical study of the 
entropy between perfect metal spheres of equal radius at smaller separations. The 
results showed that there exists a minimum separation between the spheres for which 
the negative entropies disappear. In addition, the region of negative slope of the entropy 
with separation disappears for another smaller separation. 

We have also obtained the energy and entropy for the plasma model, for which 
similar results are obtained for plasma penetration length Ap < 2i?, while for Xp > 
2R, the negative entropies disappear, and for another greater Xp, the nonmonotonic 
behavior of the entropy with separation also disappears. 

In the transparent limit [Xp ^ R), the magnetic polarization does not contribute 
to the asymptotic Casimir effect, so the problem reduces to a scalar field problem and 
we find that entropy is a positive monotonic function. 

When we study the asymptotic Casimir effect with Drude model spheres, the sys- 
tem is qualitatively different. The magnetic polarization does not contribute to the 
asymptotic Casimir effect because it depends on rf""^ instead of the dependence 
of the electric polarization term. Then also for Drude spheres the asymptotic limit 
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reduces to a scalar field problem, resulting in a positive monotonic entropy. 

Nonmonotonicities of the Casimir force between objects are not unusual. As the 
Casimir effect is not pairwise additive, the interaction between two objects is affected 
by the presence of a third, leading to nonmonotonicities of the Casimir force between 
cylinders [158] [159], or between spheres [130] in the presence of a plate. The non- 
monotonicity presented in this article has a different nature, as it does not come from 
interactions with a third object but rather from correlations with temperature, as seen 
in Eq. (8.15). 

In addition, in the assumptions of the multiscattering formalism it is implicitly 
assumed that the objects are stationary, at fixed positions in space, so all results 
shown here are valid for quasistatic processes. Therefore, we must carefully consider 
when this quasistatic assumption does not apply, because the stationary system could 
abandon equilibrium, requiring more careful study [55] [160]. 



8.7 Appendix A: Krein formula from multiscatter- 
ing formalism 

In this Appendix, we will derive the Krein formula of the electromagnetic field [25] 
from the multiscattering formula [21]. 

Krein formula factorizes the density of states of the electromagnetic field in three 
parts [25]. Here we present a derivation of such formula from the Multiscattering 
formalism of the Casimir effect [21] [22]. 

The partition function {Z = X]^o'^") '^^ in the presence of gen- 

eral dielectrics of arbitrary geometry is obtained, before regularization of the Casimir 
energy, as [22] 

where |S'„|~^ = | A + ^ is the n— Matsubara frequency contribution to the partition 
function of the thermal bath, and M. is the next N x N block matrix 

-Ma/3 = So^f^J-' + {5ap - (8.29) 
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If we multiply and divide Eq. (8.28) by |A^oo| (being Aioo the same Ai matrix as 
above, but with each object at an infinite distance from each other), the Helmholtz 
free energy of the whole system can be written as (/^J-" = log{Z)) 



n=0 



l0g\Sn\ +log|A^oo| +log 



\M\ 

l-Moo 



(8.30) 



where 



log \ Sn\ is the contribution of the thermal bath in absence of immersed objects, 
which leads to the Planck formula of the blackbody spectrum. 

log|A^oo| = X]a=i 1*^S I'^al is contribution to the Helmholtz free energy of 
each dielectric object immersed in the bath considered as an isolated object. 



log 



\M\ 
\M^\ 



log |1 — N| is the Casimir part of the Helmholtz free energy, which 



depends on the electromagnetic nature of the objects in the system, but it has a 
geometric nature, because it depends on the relative separations and orientations 
between the objects, being zero iff the relative separations between all the con- 
sidered objects become infinite. This term is the Casimir energy that we actually 
calculate in Eq. (8.1). 

The energy of a quantum system is related with the Density of States p{co) by [25] 



E 



duj-ujpiuj). 



(8.31) 



Using Eq. (8.30), the same energy in the zero temperature case can be written as 

he 



E 



27r 



log 1^,1 log |T„|+ log |1-N| 



(8.32) 



Because the dispersion relation of the massless EM field u = ck, and integrating by 
parts, we obtain 



E= dio-co-x 
Jo 



2 n' 



(8.33) 



X 



-d^ log + ^ log |T,| - log |1 - N| 



a=l 
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Then, because Eqs. (8.31) and (8.33), we can relate the Density of States with scat- 
tering properties of the system as 



p{uj) = - 

TT 



-d^ \og\SJ + J2 du. log |T,| - 9^ log |1 - N| 

a=l 

by the use of the universal relationship between trace and determinant log|y4| 
Trlog(74), this result can be written as 



8.34) 



p{u) = - 

TT 



-Tr (Sj^o'^S^) +5^Tr(T-i6'^T,) +Tr((l-N)-'a^N) . (8.35) 

As a conclusion, we can factorize the Density of States as the sum of three independent 
terms 

n 

P(w) = Po(w) + J^Pa(w) +Pc(w), (8.36) 

Q=l 

where 

• po(i^) = —-^dujlog \Suj\ is the Density of States of the free EM field. 

• Pa{oj) = ^c^oj log is the change of the vacuum Density of States of the EM field 
because each a— dielectric object considered as a isolated object in the medium. 

• pc{,<^) = ^~<9i^ log |1 — N| is the change of the Density of States related with the 
relative positions of the dielectric objects, which is zero when all the objects 
are infinitely apart from each other. pc{'^) is the source of the Casimir effect. 

Eq. (8.36) is the Krein formula [25], derived from the Multiscattering formalism of the 
Casimir effect for the EM field. In [25], Eq. (8.1) at zero temperature is obtained from 
the Krein formula removing the divergent contributions. Here we have followed the 
opposite way, a derivation of the density of states of the a thermal bath with intrusions 
and it factorization in Eq. (8.36). 

Terms po{uj) and Pai^^) diverges with u because main contribution to Weyl formula 
for 3D volumes of the EM field [107] [44], 

pvico) = -^yeico), (8.37) 



8.7 Appendix A: Krein formula from multiscattering formalism 
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but pc{<^) converges if our system consists on N compact objects, because in this case, 
N can be demonstrated to be a trace class operator for all w, then the determinant is 
well defined [23]. 
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"/ think that it is a relatively good approximation to truth - which is 
much too complicated to allow anything but approximations." 

John von Neumann 



9 

Casimir energy between non parallel 

cylinders 



Casimir energies between cylinders have been studied from long time ago [161], but 
always between parallel cylinders by the help of the axial symmetry of the system [128] 
[22]. However, the case of non-parallel cylinders has not been studied until recently 
in [162] and [163], where they studied asymptotic results. In addition to that, in [37] 
the authors performed a numerical method and studied the case of two perpendicular 
elongated objects of length L, when they incremented the length L of the objects, fixed 
at a given distance, the energy finally saturates to a finite result independent of the 
length. 

In this Chapter, we present the formalism to study the Casimir energy in systems 
with non parallel cylinders (see Fig. 9.1) with the multiscattering formalism [21] [22]. 
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This formalism will let us perform the numerical evaluation of the energy for all dis- 
tances, not only asymptotic results. In particular, we study two non parallel cylinders 
for the scalar case for cylinders subject to Dirichlet or Neumann boundary conditions 
and the electromagnetic case for perfect metals. This study can be easily extended to 
dielectric cylinders. 

Our goal is the calculation of Casimir energy between two non parallel cylinders, 
extending the study of Casimir energies to new geometries. Parallel cylinders case is 
easily reobtained as a limiting case of our formula. 

The interest of the proposed problem is to expand the study of Casimir energies for 
systems where non parallel cylinders appear. Here we remind that in the microscopic 
world, a fiber could be actually considered as an infinite cylinder. For example, Casimir 
energy could have a contribution in the stability orientation of cylinders, as in the case 
of carbon nanotubes. 

We follow this plan for the Chapter: In Sect. 9.1 we present the system under study. 
In Sect. 9.2, in order to use the multiscattering formalism, we obtain the translation 
matrices between cylindrical multipole basis when their axial axis are tilted an angle 
7 and their origin are separated by a distance d, for scalar and electromagnetic cases. 
In Sect. 9.3 we obtain the asymptotic far distance approximation of electromagnetic 
and scalar Casimir energies at zero and high temperature limits. In Sect. 9.4 we 
present some numerical evaluation of Casimir energies between perpendicular and half- 
perpendicular perfect metals cylinders at zero temperature. The mathematical details 
of the derivation of translation matrices are left to Appendix 9.6. In Appendix 9.7 
we show how to cover the limit of parallel cylinders from our formalism. Finally, in 
Appendix 9.8 we obtain the Proximity Force Approximation of the studied system, at 
the zero and high temperature limits. 

The content of this Chapter is based on the work published in [164]. 



9.1 Geometry of the system under study 
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9.1 Geometry of the system under study 

Recently, a multiscattering formalism of Casimir energies has been developed [21] [22]. 
It relates the Casimir interaction between objects with the scattering of the field from 
each object. The Casimir energy for a two objects system at any temperature is 

oo ' 

ET = kBTY^ log|l-NK)l, (9.1) 

n=0 

with Kn = 2iT^^n as the Matsubara frequencies. The tilde means that the zero 
contribution has a height of 1/2. The system is described by the N matrix. For a 
system of two objects, this matrix is N = T1U12T2U21. Tj is the T scattering matrix 
of i— object, which accounts for all the geometrical information and electromagnetic 
properties of each object. Vij are the translation matrices of electromagnetic waves 
from object i to object j, which account with all the information of the relative positions 
between the objects of our system. 

From Eq. (9.1), the quantum (T — > 0) and classical {h — )■ 0) limits are easily 
obtained as 

he f°° 

Eo = — dk\og\l-N{K)\ (9.2) 

^TT Jo 

and 

Ee/ = ^log|l-N(0)| (9.3) 

respectively. Here we denote with a sub-index T the results valid for all temperatures, 
with a sub-index the results in the quantum limit (T — )■ 0), and with a sub-index cl 
the results in the classical limit (/i — > 0), which is equivalent to the high T limit. 

In the case of cylinders, T matrices are generally known [22], as well as U matrices 
for parallel cylinders [22]. For this particular case, Casimir energy scales with the 
infinite length L of the cylinders. It is a natural question if this scaling with the length 
remains for non parallel cylinders. On one hand, for the non parallel case, the scale of 
cylinders is still the length L, but on the other hand, for the parallel case, all points 
of the cylinders were equally relevant in order to contribute to the Casimir energy and 
it is no longer true for the non parallel case, where necessarily there is a principal 
contribution of the nearest points and there are points infinitely far away. 
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The Proximity Force Approximation (PFA) applied to the electromagnetic (EM) 
Casimir energy leads to the next results for perfect metal parallel cylinders in the zero 
[129] and high temperature limits: 



EPP^ = -^hcL.I ^ (9.4) 
1920 V{d~2Rf' ^ ' 



where R is the radius of the cylinders and d is the distance between their centers. 
These energies scales with the length of the cylinders L. If we calculate the PFA for 
the non parallel case, we obtain (see Appendix 9.8) 

720|sin(7)| (c/-2i?)2' "^""-^^ 



C(3) kBT R 
''^^ 4 |sin(7)| (rf-2i?)' ^ ' 

For the non-parallel case, the energy scales with 1/| sin(7)| instead with L. This result 
suggests that the parallel and the non-parallel problems have different scales. How- 
ever, PFA is an approximation valid for short distances between the bodies. At large 
distances, the energy can be qualitatively different, as it is the case of prolate spheroids 
[165]. 

For non-parallel cylinders, Eq. (9.1) remains valid, as well as T matrices. But now 
U matrices depend not only on the relative distance d between cylinders, but also on 
their relative orientation. Then, in addition to the translation contribution, we have to 
consider the rotation around the translation axis an angle 6. For non crossed cylinders 
this can be demonstrated to be the more general transformation. 

In particular, we will study the case of a translation over a d vector of the (x, y)- 
plane, which will also be the rotation axis (see Fig. 9.1), then the change of spatial 
coordinates is 

x' = M^,-,x + d, (9.8) 



9.2 Translation matrices between non- parallel cylindrical multipole 
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with 



dcos{ 
dsin( 



(9.9) 



v 



as the translation vector over the (x, ?/)-plane, and 



■d,7 
/ 



cos(7) + cos^(6'jj)(l - cos(7)) 


sin(%) cos(%)(l - cos(7)) 


sin(6'ji) sin(7) 


sin(6'jj) cos(^^)(l — cos(7)) 


cos(7) + sin^(6'jj)(l — cos(7)) 


- cos(0ji) sin(7) 


- sin(6'jj) sin(7) 


- cos(6'ji) sin(7) 


cos(7) 



(9.10) 
\ 



V 



as the rotation matrix over the d axis an angle 7. We also will use the transformation 
in the momentum space 

( y \ ( k \ 



d,7 



(9.11) 



The derivation of the translation matrices will be given in the next Section. 



9.2 Translation matrices between non-parallel cylin- 
drical multipole basis 

In this Section, we derive the transformation matrices U between scalar and vectorial 
cylindrical multipoles of non parallel basis. In particular, we need to transform the 
outgoing cylindrical wave functions with imaginary frequency of a B' basis in terms of 
regular cylindrical wave functions with the same imaginary frequency of another basis 
B translated by a distance d and whose axes are tilted against each other by an angle 
7 , as shown in Eq. (9.8). 

We perform this analysis for scalar and vectorial cylindrical waves. 
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Figure 9.1: System of non-parallel cylinders studied in this Chapter for the particular 
case 9ji = |. The coordinate basis are also represented. 

9.2.1 Scalar waves 

We consider regular and outgoing cylindrical wave functions, which on the imaginary 
frequency axis are given by 

€iM) = In{pp)e''''e^'^\ (9.12) 

0;;XW = ^"Me^"'e^'=^^ (9.13) 

with p = \/k1 + ky and p = a/s^ + y'^. The outgoing waves at position x' can be 
expanded in terms of regular waves at x by the linear relation 

C',k(^') = E/ -2i'^n'K,nkM,l)fn:iM), (9-14) 

which defines the translation matrix Un'ki,nkz{d,'j)- The basis have a normal distance 
d and their axes are tilted against each other by an angle 7 as shown in Eq. (9.8). 



9.2 Translation matrices between non- parallel cylindrical multipole 
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The mathematical details of the derivation of this translation matrix are given in the 
Appendix A 9.6, and the translation matrix defined in Eq. (9.14) is obtained as 



(9.15) 



where ^ = ^, = ^, = H^^sM*^ ^nd k' = hi^^^^f^, 

^ p ' ^ p' ' sin(7) -L sin(7) 

For the calculation of Casimir energies, we need the translation matrices Ujj and 
Ujj, which are the translation matrices from object i to object j. For this particular 
problem, U12 = U„/fc/,„fc,(c?, 6*21, 7), and U21 = U„/fc/^,„fc^((i, vr + 6121, 7)- Notice that the 
sign of 7 is not modified. Therefore, using Eq. (9.15), it is easy to verify that 

U„.fc^,nfc.(rf,vr + ^2i,7) = (-l)"'""U„,fc.,„fc,(rf,02i,7)- (9.16) 



9.2.2 Electromagnetic waves 

For the EM field the translation matrix couples the polarizations. However, it can 
expected that, for large c?, the EM energy is given by the Dirichlet (TM) part only due 
to the logarithmic dependence on d for TM modes compared to the power-law decay 
of TE modes. However, the TM modes are coupled now by the translation matrix of 
vector waves, resulting in an attenuation of the asymptotic energy. The vector waves 
are defined in terms of the scalar cylinder waves as 

M^X = -V X , = -V X . (9.17) 

P i^P 

and equivalently for regular vector waves. To derive the transformation matrices, we 
first multiply both sides of Eq. (9.14) by z' and apply ^V'x to both sides. Then 
the left hand side gives the components of M°"^^ in the coordinate frame {x',y',z'). 
To obtain the components in the frame (x, y, z) we multiply both sides by the inverse 
rotation matrix where M^^ is defined by x' = M^^x + d (see Eq. (9.8)). This 
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yields, using V 



d,7 



V, 



MS7(x') = M,^M':S7(x') 



d,7 



dkz 
2^ 



U„'fc^,„fc^(rf,6')V X 



^ sin(6') sin(7)x ^ 



^reg 



— cos{9) sin(7)y 
y + cos(7)z 



(9.18) 



where we have used the definition of z in the initial coordinate basis given in Eq. (9.8) 



Eq (9.18) shows we need to express V x (0^'^^^(x)x) and V 



X 



Lreg 



(x)y) in terms of 



M^'^^^(x) and N^'^^^(x). One can show (by expressing V, x and y in cyhnder coordi- 
nates) that 



Vx KX(x)x) 



-i k^^ 



(9.19) 



V X Kl(x)y) = 1 (M-\,,^(x) - m:z,,M - 1 {K\uM) + KZ,,M)) ■ 

(9.20) 

After applying this result into Eq. (9.18), we need to shift the index n in order to 
express the right hand side in terms of vector waves with the same index n. This can 
be done by using 

This yields after some elementary algebra to the transformation formula for vector 
waves, 



out 



reg 



7 



(x), 



(9.22) 



^ , 27r ^n'ki,nkS'^^'y) \ 

where we have defined the EM translation matrix U,f,ff „;j^((i, 9) in terms of the trans- 
lation matrix for the scalar field given in Eq. (9.15) as 



l^n'K,nkM^)='^n'k'^,nkM,7)^ 



cos(7) + sin(7)^^ 


sin(7)^Vl + e 


-sin(7)^Vl + e 


cos(7) + sin(7)^^ 



(9.23) 



9.3 Asymptotic Casimir energy 



177 



Here we have used N = -V x M and -V x N = — M to obtain the transformation 

K K 

formula for N. 

The transformation formula for the inverse coordinate transformation is given by 
Eq. (9.22) with Vn'k'„nk, replaced by (-l)"+"'U„/fc/ ,„fc^, see Eq. (9.16), then 

K'KMM7r + 9,,,j) = (-l)"'-"Uf,tt„,^(t/,%,,7). (9.24) 



9.3 Asymptotic Casimir energy 

In this Section we obtain the approximation at large distances of electromagnetic and 
scalar Casimir energies at zero and high temperature limits. The approximation per- 
formed in this Section is valid for small cylinder radius or for large distance between 
the cylinders. 

To perform this limit we first apply the change to adimensional frequencies Qz = dkz, 
and q = d k, ioi the zero temperature case. As a consequence, U matrices depend 
on adimensional variables, while the Tj matrix of each cylinder scales with the ratio 
between the radius of this cylinder and the distance between their centers = In 
this Chapter, we assume that both cylinders have the same radius R, but the obtained 
results can be easily extended to the most general case without difficulty. 

After that, we apply the relation log \ A\ = Tr log{A) and perform the Taylor expan- 
sion to log(l — N) = — Yl'^=i ^'^^^''^ to Eqs. (9.2) and (9.3) in the adimensional variable 
r = ^, then it is an expansion on large distances, but also on small radius. Then the 
integrand is now log(l — N) ~ — Tr (N), but with the T matrices in N substituted by 
the tree level of their Taylor expansion over r. 



9.3.1 Scalar field 

For Dirichlet and Neumann boundary conditions, the scattering amplitudes of a cylin- 
der of radius R are given by the known expressions, 
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To obtain the asymptotic energy for Dirichlet boundary conditions cylinders at large 
distances d (small radius R), we need to consider only the terms with n = n' = 
and Pmax = 1 of N matrix of Eqs. (9.2) and (9.3). Then the main contribution 
at large distance of the T matrix with argument z is = — (7 + los(i)) ^ + 



O [z ] ~ — log^ (2;), with the Euler constant 7. For z = r^cp' + q^, we have ^^^^^ 
— log~^(r) up to logarithmic corrections. Then the Asymptotic Casimir energy at zero 
temperature between non parallel Dirichlet cylinders is 

^° = "8ci|sin(7)|log^(f)- ^^-^^^ 
For parallel cylinders of length L ^ (rf, i?), one has [22] 

^° = "8vrrf2 1og^(f)' ^^-'^^ 

so that one could identify | sin(7)| ^ tt-^. 

In order to study the classical limit we just need to consider the zero Matsubara 
frequency contribution to obtain the energy, but a simple scaling analysis shows that 
MoA;^,ofcz('^ = 0) ~ for /c^ — > with logarithmic corrections. Hence the trace 
of N is not well defined, i.e., ^Qk^fikX^ — 0) is not a trace class operator, so that 
|1 — ^Qk^fikX^ = 0)1 is not well defined. Kenneth and Klich showed that N is a trace 
class operator for compact objects [61]. While parallel cylinders are also non-compact 
objects, they constitute effectively a 2D problem with two compact discs. It appears 
that two tilted infinitely long cylinders provide the first example of non-compact objects 
whose geometry cannot be reduced to a lower dimensional one. 

However, we expect the force to be well defined, which means that the operator 
dt^Qk^idkX^ — 0) should be a trace class operator. This is indeed the the 
following calculation shows. We have that the force in the high T limit is 



4ci|sin(7)|log2(f)' 



(9.29) 
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and once we have the force, it is a straightforward exercise to obtain the energy as 
^ci = dxFci{x), then it is 

= TTTv- 9.30 

4|sin(7)|log(f) ^ ' 

Note that the sign of the energy is negative because the power of the logarithm has 

been reduced from two to one and R <^ d in this limit. This is the only case when 

we find a non trace class operator for scalar fields. When frequency increases or the 

boundary conditions of one of the cylinders is not of the Dirichlet type, all N operators 

are trace class operators, then they give finite contributions to finite energies. 

For example, for Neumann boundary conditions we have to consider the terms with 

Pmax = 1, |?^| < 1 and \n'\ < 1 since ^ftr) = + I'^l < 1- With this 



expansion we get to lowest order in ^ the energy 



d 



hcR'^ 



and the classical limit for this case is directly the finite result 
9.3.2 Electromagnetic Casimir energy 

For perfect metals, the scattering amplitudes of a cylinder of radius R are given by 
expressions 

^n'%^,nk^ = (pi?) ^"'""^^^^ ~ (9.33) 

where non diagonal terms are zero. We calculate the energies with the expansions of 
Eqs. (9.2) and (9.3). The leading part of the energy in the far distance approximation is 
given by the term that is quadratic in the T-matrix for the TM waves. The polarization 
mixing in the translation matrices given in Eq. (9.22) plays now a role, because the 
part of the polarization which is transformed from TM to TE does not contribute to 
the lowest order of Casimir energy, so the asymptotic energy must be lower than in the 
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scalar case, and equal to the scalar case just in these cases when polarizations do not 
mix. It is for parallel cylinders, as already known, and for the classical limit, because 
mixing terms are proportional to sin(7) and n. Then the asymptotic Casimir energy 
for perfect metal cylinders can be written as 

he 

8ci|sm(7)|log (f) 
where the function ^2(7) is defined as 

'"^'^ (cos(7) sin(y9) + sin(7) cos(v9) sin(0))^ 



^(7) = 7r~ d(psm{ip) . 2 2' 

^TT Jo Jo (sin(7) cos(v?) + cos(7) sin(v9) sin(0)) + sin (y?) cos2(0) 

(9.36) 

We have Q{0) = 1 so that for 7 — )■ the Dirichlet result given in Eq. (9.27) is 
recovered as expected. For orthogonal cylinders one has ^^(|) = 1 — log(2). We have 
not been able to perform the integral definition of ^(7). Anyway, a Fourier series 
expansion can be performed numerically giving us a good enough solution. Because 
the symmetry of the system, the series just contains cosines of even terms, then 
^(7) ~ X]^o ^2n cos(2n7). We find numerically that = 0.6137, ^2 = 0.3333, 
^4 = 0.0333 and = 0.0096. The function ^(7) is plotted in the Fig. 9.2. 

As in the classical limit the polarization mixing disappears, we have actually two 
scalar problems. One for TM modes, equivalent to a scalar field between cylinders 
subject to Neumann boundary conditions because of the perfect diamagneticity, which 
leads to a well defined trace class operator, and another one for TE modes, equivalent to 
another scalar field between cylinders subject to Dirichlet boundary conditions because 
of the perfect refiexibility. This part of the system is the dominant one at large distances 
and presents the same problem in the classical limit as its scalar analog. 

When we try to calculate the classical limit of the energy just considering the 
zero Matsubara frequency contribution, we find again the same problem as in the 
scalar field with Dirichlet boundary conditions. A simple scaling analysis shows that 
f^ok^,ok:,{i^ = 0) ~ 1^ for kz ^ with logarithmic corrections also for the leading 
order of classical limit of the energy for perfect metal cylinders. Also in this case 
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1.4 7 



Figure 9.2: Function $^(7) given in Eq. (9.36), as a function of the tilted angle 7. 



^ok^fik^i^ = 0) is not a trace class operator, so that |1 — Nok^.fik^'^ = 0)1 is not well 
defined [61]. It is another more physical case when non compact objects lead to a non 
trace class N operator. 

However, as in the scalar case, we are able to obtain the classical limit of the force, 
because ddNok^fik^i^ = 0) is a trace class operator as shown before. Then the far 
distance approximation of the classical limit of the force is 



" 4rf|sin(7)|log2(f)' 
and the corresponding energy is obtained with the integral Ed = dx Fci{x) as 



(9.37) 



(9.38) 



4|sin(7)|log(f)' 

Then we obtain the same result as in Eq. (9.30) of the scalar CclSG, clS expected. Once 
again, this anomalous behavior of the high temperature limit just can appear when the 
leading contribution of cylinder T matrices of both cylinders presents a log-divergence 
when the frequency k goes to zero. As in the scalar case, if just one cylinder presents 
this divergence, N matrix is still of trace class. In particular, for two dielectric cylinders 
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(or even for the case of a metallic cylinders and a dielectric cylinder) this anomalous 
behavior is not present. 

9.4 Numerical results 

Following the analytical asymptotic results of Casimir energies, we perform a numerical 
evaluation of Casimir energies between non parallel cylinders in order to study inter- 
mediate regimes and the convergence to asymptotic results. We center our study in 
the zero temperature case, then we perform a numerical implementation of Eq. (9.2). 
Finite temperatures could be performed on the same way, just being careful with the 
zero Matsubara frequency contribution for two perfect metals cylinders. This proce- 
dure is simple, but with some subtleties. If we observe the expression of the N matrix, 
we see that there is an internal integration loop over k± in order to obtain the matrix 
in {nkz, n'k'^) labels. After that we have to perform the logarithm of a linear operator 
over one set of discrete variables (n, n') and over another set of continuous variables 
{kz,k'^), because now the N matrix is not diagonal in kz- Obviously, the number of 
elements of N must be finite, which means we will use cylindrical multipoles with label 
n G [— "n^max; ''^max]- The cffcct of this cutoff is the transformation of the exact formula 
given by Eq. (9.2) into an asymptotic result, valid at a given error until a characteristic 
distance between the cylinders which depends on the rimax used as a numerical cutoff 
to obtain a numerical tractable method. In addition to that, we have to calculate the 
determinant of an operator defined over the continuous variables [k^, k'^). One possibil- 
ity would be the approximation of the determinant by a series of traces of powers of N, 
but this procedure introduces another new asymptotic approximation to the problem. 
The other possibility is to apply a discretization to the {kz,k'^) space and evaluate 
the determinant as usual with enough discretization points in order to minimize the 
discretization error. We expect that this procedure would be some kind of trapezium 
rule of integration for determinants of continuous operators. It is important to note 
that the evaluation of each N matrix is now several orders of magnitude harder than 
the parallel cylinders case. The final reason is the symmetry breaking of the system. 
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which not only produces that the N matrix is not diagonal in fc^, making the evaluation 
of the determinant harder, it also means that an internal integration over k± must be 
performed numerically, because it depends on the T matrix of one of the cylinders, 
while in the parallel case this integration is divided into two integrals independent of 
the T matrix and which can been performed analytically (see Appendix B 9.7). As a 
consequence, we just have obtained numerical results for perpendicular cylinders 7 = f 
and for 7 = | for all distances r = ^ < 0.45. Then we are far of the range of validity 
of PEA. We have used cylindrical multipoles just until \n\ < 3. It is enough to observe 
corrections to the small radius approximation and a soft convergence to PFA results. 

In Fig. 9.3 and in Fig. 9.4 we show the energies scaled with the PFA energy for 
cylinders with 7 = f and with 7 = | respectively. In the two figures we observe 
common characteristics with the parallel cylinders case. It looks like the energies (red 
curves) do not converge to their asymptotic result (blue curves) at large distances, 
i.e. small r, but it is not the case. As shown in Eq. (9.35), the asymptotic distance 
Casimir energy Eass decays with log~^(r), then the convergence is reached at extremely 
small distances and it is not observable in out plots. At intermediate distances, higher 
order cylindrical contributions get to produce relevant contributions to Casimir energy, 
and the energy became to be an appreciable fraction of the PFA result (black curve), 
until the approximations of the numerical method become relevant and we loose the 
convergence to PFA at shorter distances. The use of cylindrical multipoles with \n\ < 
nmax is valid for a given error until a characteristic distance between the cylinders 
which depends on the n^ax used as a numerical cutoff to obtain a numerical tractable 
method. We have checked numerically that the results obtained at r = ^ < 0.3 are 
well described with rimax = 3. To study even shorter distances and the convergence to 
PFA result, we should include more multipoles into our study. 

We are also interested in the dependence of the energy with the tilted angle 7 
between the cylinders at a given distance. As seen in Eq. (9.35) and in Eq. (9.61), 
the asymptotic distance approximation and PFA energies are inversely proportional to 
I sin(7)|, but they have different angular dependence. To study the angular dependence 
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0.1 0.2 0.3 0.4 r 0.5 



Figure 9.3: Energies compared with PFA as a function of r = ^ for 7 = f • The black 
curve is the PFA energy itself. The blue curve is the asymptotic distance energy and the red 
points are the numerical Casimir energies. The red curve joins the points r < 0.3, where our 
numerical method converges properly to the correct result. A convergence to the PFA result 
is expected when r — ^. 

of the energy at intermediate distances, we represent in Fig. 9.5 how is modified the 
product ^(r, 7) = i?(r, 7)! sin(7)| for all angles compared with the same function at 
7 = |, so we eliminate the dependence with the distance and we are able to study 
the relative dependence of the energy with the angle. We can observe that, at large 
distances, the dependence of Ci;(r, 7) with 7 is the one of the asymptotic results, but 
when we reduce the distance between the cylinders, a;(r, 7) becomes more and more 
independent of the angle, showing that it tends to converge to the constant value of 
the PFA result. 

9.5 Conclusions 

In this Chapter we have expanded the study of Casimir energy to non parallel cylinders, 
for the scalar field and for the electromagnetic field. 

In order to study this system, we have obtained the translation matrices of outgoing 
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Figure 9.4: Energies compared with PFA as a function of r = ^ for 7 = f • The black 
curve is the PFA energy itself. The blue curve is the asymptotic distance energy and the red 
points are the numerical Casimir energies. The red curve joins the points r < 0.3, where our 
numerical method converges properly to the correct result. A convergence to the PFA result 
is expected when r — )• |. 

cylindrical multipoles in terms of regular cylindrical multipoles of a translated non 
parallel basis in Sect. 9.2 for scalar and vectorial cases. 

We have studied the asymptotic far distance regime in Sect. 9.3 and near contact 
regime in Appendix 9.8, at the zero and high temperature regimes. 

We have found that the zero Matsubara frequency case leads to a non trace class N 
operator for the scalar field when both cylinders impose Dirichlet boundary conditions 
and for the electromagnetic field when both cylinders are perfect metals. The reason 
is that, contrary to the case of parallel cylinders, non parallel cylinders system is not 
reducible to an effectively 2D system of compact objects, then the theorem of Kenneth 
and Klich [61] is not applicable to this system, and N matrix could be a trace-class 
operator or not. Anyway we can recover an energy towards the integration of the 
force, because ddN{K = 0) is a trace-class operator. We should remark that the zero 
Matsubara frequency contributes to Casimir energy at any non zero temperature, then 
this result is relevant not only in the high temperature limit. 
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Figure 9.5: Ratio 



where c^(r,7) = ^(r, 7)] sm(7)| for angles 7 G [fg>f]- The 
asymptotic distance result and the PFA are also plotted. Curves are smoothing from the 
asymptotic result to the constant value of the PFA result when the distance is reduced, as 
expected. 



We have studied numerically the zero temperature Casimir energy between per- 
fect metal cylinders for all distances and angles, showing how asymptotic large and 
short limits are connected at intermediate distances, and the variation of the angular 
dependence at a constant distance. 

We have to remark that all the study made for perfect metal cylinders could be 
easily repeated for dielectrics. In this case at large distances the energy decays with a 
power law behavior, the zero Matsubara frequency leads to a trace class operator and 
the energy is generally smaller than in the perfect metal case. 

The study of Casimir energy between cylinders as a function of their relative angle 
7 could be experimentally relevant. We could use the different scales of the energy 
when cylinders are parallel, when E\\ oc instead oc ^|^.^^^^| for the non parallel 
case. It turns from an irrelevant value for non parallel cylinders to a much greater value 
proportional to their length L, which could be experimentally studied. In addition to 
that, Casimir effect could be relevant in nematic order of elongated fibers, because it 
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flavors the parallel configuration of them. 



9.6 Appendix A: Derivation of translation matrices 
for scalar cylindrical multipoles between non- 
parallel system of coordinates 

In this Appendix, we provide the mathematical details of the derivation the represen- 
tation of Eq. (9.15). It is the representation of 0°^*^, (x'), outgoing scalar cylindrical 
multipoles of a given basis in terms of 0,^'^^^(x), regular cylindrical multipoles of another 
basis translated a distance d and rotated around the axis defined by the translation 
vector an angle 6', as shown by the coordinate transformation given in Eq. (9.8). Then 
the translation matrix U„/fc/.,„fe^((i, 7) is defined by the linear relation (9.14). We start 
from the 2D Fourier transform of the outgoing waves, 



00 



Jn9 _ / "''^x / "So^^ „-\n '^x+^l^^y \ c 



KMe'- ^ ^ ^M-,r (.39) 
First of all, we apply a change of variables from {k^, ky) to (/c||, kyj^) defined as 

k\\ \ I cos(9ji) sin(6',j) \ ( k^ \ 
k± J y -sm(%) cos{e ji) ) \ h ) 
and the same change of variables to the conjugate variables in the position space as 



cos( 


Oji) 


sin( 


Oji) 


— sin( 


Oji) 


cos( 


Oji) 




(9.41) 



therefore, the integral of Eq. (9.39) can be written in terms of the new integration 
variables as 

r°° rlh, r°°rlh. / h,i -\- ik , \ ^ „i{k»r»+kj_rj_) 



Kn{pp)e'^'= ^27r(-2)" ^ i e*"'^--^— 2 ^. (9.42) 

7-00 27r 27r ^ ^ \ p J k^ + kl+p^ 

The integration over k\\ can be easily performed using the residue theorem. The 



integrand has the poles fcy = -^i^p^ + k\. Hence we obtain 

= e-.. / I j (9.43) 
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where the minus (plus) sign apphes to ry > (r|| < 0). After multiplying by e*^^^, we 
obtain 

J-oo \ P J 

Writing the same formula in the prime coordinate system: 



/oo 
dk 
■oo 



(9.45) 



2 



and using p = a/kM-A^I, = a/k^ + A;'^, it is easy to obtain A;^ + A;^ = A;^^ + A;' 
A;_Lr_L + fc^z = k'j_r'j_ + k'^z' and k\\ = k'\\ = a/ k'^ + p'"^ by the use of Eq. (9.11). Then, 
applying the change of coordinates given in Eq. (9.8), we have the outgoing cylindrical 
functions of one basis in terms of outgoing plane waves of the translated and rotated 
basis as 

KAp'py'^'e^'^^' = e-'^^' / dk'^ k^ + Vk^+P e V 

(9.46) 

Thanks to the relationship {xk^. + t/A;^^) = (r||A;|| + r_|_A;_|_), we can use the expansion of 
2D plane waves in cylindrical waves as. 



gi(r||fc||+rxfcx) 



^i{xk^+yky) ^ J^i'^Jn {p\Jkl + kl^ e^"^e Vy^/^ (9.47) 



where p = ^Jx^ + y'^. Obviously we have cos ^ ( , ^ ) = sin ^ ( , ^ ) . Apply- 
ing the relation p = i^yk^ + A;^, we transform the equality to 



J2 ^"^n (^PP) e^"^e-™^"' U; , (9.48) 

neZ 

which can be written as 



(9.49) 
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Using the relation p = i^k'^ + A;^, it is easy to obtain that a/^^ + = "^^x- Therefore, 
we can write {ky + a/^J|~+~P^) = {ky — ik^) = e''-^^'{k± — ik\\) = e^^^'{k± + a/pM-AJJ) and 
the series acquire an additional phase proportional to 9ji 

^Hrfy+r^k^) ^ (pp) e'ne^-ine,^ fk^±VK±Z\ . (9.50) 

nez V ^ / 

By the insertion of this sum in the integrand, we are able to bring the result as a sum 
of modes and a integration over k'j_ frequencies as 



"n. / I — \ —n 

p2 



X 



2^k'\+p''^ 

In order to bring the result into the form of Eq. (9.14) we change the variable of 
integration from fc_L to k^ with the help of the relation k^ = cos{'~f)k'^ + sin(7)fc^, then 
^^'± ~ [slr^7)l' expansion of outgoing cylindrical function in terms of regular 

waves of a non-parallel cylindrical basis is 



^ 27r L |sm(7)| \^ p' 



2^ " 



X 



where 



X 1 ^^ + ^^^^^ ) '~'f^, Inipp)e-'e^'^^ (9.52) 
P J 2^/k'l+p^ 



k, = !^^^!M^, k', = (9.53) 

sm(7) sm(7) 



Then the translation matrix defined in Eq. (9.14) is: 



^ |sm(7)| V ^ y V ^ / 2v/A;'l+p'2 

with ^ = ^ and = which is the result given in Eq. (9.15). It is interesting to note 
that an additional translation over the z axis, this is to say, for cases with dz 7^ 0, the 
generalization of the translation matrices showed here is straightforward and results in 
just an additional phase factor e"*'^"'^"' (which appears in the coordinate transformation 
of z' showed in Eq. (9.46)) in the expression of Un'k'^,nk^{d,'j) given in Eq. (9.54). 
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9.7 Appendix B: Translation matrices for scalar cy- 
lindrical multipoles between parallel system of 
coordinates 

In this Appendix, we reobtain the translation matrix between parallel cylindrical multi- 
poles as a limiting case when 7 — )■ of Eq. (9.15) for scalar and Eq. (9.22) for vectorial 
multipoles. 

In the parallel case, the cross polarization matrix which appears in Eq. (9.22) is 
reduced to the identity matrix. Then for the parallel case polarizations do not couple. 
As a consequence, Eq. (9.22) is reduced to two scalar decoupled translation problems, 
one for each independent polarization. So it is enough to obtain the translation matrix 
for scalar multipoles in order to obtain the result for the vectorial case. For the parallel 
case we have kj_ = k'j_ and = k'^. As a consequence, ^ = 1 and ^' = ^. Then we 
can follow the derivation of the scalar translation matrix until Eq. (9.51), where for 
the parallel case we have 

^n'k',,nk.id^1 = 0) = 

, , , 27r / , \"-"' p-'^a/«^+^I+^ 

(-ire-('-*y_^*,_.(fe-.;,(^ + /rTF) (9^55) 

Instead of applying the change of variable from k± to k^ , which would be singular in the 
parallel case, we carry out the integration over k±, resulting in the known translation 
matrix between parallel cylindrical scalar multipoles and for each polarization of vector 
multipoles as 

Un'ki,nkMl = 0) = {-ir'e'^^'-'''^'^^Kn-n'{d./^^+kl)—6{k, - k',). (9.56) 
Then we recover the parallel translation matrices [22] as a particular case. 
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9.8 Appendix C: Proximity Force Approximation 
of Casimir energy 

Small Radius Approximation is valid for large distance. In the short distance regime 
Proximity Force Approximation (PFA) is the valid one. We just perform the calculus 
for perfect metal cylinders. The scalar case is just a half of the perfect metal result. 

The area across which the two cylinders overlap, viewed along the axis that is 
perpendicular to the two cylinder axes and that intersects the axes in their crossing 
point, forms a parallelogram of edge length 2R/\ sin(7)|. Let us denote the coordinates 
along the edges of this parallelogram as u and v (see Fig. 9.6). Then the local distance 
h{u, v) between the two cylinder surfaces, measured normal to the plane that is spanned 
by the cylinder axes, is given by the function 

h{u, v)=d-\lR'^-{u\ sin(7) \ - Rf - sj ~ {v\ sin(7) | - Rf . (9.57) 

where d is the distance between the cylinder axes. Taking into account that a surface 
element of the parallelogram is given by si\i['^)dudv , the PFA energy can be written as 



2R 



This expression has the advantage that we can expand the square roots for small l/R, 
where I = d — 2R is the surface-to-surface distance, which leads to 



E^^ = - , \ „ g / / dsdt 
' 720 |sin(7)|/2 7_ 



R/l 



(9.59) 



The energy can be written after introducing new integration variables as 



hcTT^ 1 R /-V^ /-V^ dsdt 



In this expression we can extend the integration limits to infinity to obtain the limiting 
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Figure 9.6: Region of integration for PFA. The unit vectors u and v are shown too. 



behavior for small l/R. This yields 

1^0° 720 |sin(7)|(2/, ^ J, (i + dV 720|sin(7)|P ^ ' 

For small l/R this approximation deviates from the exact integral of Eq. (9.59) by less 
than 1% for l/R < 0.01. It is instructive to compare this PFA energy to the one for 
two spheres of radius R and surface-to-surface distance / which is E^^^ = — fH^^ [21] 
and hence reduced by a factor of two compared to the case of perpendicular cylinders 
(7 = 1). 

With exactly the same procedure, it is possible to obtain the PFA energy for the 
classical limit, performing the same integration, but now over the classical limit of the 
energy instead of the zero temperature case. Then the classical limit of the PFA energy 
is 

^PPA TnC(3) f f |sin(7)l | SlU (7) | ^Mfif 



= -ksT^-p- / / ' , : . (9.62) 
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Performing the same transformations and approximation as in the quantum case, we 
find that in the hmiting behavior for small Z/i?, the PFA energy is 

i,^o ^ ^ 87r|sm(7)|Uo 7o (l + f ) ^ |sin(7)M 

(9.63) 
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'So Long, and Thanks for All the Fish. " 

Douglas Adams 




Conclusion 



This Thesis is focused on the study of Casimir effect as a consequence of fluctuations 
of fields and its peculiarities and properties. 

We have obtained results in several categories. 

• On the conceptual side, we have proposed a dynamical formalism of Casimir 
effect, which let us define the Casimir effect as the response of a fluctuanting 
medium to the breakdown of the translation symmetry because of the presence of 
intrusions in that medium. Such formalism let us focus on the Casimir effect as 
an universal feature of steady states. We have defined the stochastic variable 
" Stochastic Casimir force over an a body" whose mean value is the Casimir force 
itself. 

• The formalism presented here generalizes the Casimir effect to non-equilibrium 
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steady states systems, and it is demonstrated that the equihbrium Casimir effect 
is a particular but important case. In fact, we obtained the exphcit equivalence 
between the formalism presented here, the Stress-Tensor formalism and the par- 
tition function formalism of equilibrium Casimir effect. Even the important case 
of Casimir effect generated from quantum-thermal fluctuations of the EM field 
is covered as a particular case with this formalism, via the Parisi-Wu formalism. 

• As the Casimir force is a stochastic variable, it must follows some probability 
distribution function, then it has interest the study of its moments, in particular, 
we defined the variance of the Casimir force, which let us study the fluctuations 
of these Fluctuation-Induced forces. In fact, we have obtained the variance of 
the Casimir effect between two perfect metal pistons of arbitrary cross section, 
at any given temperature at the first time to our knowledge. 

• On the mathematical side, we have obtained the two point correlation function 
of steady states generated by a kind of generalized Langevin equation when the 
temporal evolution operator is generalized to a polynomial of any degree. We 
obtained the Casimir force between two parallel Dirichlet plates in such steady 
states. 

• On the applied side, we used the recent developed version of the multiscattering 
formalism to study the physical peculiarities of the EM Casimir effect: Nonmono- 
tonicities of Casimir forces between two objects because of the presence of a third 
one, or because of the appearance of intervals of negative entropies in the sys- 
tem; the obtention of the PSA from first principles for the complete EM field, the 
generalization of the N points formalism of Power and Thirunamachandran, the 
Casimir energy between non-parallel cylinders, and the Casimir effect between 
Topological Insulators with their anomalous behavior. 
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10.1 Specific outcomes 

• Presentation of the dynamical formalism of Casimir effect. The Casimir effect is 
the response of a fluctuanting medium to the breakdown of the translation sym- 
metry because of the presence of intrusions in that medium. 

1. Definition of the Stochastical variable "Stochastic Casimir force over the 
object a" 

2. Generalization of Casimir effect to non-equilibrium steady states via the 
dynamical formalism. 

3. Reobtention of the equilibrium Casimir effect as a particular case of the 
dynamical formalism. 

4. Explicit equivalence of the Stress-Tensor formalism, the partition function 
formalism and the dynamical formalism of the Casimir effect in the equilib- 
rium. 

5. Reobtention of the originally proposed EM Casimir effect as a consequence 
of the quantum fluctuations of the EM field via Parisi-Wu formalism. 

6. Definition of the variance of the Casimir force, and result for perfect metal 
pistons of arbitrary section. 

7. Generalization of the Langevin equation to temporal derivatives greater than 
one, and obtention of Casimir force in such steady states. 

• Derivation of the Pairwise Summation Approximation (PSA) for the EM Casimir 
effect from first principles from the multiscattering formalism of Casimir effect, 
including next order perturbation terms. 

1. New results of PSA for dielectrics with electric response, magnetic response 
and general magnetoelcctric coupled responses. We include indications on 
how to proceed when the electromagnetic response of the objects are not 
homogeneous. Results valid for any given temperature. Recovery of the su- 
perposition principle in the PSA, and breaking of it with the first correction 
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term to PSA. 

2. Reobtention of the N points potential formalism of Power and Thiruna- 
machandran as the large distance limit of PSA, and extension of such for- 
malism to generalized (linear) dielectrics at any given temperature. 

3. Use of PSA for the complete EM field for the study of the Casimir effect 
between Topological Insulators. Evaluation of the adequacy of the approx- 
imation as a good qualitative and quantitative approximation. Criteria for 
obtaining the appearance of the different regimes of behaviors of TI depend- 
ing on the magnitude and sign of the magnetoelectric polarizabilities. 

• Nonmonotonicities of the Casimir force between two objects. 

1. Nonmonotonicities of the Casimir force between two objects because of the 
presence of a third one. 

2. Study of the entropy of a system of two metallic spheres originated because 
of the Casimir effect. Appearance of intervals of negative entropy at an 
interval of given temperatures and distances for Perfect metal model and 
plasma model with lower penetration length of a critical one. No appearance 
of this range for all other cases of the plasma model and for the Drude model. 

3. Thermodynamic consequences of the appearance of negative entropies of 
Casimir effect. The principles of Thermodynamics continue to be satisfied, 
of course, but the Casimir force is not monotonous with temperature pro- 
vided that the entropy is non-monotonic with distance. 

• Casimir energy between non-parallel cylinders. 

1. The energy scales with the cosecant of the angle between their axis, nor with 
their length L. 

2. New translation matrices between translated non-parallel cylindrical coor- 
dinate systems. 



10.1 Specific outcomes 



199 



3. First counterexample of the theorem of Kenneth and Khch [61] in the high 
T hmit for Dirichlet cyhnders and for perfect metal cylinders, because the 
condition of compactness of the intruders is not fulfilled. 
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